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ABSTRACT 


This research was directed towards an experimental investigation of 
the performance of the immiscible carbon dioxide Water-Alternating-Gas 
(WAG) process, which holds considerable promise for the recovery of moder- 
ately viscous heavy oils from thin and/or deep formations, for which reser- 
voir conditions do not favour the application of any of the thermal recovery 
techniques. 


The experimental work was conducted to study the diffusion and grav- 
ity segregation phenomena occurring during the process. Work was also 
done to study the non-isothermal immiscible carbon dioxide WAG process. 


A number of measurements were made to study the diffusivity of a gas 
phase into a liquid phase. In this study, two gases: carbon dioxide and 
methane were used as the diffusing gas phase. Oils of different physical prop- 
erties were used to represent the liquid phase. The measurements were made 
at various pressure and temperature conditions. It was observed that, based 
on the results obtained, the diffusivity of gaseous carbon dioxide or methane 
increased with increasing pressure and temperature and decreased with in- 
creasing oil viscosity and oil molecular weight and that carbon dioxide dif- 
fused faster than methane. Using the data collected, an empirical correlation 
was also developed. 


Several displacement experiments were performed to determine the 
possible application of a carbonated waterflood in place of an immiscible car- 
bon dioxide WAG flood. A carbonated waterflood was found to be inferior to 
an immiscible carbon dioxide WAG flood, because the carbon dioxide re- 
quirement was too high. 


Gravity segregation of carbon dioxide and water was investigated by 
conducting a number of vertical and horizontal displacement experiments. It 
was observed that gravity segregation affected the displacement efficiency of 
the immiscible WAG process. As well, it was noted that transverse diffusion 
of carbon dioxide in the horizontal direction normal to the vertical longitudi- 
nal direction helped delay the upward gravity channelling of carbon dioxide. 
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In addition to the studies mentioned above, a non-equilibrium math- 
ematical model with phase change for the non-isothermal immiscible carbon 
dioxide WAG process was developed, using non-equilibrium thermodynam- 
ics theory. Two sets of scaling criteria for the non-isothermal carbon dioxide 
process were derived using the mathematical model. These were employed 
to design and perform non-isothermal experiments. 


Moreover, the experimental results obtained in this investigation, as 
well as those obtained in previous investigations, were correlated using the 
scaling criteria derived, thus demonstrating the usefulness of the criteria. 
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NOMENCLATURE 


area open for flow in an injection or production well, m2 
scaling factor 

interfacial area between phase j and 1, m* 

concentration of the i component in the j phase, mass fraction 
molecular diffusion coefficient, m?/s 

longitudinal hydrodynamic dispersion coefficient, m*/s 
transverse hydrodynamic dispersion coefficient, m?/s 
convective dispersion tensor of the j phase, m*/s 


molecular diffusion tensor of the i component in the j phase, m*/s 
hydrodynamic dispersion tensor of the i component in the j phase, 
m?/s 

rate of total energy of the j phase injected at the injection well, kJ/s 
formation electrical resistivity factor 

buoyancy force exerted on the j phase per unit volume, N/m? 
drag force exerted on the j phase per unit volume, N/m?® 

fugacity of the i component in the j phase, kPa 

gravitational acceleration, m/ 52 

height or thickness of the reservoir or model, m 

enthalpy of the j phase per unit mass, kJ/kg 

local mass transfer coefficient of the i component in the j phase 
related to the jl interface, (m*-s)"! 

effective permeability tensor of the j phase, m? 

relative permeability of the j phase, dimensionless 

thermal conductivity of the j phase, kW/m-K 

absolute permeability, m? 

Boltzmann’s constant (1.38E-23 J/K) 

length of the reservoir or model, m 

rate of momentum of the j phase at the injection well, kg-m/s* 
rate of water evaporation or condensation per unit volume, 
kg/m*-s 

mass transfer rate of the i component in or out the j phase through 


jl interface per unit volume, kg/m°-s 
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P. capillary pressure, kPa 


Pj pressure of the j phase, kPa 

@) volumetric injection rate, m?/s 

R radius of molecule, cm 

Sj saturation of the j phase, fraction 

8; entropy of the j phase per unit mass per absolute temperature, 


kJ/kg-K; or rate of entropy created at the injection well, kJ/s-K 


T temperature, °C or K 

‘ time, day or second 

Uj internal energy of the j phase per unit mass, kJ/kg 
9; superficial velocity vector of the j phase, m/s 
9; superficial velocity of the j phase, m/s 

V total volume of fluid injected, fraction 

Vco, slug volume of carbon dioxide, fraction 

V, slug volume of water, fraction 

W width of the reservoir or model, m 

W; mass injection rate of the i component, kg/s 
Ve cartesian coordinates 


Greek Symbols 


i) porosity, fraction 

Pj density of the j phase, kJ/kg 

T; shear stress of the j phase, N/m? 

Ne chemical potential of the i component in the j phase, kJ/kg 
Lj dynamic viscosity of the j phase, mPa.s 

or interfacial tension between the j and | phases, kg/s? 
Subscripts 


A solute 

B solvent 

D dimensionless quantity 

g gas phase 

CO2,g CO2 component in the gas phase 
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CO2z,0 | CQOz2 component in the oil phase 
CO2z,w CO? component in the water phase 


inj located at the injection well 

¢ liquid 

M model 

O oil phase 

0,0 oil component in the oil phase 
Ne prototype 

prod located at the production well 
R reference quantity 

r rock matrix 

e vapour 

Ww water 

wW,w water component in the water phase 
x x-direction 

y y-direction 

ve z-direction 


Superscripts 

I located at the interface 

Ijl located at the interface between the j and | phases 
Abbreviations 

Abs. absolute 

BD blowdown 

BT breakthrough 

d darcy 


GOR gas-oil ratio 

HCPV hydrocarbon pore volume 
Inj. injected 

LC linear core 

MW, __ oil molecular weight 

Rec. recovery 
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two-dimensional 
velocity 

viscosity 
water-alternating-gas 
waterflood 
water-oil-ratio 
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INTRODUCTION 


The immiscible carbon dioxide WAG (Water-Alternating-Gas) process 
holds promise for the recovery of moderately viscous heavy oils from thin 
(less than 10 m) and/or deep (greater than 1000 m) formations where thermal 
recovery methods are ineffective due to high heat loss and other limitations. 


Many studies have been conducted to date to examine the effectiveness 
of the process, as well as its mechanism. Laboratory studies on the application 
of the process to the recovery of oil from thin and deep reservoirs showed 
that a substantial volume of oil, about 10 to 30% incremental oil compared to 
a waterflood, can be recovered by this process!*. Field studies and projects 
conducted in the United States have produced good results?!9. In Alberta, the 
effectiveness of this process was demonstrated by the commercial immiscible 
carbon dioxide project in Retlaw, started in January 1991. Within six months 
after production began, the cumulative oil production exceeded 120,000 sm?.!! 


The basic idea of this recovery process is to alternately inject small slugs 
of carbon dioxide and water. When carbon dioxide is injected into a reser- 
voir, it will dissolve in the oil, causing swelling of the oil and lowering of the 
oil viscosity. Following this, water is injected to displace the carbon dioxide- 
swollen oil. The reason slugs of carbon dioxide and water are alternatively 
injected is to control the mobility of the gas phase. Carbon dioxide gas is pre- 
ferred in this process because it has a very high solubility in oil, compared to 
other gases such as methane, propane, ethane, and nitrogen. In order to 
make the process more effective, carbon dioxide must be injected as gas rather 
than liquid because liquid carbon dioxide may only be sparingly soluble in oil. 


Even though the process is quite successful in recovering moderately 
viscous heavy oils from thin and deep reservoirs, there are still a few prob- 
lems associated with it. Two of the few problems existing in the process, 
which are also the focus of this current study, are non-equilibrium phenom- 
ena and gravity segregation of the injected carbon dioxide gas. The non-equi- 
librium phenomena arise because the concentration of carbon dioxide near or 
around the injection well is higher than that away from or near the produc- 
tion well. Consequently, mass transfer has to take place to even out the car- 
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bon dioxide concentration. One of the well known mass transfer processes is 
diffusion which occurs when a system is out of equilibrium. Diffusion is the 
mass transport process which involves the movement of molecules from one 
point to another and is known to play an important role in the immiscible 
carbon dioxide WAG process. Therefore, one of the objectives of this research 
is to measure the diffusivity of carbon dioxide in different heavy oils at vari- 
ous pressures and temperatures. A second objective is to develop a mathe- 
matical model which includes the non-equilibrium phenomena. 


As mentioned above, because carbon dioxide is injected as gas rather 
than as liquid, its density is much lower than that of the reservoir oil. The 
density difference between the injected carbon dioxide and the reservoir oil 
results in segregation of the injected carbon dioxide, causing carbon dioxide to 
rise to the top of the reservoir and then to finger through the oil. 
Furthermore, the solution of carbon dioxide in oil increases the oil density, 
making the density difference worse. 


Studies conducted to-date have been focussed on investigating gravity 
segregation in miscible gas displacements; none has been conducted to inves- 
tigate the same phenomenon in immiscible gas displacements. Thus, there is 
a need to study gravity segregation in the immiscible carbon dioxide dis- 
placement process. 
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STATEMENT OF THE PROBLEM 


Gravity segregation of the injected fluids and non-equilibrium phe- 
nomenon may have detrimental effects on the performance of the immisci- 
ble carbon dioxide WAG process. In this study, the effects of gravity segrega- 
tion and non-equilibrium phenomena in the immiscible carbon dioxide 
WAG flooding process are to be investigated. The main objectives of this 
study can be summarized as follows: 


1. To study the diffusivity of carbon dioxide into oil at various 
pressures and temperatures. 


2. To design a mathematical model for the non-equilibrium and non- 
isothermal immiscible carbon dioxide WAG process, which is used to 
derive scaling groups. 


3. To correlate the experimental results using the scaling groups de- 
rived in this study. 


4. To investigate the effect of gravity segregation on the immiscible car- 
bon dioxide WAG process and to determine whether gravity segrega- 
tion increases or decreases the efficiency of the process, thus observing 
how it affects oil recovery from thin reservoirs. 
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REVIEW OF THE LITERATURE 


The idea of using carbon dioxide as an enhanced oil recovery agent is 
not new. As early as in the 1950’s, Martin!?"!4 conducted the first study on the 
possibility of flooding an oil reservoir with carbonated water rather than with 
plain water. He found that carbon dioxide has a strong tendency to improve 
the mobility of the oil by lowering the viscosity of the oil after it dissolved in 
the oil. Following the pioneering work of Martin, other investigators!>-?? 
also recognized that viscosity reduction due to the solution of carbon dioxide 
in oil led to improved sweep efficiency, thus increasing oil recovery. Holm 
and Josendal*? experimentally investigated the miscible displacement of a 
very light oil by carbon dioxide and discovered that using carbon dioxide to 
displace oil under miscible conditions could recover up to 95% of the in-place 
oil. However, at that time, due to relatively high primary oil production and 
the very high cost of carbon dioxide compared to the price of oil, the use of 
carbon dioxide in enhanced oil recovery was not attractive to most oil pro- 
ducers. Today, due to low primary oil production, high oil demands, and 
plentiful sources of carbon dioxide available at relatively low prices around 
the globe, there is an increasing emphasis on using carbon dioxide as an en- 
hanced oil recovery agent. 


The objective of this chapter is to provide a review of past work, to 
study the displacement of oil by carbon dioxide, and its relationship to the 
current study, which was directed towards an investigation of the immiscible 
displacement of moderately heavy oils by carbon dioxide and water. 


3.1 - Diffusivity of Carbon Dioxide into Oil at Reservoir Conditions 


In the immiscible carbon dioxide process, the four mechanisms con- 
tributing to oil recovery, according to Rojas and Farouq Ali’s findings*4, are 
oil swelling, viscosity reduction, interfacial tension reduction leading to the 
formation of water-in-oil emulsions, and solution gas drive resulting from 
blowdown at the end of the displacement. The so-called oil swelling occurs 
when carbon dioxide dissolves in oil. There are two mass transfer processes 
which take place when carbon dioxide dissolves in oil: solution and diffusion. 
It is known that solution is a fast process whereas diffusion is a slow process. 
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In the current research, only the diffusion process is studied. Hence, the next 
section provides a review of the work done by other researchers on the 
diffusivity of a gas phase in a liquid phase. 


3.1.1 - Measurement of Diffusion Coefficient 


As mentioned in the above paragraph, diffusion is another means of 
mass transfer between carbon dioxide and oil, besides solution. It involves 
the movement of molecules from one point to another due to concentration, 
temperature, and pressure gradients; and it is independent of any convective 
forces in the system. It can also be defined as dispersion at the molecular 
level. Dispersion is the mixing of fluids which occurs when one fluid dis- 
places another. 


According to Grogan and Pinczewski”’, molecular diffusion plays an 
important role in the recovery of residual oil at the pore level since it is an 
important rate controlling mechanism in a carbon dioxide flood. Several 
authors*4-?7 have pointed out the advantages of diffusion. It helps carbon 
dioxide to penetrate oil, inhibit viscous fingering, delay gas breakthrough and 
increase the oil rate. 


Since 1933, efforts have been made by several investigators to measure 
experimentally molecular diffusion coefficients at elevated pressure and tem- 
perature conditions. Pomeroy, Lacey, Scudder, and Stapp?® were the first to 
design an unsteady-state method which allowed the determination of diffu- 
sion constants at pressure up to 3 MPa. The method was based on Fick’s sec- 
ond law of diffusion. They reported values in the range of 5x10”? to 14x10? 
m?/s for the diffusivity of methane in isopentane for pressures in the range 
of 0 to 2 MPa and at 30°C. Employing the same experimental apparatus, Hill 
and Lacey”? measured the diffusivity of methane in different oils for tempera- 
tures in the range of 30°C to 60°C and at 2 MPa. The reported diffusion coeffi- 
cients are in the range of 3.8x107!° to 1.3x10°° m*/s. Diffusion was noted to 
increase with increasing temperature. Hill and Lacey*° also presented a set of 
data for the diffusivity of propane in liquid hydrocarbons for pressures in the 
range of 0.2 to 1.4 MPa and temperatures from 30 to 60°C. They reported val- 
ues in the range of 1.7x10° to 6.3x10° m*/s. The diffusivity of propane was 
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found to increase with increasing temperature and pressure. Other investiga- 
tors®!-5 have also reported values for the diffusivity of methane and propane 
in various liquid hydrocarbons over various pressure and temperature 
ranges. 


There have been many measurements of the diffusivity of methane in 
hydrocarbon solvents at elevated temperatures and pressures, but very few 
for the diffusivity of carbon dioxide. The only measurements reported are 
those of Schmidt, Leshchyshyn, and Puttagunta®®; Denoyelle and Bardon?’; de 
Boer, Wellington and Tschiedel®®; Renner®?; and Grogan, Pinczewski, 
Ruskauff, and Orr*®. Schmidt et al.°° performed the first work on determin- 
ing the diffusivity in reservoir fluids at reservoir conditions. They measured 
the diffusivity of carbon dioxide in a 362,000 mPa.s bitumen at 5 MPa and 
temperatures ranging from 20 to 200°C. Their results demonstrated that the 
diffusivity of carbon dioxide in bitumen increased with increasing tempera- 
ture. Denoyelle and Bardon*’ reported that the diffusivities of carbon dioxide 
in oil at reservoir conditions were about 5 to 10 times higher than those mea- 
sured at atmospheric conditions. They concluded that the diffusivity coeffi- 
cients at atmospheric conditions could not be used to estimate those at reser- 
voir conditions. de Boer et al.°® studied the diffusivity of carbon dioxide 
through an interface (water layer) into an oil droplet. Their observations are 
at variance with the work of Denoyelle and Bardon®’. They observed that the 
diffusion rates of carbon dioxide in crude oil at elevated pressures were con- 
sistent with calculated rates based on diffusion coefficients measured at atmo- 
spheric conditions, provided that there was no precipitation of asphaltenes. 
They also observed that the precipitation of asphaltenes at the oil-water inter- 
face formed a highly resistive layer which retarded carbon dioxide mass trans- 
fer. 


Employing the method of Reamer et al.3*, Renner®? measured the dif- 
fusivity of carbon dioxide in decane and brine in consolidated porous media 
at pressures from 0.9 to 5.8 MPa and at 38°C. In this experimental study, he 
noted that the diffusivity of carbon dioxide increased with increasing pressure 
when the core was placed in the horizontal position and was independent of 
pressure when the core was in the vertical position. Grogan et al.*° per- 
formed their measurements for carbon dioxide diffusion in pentane, decane 
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and hexadecane. The reported data indicated that the diffusion coefficient de- 
creased with increasing solvent molecular weight and increased with increas- 
ing carbon dioxide concentration. Furthermore, they also concluded that cor- 
relations for the diffusivity of carbon dioxide in hydrocarbon solvents in 
terms of solvent viscosity, developed from measurements at atmospheric 
pressure, provide realistic estimates for the diffusion coefficient at reservoir 
conditions. This is opposite to Denoyelle and Bardon’s conclusion?’. 


3.1.2 - Correlations of Diffusion Coefficients 


In addition to experimentally determining the diffusion coefficient, ef- 
forts have been made to propose correlations for estimating the diffusivity of 
carbon dioxide. The first effort was made by McManamey and Wollen*!, who 
presented an empirical correlation for the diffusivity of carbon dioxide in or- 
ganic liquids as a function of solvent viscosity only. This correlation was later 
shown by Denoyelle and Bardon?’ to give poor estimates of carbon dioxide 
diffusivity at pressures higher than atmospheric pressure. By fitting the data 
reported by Schmidt et al.°° to the Nernst-Einstein equation (which is 


#2 ~ €KT - 4,49 Aa - : as 
Ta aRGEE, Schmidt** presented an empirical correlation for the diffusivity 


of carbon dioxide in bitumen as a function of viscosity and temperature. 
Renner?’ also correlated his data and others’ to come up with a correlation for 
predicting the diffusivity of carbon dioxide, methane, ethane, and propane in 
liquid hydrocarbons. His correlation showed that increasing pressure de- 
creased the diffusivity of a gas phase in oils. Mehrotra, Garg, and Svrcek*? 
used Schmidt et al.’s data*® to evaluate several correlations. They found that 
Umesi-Danner’s correlation, together with the use of the corresponding states 
theory of Teja*4, could be used to predict the diffusivity of carbon dioxide in 
bitumen at reservoir conditions. The disadvantage of using this technique 
for field applications is that Teja’s method is complicated; that is, it requires 
too many calculations and thermodynamics data. 


3.2 - Carbonated Water Injection 


The carbonated water injection process is quite different from the 
WAG injection process. Water is saturated with carbon dioxide at the surface 
before being injected into the reservoir. This process was first studied by 
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Martin!2"4 and Saxon, Breston, and Macfarlane* in the early 50’s, and re-in- 
vestigated by Holm?° and de Nevers*!/46 in the 60’s. Based on the results 
gathered, Saxon et al.*° concluded that a carbonated waterflood would not re- 
duce the oil saturation any further than a waterflood. Holm? reported 
higher carbon dioxide requirements and lower recoveries for a carbonated wa- 
terflood, as compared to a carbon dioxide slug flood. 


In the 50’s and 60’s this oil recovery technology was tried in oilfields in 
New York and Oklahoma*®48. After more than a year of operation, no signif- 
icant oil was recovered. It was tried again in Texas by Amoco in 1987. After 
only 7.5 months, injection was terminated because of extensive corrosion and 
plugging problems. There was no mention of any incremental recovery??. 


Perez, Poston, and Sharif? studied carbonated water imbibition in core 
plugs at various pressure and temperature conditions. They concluded that 
this method may hold promise for increasing oil production rates from frac- 
tured, low matrix permeability, and low gas-oil ratio oil reservoirs. 


3.3 - Gravity Segregation 


According to Morrow and Hornof”!, gravity segregation occurs when 
capillary and viscous forces are insufficient to overcome the effect of buoy- 
ancy forces. The first laboratory study of gravity segregation performed by 
Craig, Sanderlin, Moore, and Geffen*? in 1957, showed that in linear gas or 
water injection operations in horizontal formations, segregation of the fluids 
due to the gravity effect could result in oil recoveries to breakthrough as low 
as 20% of those otherwise expected and that, in 5-spot injection operations in 
such formations, the oil recoveries at breakthrough could be as low as 40% of 
those predicted by methods which assume negligible gravity effects. Their re- 
sults also indicated that the magnitude of segregation of the fluids due to 
gravity was influenced by the average injection rate rather than temporal 


variations. 


Richardson and Perkins’? studied experimentally the effect of rate on 
oil recovery by waterflooding. They concluded that decreasing the rate at 
which water is injected into a thick, homogeneous reservoir sand increases 
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the tendency of the water to underrun the oil. Miller°4 also observed that 
water underran the oil in his study. 


Blackwell, Terry, Rayne, Lindley, and Henderson”? reported an investi- 
gation of the efficiency of water-solvent mixtures for oil recovery. They 
observed that, while flowing through sands the water and gas segregated into 
a gas layer on the top and a water layer at the bottom and that gas occupied a 
much smaller fraction of the vertical cross section than water. As a result, 
only the thin top layer was miscibly flooded by gas, whereas the bottom por- 
tion was waterflooded. This phenomenon was also observed by other inves- 
tigators-4°6>”, 

Stalkup°® hypothesized that injection of solvent with low viscosity 
into a vertical cross section that contained oil more dense than the solvent 
led to the formation of a gravity tongue if the flow velocity was low enough 
so that vertical transport by gravity segregation dominated the viscous forces 
that induce instability. As the flow velocity increased, gravity forces played a 
diminishing role, and eventually viscous forces dominated. 


Van der Poel®’, Miller°*, and Warner°® studied miscible displacements 
in horizontal reservoirs. They all observed that when oil was displaced from 
a horizontal reservoir by a solvent of lower density, the latter tended to over- 
ride the former in the shape of a tongue owing to gravity segregation. Fayers 
and Newley”’ noted that the tendency for solvent to rise became more pro- 
nounced when the flow rate was reduced. Miller°* observed that the gravity 
override tongue pushed the reservoir oil down into the center of the reser- 
voir which was then slowly displaced immiscibly by the water injected along 
with the gas. Based on this observation, he concluded that, in field opera- 
tions, when the economic limit is reached, a substantial amount of the oil 
displaced by the miscible fluid early in the reservoir life is left unswept at 
abandonment. 


According to Miller°4, there are two important mechanisms causing 
the formation of an overriding gravity tongue. They are the density differ- 
ence between the two miscible fluids, which causes the injected lighter fluid 
to override the in-place fluid, and the vertical counterflow segregation result- 
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ing from a density difference between two immiscible fluids which results in 
the migration of the lighter hydrocarbon fluid to the top of the reservoir. 


Miller°4 also investigated the effect of reservoir thickness on gravity- 
controlled displacements. He concluded that if a gravity finger was formed by 
density difference between two miscible fluids at the interface when a lighter 
fluid was injected, an increase in reservoir thickness would lead to more 
overriding because there was a greater gravity potential gradient at the top of 
the reservoir than at the bottom of the reservoir. However, if the gravity fin- 
ger was caused by counterflow segregation, then a decrease in reservoir thick- 
ness would lead to poorer performance since the fluids would have a shorter 
distance to travel to segregate. He also suggested that capillary forces could 
prevent segregation, or at least retard the rate of segregation, in very thin 
sands. 


Spivak®? investigated gravity segregation in the two-phase displace- 
ment process, using a three-dimensional, three-phase and incompressible 
simulator, and found that gravity segregation increased with increasing per- 
meability (either horizontal or vertical), density difference and viscosity ratio, 
but decreased with increasing rate and level of viscosity for a fixed viscosity 
ratio. In contrast to Miller’s conclusion’, Spivak found that gravity segrega- 
tion in a 6 m thick reservoir could be as severe as in a 60 m thick reservoir. 
Based upon this finding, he concluded that formation thickness alone did not 
determine whether gravity segregation would be a problem. 


Araktingi and Orr ©! used a simulator to study gravity segregation in 
miscible displacements in a vertical cross-section. Their study confirmed that, 
in homogeneous porous media, better displacement performance could be ob- 
tained at high viscous-to-gravity force ratios for any mobility ratio. As a re- 
sult, high recoveries were reported at increasing viscous-to-gravity force ra- 
tios. 


62 conducted six horizontal 


Thomas, Bergins, Monger and Bassiouni 
cyclic carbon dioxide corefloods to investigate the influence of gravity segrega- 
tion on oil recovery. They concluded that any effect of gravity segregation on 
the process was caused mainly during the huff (injection) stage and that grav- 


ity override benefited process performance. A few vertical cyclic carbon diox- 
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ide corefloods were carried out also. In these experiments, carbon dioxide was 
injected at the bottom. Gravity segregation was found to promote oil contact 
by facilitating the deeper penetration of carbon dioxide into the core resulting 
in improved oil recovery. 


Using a compositional simulator to investigate the effect of solvent 
composition on displacement performance in two-dimensional (x-z) flow 
where both viscous and gravitational forces cause transverse transport or 
cross flow of fluids, Pande®* concluded that the magnitude of the crossflow is 
affected by the oil-solvent density ratio and that the driving force for gravity 
crossflow is provided by the difference in fluid densities. 


Gravity segregation has not only been studied in horizontal flood ex- 
periments, but also in vertical flood experiments. Hill®*, by conducting exper- 
iments on the displacement of sugar liquor by water from columns of granu- 
lar bone of charcoal (sweetening off), came up with a critical velocity expres- 
sion which predicts the rate above which instabilities can occur and viscous 
forces can dominate for vertically downward displacements. At rates less 
than this critical rate, gravity forces dominate and instabilities will not occur. 
The equation, defined in terms of the viscosities and densities of the two flu- 
ids, also accounts for the channelling which sometimes occurs when one 
fluid follows another along a uniformly packed column. 


Using the same idea as did Hill64, Dumore’® developed an analytical 
model for predicting a gravity stable rate for a miscible displacement which al- 
lows mixing between the solvent and oil. 


Slobod and Howlett®® used the critical velocity concept to determine 
the critical flow velocity for their vertically downward miscible displacement 
experiments. Their experiments were performed using fluids of various den- 
sities and viscosities. The results indicated that gravity segregation could act 
to shorten the mixing zone developed between the displaced and displacing 
phases when the displacing material was the less dense phase, and lengthen 
the mixing zone when the displacing phase was the more dense phase. 


Experimental studies performed by Tiffin and Kremesec®” showed that 
oil recovery at breakthrough was nearly doubled when carbon dioxide was in- 
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jected at the top rather than at the bottom in vertical miscible displacement 
experiments. A comparison of vertical and horizontal displacement experi- 
ments was also made in their studies. Downward gravity-assisted displace- 
ments, even at velocities much greater than the critical rate, were shown to be 
more efficient than horizontal carbon dioxide displacements at comparable 
rates. Fong, Tang, Emanuel, and Sabat®® recently conducted the same experi- 
mental studies as did Tiffin and Kremesec. Their results indicated that, at the 
same oil recovery, the producing gas-oil ratio of a vertical core flood experi- 
ment where gas was injected at the bottom, was twice that of a horizontal 
displacement. 


3.3.1 - Effect of Gravity Segregation on Saturation Distribution 


In 1951, Martin®? performed a theoretical study on gravity segregation 
problems during two-phase flow in porous media. He concluded that only a 
short time would be required for an initially uniformly distributed gas 
saturation to reach a distribution in which most of the gas, which is free to 
flow, is segregated at the top of the reservoir. 


During the course of their experiments, Richardson Perkin®? observed 
that the volume of gas concentrated at the top of the sand due to gravity seg- 
regation was displaced by oil which in turn was displaced by water. 


Walsh and Moon” theoretically analyzed the immiscible displacement 
of oil by water injection in a down dip reservoir. Their results showed that 
the final oil saturation increased as the gravity number increased. In their 
(Po—Pw )kgsina 

oT 
Pw are, respectively, the densities of oil and water; and k, g, Uo and ur are, re- 


study, the water-oil gravity number was defined as , where Pp, and 


spectively, the absolute permeability, gravity, oil viscosity, and total velocity. 


Hovanessian and Fayers’! studied gravity effects in linear waterfloods 
and concluded that the gravitational forces had a pronounced effect on the 
saturation profiles. 


Recently, Pande® and Belgrave and Win” used a three-phase, three- 
dimensional black oil simulator to study gas saturation distribution in homo- 
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geneous reservoirs at various viscous to gravity force ratios. They showed 
that, on moving upward from the bottom to the top of the reservoir, the gas 
saturation increased and decreased on moving outward to the production 
wells. 


Cook’, analyzing gravity segregation performance during natural de- 
pletion, concluded that gas liberated from solution can segregate to the top of 
sand developments in such a manner that a zone of high gas saturation 
forms along sand tops and is underlain by a zone of high oil saturation where 
flowing gas-oil ratios are maintained very low. 


3.3.2 - Effect of Segregation on Sweep Efficiency 


It is well established°?,©9.74-75 that gravity tends to decrease vertical 
sweep efficiency in horizontal homogeneous reservoirs because of gravity 
tonguing. Belgrave and Win” observed that, for homogeneous reservoirs, 
vertical sweep efficiency increased with increasing gas injection rates. Re- 
cently, Ekrann’® proposed an analytical model to predict the oil recovery at 
breakthrough with gravity effects for a waterflooding process, using vertical 
equilibrium theory. Using his model to generate results, he concluded that, 
in two-layer horizontal reservoirs, gravity may have a profound effect on the 
vertical sweep efficiency of a waterflooding process. 


Stone’” hypothesized that it takes some time for gravity segregation to 
occur; thus, there is a region near the injection well where vertical confor- 
mance is good, and the size of this region is determined principally by the in- 
jection rate, the vertical permeability, and the density difference between wa- 
ter and gas. Using this hypothesis, he proposed an analytical model to predict 
the size of this zone and, hence, the vertical sweep efficiency of a WAG flood. 
In deriving this analytical model, Stone assumed that a mobile gas zone was 
formed at the top of the reservoir and a mobile water zone at the bottom. 


Craig et al.°* measured vertical sweep at breakthrough in homoge- 
neous and isotropic cross-sectional models. Their experiments showed a re- 
duction in vertical sweep for increasing values of mobility ratio and decreas- 
ing viscous-to-gravity forces ratio. 
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3.3.3 - Effect of Reservoir Heterogeneities 


Studies have shown that gravity segregation of the fluids is sensitive to 
reservoir heterogeneities®?°°57,”2,_ A simulation study on gravity segregation 
on a two-layered reservoir conducted by Spivak®? showed that gravity effects 
were diminished when the low-permeability layer was on the bottom and in- 
creased when the high-permeability layer was on the bottom. Warner?® ex- 
plained that if the high-permeability layer is located on the top of the reser- 
voir, it provides a high-conductivity path to transport the carbon dioxide 
rapidly to the producer, once carbon dioxide has segregated into this layer . 


Spivak®® and Warner°® also investigated gravity segregation in ran- 
domly heterogeneous reservoirs. Warner®6 reported a nearly doubled oil re- 
covery for the carbon dioxide slug process when random heterogeneity was 
used in place of homogeneity in his simulation study. They both agreed that 
the random heterogeneity tended to reduce gravity segregation. 


Araktingi and Orr®! also showed that for injection of a light solvent 
into a layered reservoir with high permeability low in the reservoir, better 
displacement efficiency was observed at intermediate values of viscous-to- 
gravity forces ratio than at high or low values. If the high permeability was at 
the top, a high value of viscous-to-gravity force ratios was preferred. 


3.3.4 - Field Observation of Gravity Segregation 


Wilcox, Polzin, Kuo, Saidikowski, and Humphrey’® noted the overrid- 
ing of the injected gas in the Prudhoe Bay Miscible Gas Project in Alaska due 
to high vertical permeability and the large difference in density between the 
gas and reservoir fluids. The injected gas rose to the top of the reservoir or 
underneath shales, forming cone-shaped swept intervals around WAG Injec- 
tors. 


Pritchard and Nieman’? reported segregation of the injected gas after 
examining the dynamics of a hydrocarbon miscible flood in the Judy Creek 
field through the inspection of saturation cross-section between the injector- 
producer pair. Early in the solvent injection life, the solvent displaced an oil 
bank from the swept zone out into the reservoir. With further miscible injec- 
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tion, the solvent tended to override the water flowing interval and traveled 
along the top of the reservoir to the producers. This solvent then displaced 
oil that resided in the upper portion of the reservoir. As long as injection of 
the miscible fluid was continued, oil continued to be produced from either 
within or below the gas flowing zone. When water was injected, the oil bank 
also migrated upward in the reservoir due primarily to density differences be- 
tween oil and water. The oil would then move into the previously oil swept 
gas flowing zone and re-saturate this region back to the residual oil saturation 
to waterflooding. This oil would become immobile and could not be recov- 
ered: 


3.4 - Background to Non-Equilibrium Processes 


Non-equilibrium thermodynamics, also known as irreversible thermo- 
dynamics, considers systems which are not in thermodynamic equilibrium, 
or which are approaching equilibrium. Historically, its roots are found in the 
phenomenological laws of viscous flow (Newton), heat conduction (Fourier), 
diffusion (Fick), and electrical conduction (Ohm). Its theory restricts itself to 
large systems that can be treated as continuous media and can be assumed to 
be in local equilibrium. That is to say, a non-equilibrium system in local equi- 
librium can be divided into cells small enough so that any thermodynamic 
properties — such as mass, density, pressure, and temperature — which in 
non-equilibrium situations can be functions of space and time, vary slightly 
over one cell. On the other hand, these cells must be large on a microscopic 
scale — large enough so that they can still be treated as thermodynamic sub- 
systems in contact with their surroundings. It should as well be kept in mind 
that local equilibrium cells must be open for energy and/or mass transport in 
order to account for the overall macroscopic time evolution of the system®. 


Non-equilibrium thermodynamics has many applications in petrole- 
um engineering. The idea of petroleum reservoir engineering simulation is 
analogous to the theory of non-equilibrium thermodynamics. The reservoir 
under study is divided into many blocks. In each block, the fluid properties 
and saturations, pressure, and temperature are everywhere uniform, but vary 
from block to block. In the immiscible carbon dioxide injection process which 
is the subject of this research, when carbon dioxide is injected into a heavy oil 
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reservoir, initially, because the concentrations of carbon dioxide in the in- 
jected gas stream and those in the reservoir fluids are not equal to the ther- 
modynamic equilibrium values, mass transfer occurs, tending to even out the 
concentrations. Furthermore, as was already mentioned in Chapter 1, the 
concentration of carbon dioxide near or around the injection well is higher 
than that away from or near the production well. Thus, mass transfer has to 
take place to even out the carbon dioxide concentration. Consequently, it is 
important to know the conditions and the time required for thermodynamic 
equilibrium to be reached. The time at which equilibrium is reached is 
known as the equilibrium time. Goss and Exall®! performed an experimental 
study on the equilibrium time of a carbon dioxide-bitumen mixture. An 
18,000 mPa.s sample of bitumen was left exposed to carbon dioxide at 6.8 MPa 
and 50°C. After 12.5 days, they observed that concentrations of carbon dioxide 
in the gas phase equaled those in the bitumen phase. 


Martin, Combarnous, and Charpentier®* conducted an experimental 
study on mass transfer and phase distribution in two-phase flow through 
porous media under conditions similar to hydrodynamic reservoir condi- 
tions. The gas phase was carbon dioxide and the liquid phase iso-octane. 
Based upon the experimental results and observations, they concluded that 
the mass transfer process was not instantaneous while the phases are flowing 
co-currently inside a porous medium, that the equilibrium lag varied as a 
function of the ratio of gas to liquid velocity, and that mass transfer variations 
were mainly due to specific interfacial area variations. 


3.4.1 - Background to Mathematical Modelling of Non-equilibrium 


In order to model the mass transport phenomena under non-equilib- 
rium conditions, a mass balance for each phase has to be written. A set of 
macroscopic mass balances with interfacial mass transfer phenomena were 
presented by Bird, Stewart, and Lightfoot®’. By considering each phase sepa- 
rately in a two-phase flow system, interfacial balance equations, and interfa- 
cial boundary, and using a time averaging procedure, Ishii®+ developed a 
mathematical model which consists of two sets of equations which govern 
the mass, momentum, and energy balances for each phase. This so-called 
two-fluid model was later shown by Kataoka and Ishii® to predict accurately 
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mechanical and thermal non-equilibrium between phases, if interfacial trans- 
fer terms are modelled accurately. 


Farouq Ali, Chakma and Islam®® proposed a non-equilibrium mathe- 
matical model for the simulation of alkaline/polymer injections. The model 
consists of mass balances for water, oil, gas, and chemical transport in the 
porous medium. A numerical simulation study was performed to investi- 
gate the time effects of absorption, dispersion, and interfacial tension. Bree®” 
developed a non-equilibrium thermodynamic theory of continuous media 
for simple mixtures. Global and local entropy inequalities were derived for 
the mixture and for each of its constituents. He found that the set of entropy 
inequalities for the constituents was not equivalent to the single entropy in- 
equality for the mixture. Mathieu and Lebon®® presented a non-equilibrium 
thermodynamic model for describing the flow of a two-phase fluid with 
phase transition. A non-zero slip velocity, as well as a temperature jump be- 
tween the two phases, was taken into account. They also devised a method 
for calculating the velocity of condensation which is based on the Knudsen 
number, which is the ratio of the molecular mean free path to the mean 
diameter of the particle. Lozada and Farouq Ali®? proposed a set of mass 
balance equations which governed non-equilibrium transport phenomena in 
an immiscible carbon dioxide flood. Interfacial mass transfer effects were 
taken into account. The two-film theory was employed to model the 
interfacial mass transfer effect. They obtained six sets of scaling criteria, which 
were used to design the experiment for a non-equilibrium immiscible carbon 
dioxide drive”?. Their equations were not complete because momentum, 
total energy , and entropy balances were not included. 
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4 - DEVELOPMENT OF THE MATHEMATICAL MODEL FOR THE 
NON-ISOTHERMAL AND NON-EQUILIBRIUM IMMISCIBLE 
CARBON DIOXIDE WAG PROCESS 


To-date there has been no mathematical model for an immiscible car- 
bon dioxide flood which includes non-equilibrium and phase change phe- 
nomena for non-isothermal cases, where the reservoir temperature is differ- 
ent from the surface temperature. The mathematical model previously pre- 
sented by Lozada and Faroug Ali®’, which consists of mass balances for each 
component i in each phase j, includes the non-equilibrium phenomena 
amongst the oil, gas, and water phases only. Their model is valid for isother- 
mal cases, where the reservoir temperature is the same as the surface temper- 
ature. A phase change was not considered in their model either, for the gas 
phase was assumed to consist of pure carbon dioxide only. Consequently, 
there is a need to develop a non-isothermal and non-equilibrium mathemat- 
ical model including phase changes and interfacial mass transfer, which was 
one of the objectives of this research. 


4.1 - Derivation of the Mathematical Model 


The following assumptions were made when developing the model. 
The process was treated as three-phase flow (oleic, aqueous, and vapour). The 
gas phase was assumed to consist of carbon dioxide and water vapour, while 
the oil phase consisted of carbon dioxide and heavy oil, and the water phase 
consisted of carbon dioxide and water. Mass transfer was assumed to occur 
from the gas phase to the oil phase and from the oil phase to the water phase. 
Mass transfer was assumed to result from diffusion, convection, and disper- 
sion. Interfacial mass transfer between the oleic phase and the vapour phase 
and/or the water phase for any component i was modelled using an interfa- 
cial film model. Darcy’s, Fourier’s, and Fick’s laws were assumed to be valid. 
The system was assumed to be in local thermodynamic equilibrium. Chemi- 
cal reactions were neglected. Negligible potential energy was assumed. Gas 
bubbles were assumed to be of a single and constant size. The velocity of the i 
component in the j phase was assumed to be that of the j phase. In- 
stantaneous phase equilibrium was assumed to exist at the interfaces between 
the oil and gas phases and the oil and water phases. No accumulation at the 
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interfaces was involved. Oil, gas, and water were considered to be Newtonian 
fluids. The viscosity of the gas phase was that of carbon dioxide. The above 
assumptions lead to the following mathematical model for a non-isothermal 
and non-equilibrium immiscible carbon dioxide flood. The diagram below 
presents a schematic of the system under consideration. 


water [---- 


In the mass balance for each component, the following terms are con- 
sidered: convective mass flux, diffusive mass flux, and interfacial mass trans- 
fer. 


Mass Balance for Carbon Dioxide in the Oil Phase®? 
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Mass Balance for Carbon Dioxide in the Gas Phase®? 
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Mass Balance for Carbon Dioxide in the Water Phase®? 
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(4.3) 
Mass Balance for Oil in the Oil Phase®? 
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| for Water in the Water Ph 


Pa) = =* O 
S| CuPous + Sy (D. +D, (Cue ) “ 
rf) —— =* O 
P| Coa ny ae Sy (5. = BE w } (Coma ] a 
2 ) 
2 | CaP ar +S, (D. +t: D, whe SE, (es oe | + Nee NEN nays 


+N gy ole = ($8 uP wCw.w) (4.5) 
Mass Balance for Water in the Gas Phase 
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where, 
N 7 = tate of mass flux of water from vapour phase to liquid phase, 
wWw,et- 


N /-s¢ 7 Tate of mass flux of water from liquid phase to vapour phase. 


It is assumed that when water vapour condenses it has to give up a 
certain amount of latent heat which is enough to vaporize the same amount 
of water to the vapour phase. Therefore, NN Ses z can be equated to NY Tape 
Adding these two mass balances for water will eliminate these two terms and 
result in one single mass balance for water, which is as follows: 
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Note that in the case where the condensing gas phase is carbon dioxide 
and the vaporizing liquid phase water, the assumption as stated above is not 
valid because the latent heats of carbon dioxide and water are not the same. 
In such a case, in order to determine the amount of water that vaporizes for 
each mole or mass of carbon dioxide that condenses and/or diffuses, one 
needs to perform a heat balance at the interface using the latent heat as shown 
below. 
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Knowing the latent heats, the relationship between the two fluxes can 
be determined from the above heat balance. 


In the momentum balance for each component, the following terms 
are considered: momentum gain by convection, diffusive momentum flux, 
gravitational force, pressure force, drag force, buoyancy force, momentum 
gain by viscous transfer, and momentum gain by interfacial mass flux. 
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Momentum Balance for Carbon Dioxide in the Gas Phase 
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In the total energy balance for each phase, the following terms are con- 
sidered: convective transport of kinetic energy along with the fluid across the 
system boundary, kinetic energy of diffusion, work done by pressure forces, 
work done by viscous forces, work done by mass forces, reversible conversion 
of kinetic energy to internal energy due to viscosity, irreversible conversion 
of kinetic energy to internal energy due to viscosity (i.e., the energy that is dis- 
sipated), work done by drag forces, work done by buoyancy forces, convective 
heat and kinetic energy transport due to interfacial mass transfer, and convec- 
tive/conductive heat flux. 
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The following balance equations express the changes of the entropy of 


each component with time, which are due to the flows of entropy into the 
volume element and the presence of an entropy source due to irreversible 
phenomena inside the volume element. The following terms are considered 
in the entropy balance for each component: convective entropy flux, diffusive 
entropy flux, entropy created due to convective heat flux, entropy transport 
due to temperature gradient, entropy transport due to chemical potential 
gradient, entropy created by work done by buoyancy forces, entropy created by 
work done by drag forces, entropy created by work done by mass forces, 
entropy flux due to velocity gradient, and entropy flux due to interfacial mass 
transfer. 


Entropy Balance for Carbon Dioxide in th il Ph 


3 | So (DB. + Deo,.} 


2 (Coo, oPoPox8co, } Moe Te ii (Coos 9Pe) 


A a 08,(D. + Deo, a 
+57 (Cco.PoPey8co. ae amen era 2107.0 5y (Coon. o) 


a 0s.(B. +Doo,.0) 


> = Nco,,0 “ (Cco,,0Po} 


e) 
a ay (Cco,,oPo%0z8co, + 


0 hco Pa) hco ah hco, 
ax P (Coo, "2.94, + 2 (Ceo, aeeraaeacy a Cco,, oPo 44 ~ Boz 


1 0 (1 Ona. 
}+Coo,, oPolico, Voy =— ay (= + Coo, oPohco, Voz =— os (= 


) 
+Cco, Poco, Vox =— ax (= 


ee Bi f : 


ci “ae > a ¥) > 6 i 
: ae pete: gp! Wi agai Ae a sali pr * ruee Ms 
woes 
rpg) @ (7) i (citi fay 4 snes ~ i” Ns ait ib ava Asin a ean Vi R F In, io 5 fo . 
; . —T vi 
whined’ é ah oan Qe ny re ie to poem Pr) San kivectnty: niguo 
. aisha ‘ - rsftvé de : : 
by get inae AM PAL: 4 bait oni, 6 $14 9 0ef4) = Auiteniae: cat oti ws ig 
- 
ithe) rar ad ky whut 1? wee: fie AON priya. rie 13 ‘nolo yore a. 


~9t' au YoOWh mss. | work AVI Vrs IS SEes fess oe i ; 
i 17 FY : f wit ah ji Wendie V ; Me) Tete ) site th q . oie i aos 


° 
' 


oe Ao we) ot rey eee Water ety OR aad. ny vil GA 
Gon ate ee Onbb Mew vd We ay WA ® ous on yay 
pre! i peut Fy i; + AL Ss Wy Puar’y Ti 7 Polis: ; Th Sante " ous nh 


‘= 


ae | = 
Ms y & yo oa 
: , 
l 
tof J ar ’ anit saat 
{ 7 
7 LP 
, » 7 yy 
' g°.ot f 
p 
oA Ss ty i 7 a be / he 
ir 2 ‘ ’ Fie, cr iz N 7 ha 7 


ize er) 


r 


Ae, A er ce eens 


re?) 


a cr 3 

+98.(D. ots Deo,.0 }5-(Cc0,.oPo)- ra ae 
NCO,,0 
“| 


Nco,,0 3 Ceo, Po 
sk Jk 


a 
H}O 
No 
° 
Santi 


+S, (B. a cece \2 (Coo, ,oPo 


+05.(D, ti De0,.0 }=(Coo,.0Po) = (80.8% a DoyBy + Wed 


1 D D 1 ey 2 
a No, ogx®CO, 4 Nop, ogx —(95, = Dee) F Nco,,ogy8CO, $e N co, ogy 5a (Oxy aA dy) 


a 
ee? Dire? 2 
me Noo,, ogz®CO, + Noo,, ogz a =e -82,) + Nco,,owx®C0, + No, ,owx Sq (Own a ox) 


2 2 
,OWy »,OWY wy oy ,OWZ 9) 
SF Neco, Sco, ta Noo, x oe -% ar Noo, SCO 


1 ) 
"s Nco,,owz sq (Owe - 02, = 5 (®SoPsCc0,,08C0, ). (4.19) 
ntr Balance for Carbon Dioxide in th Ph 


) 08,(D, + Deo.) A) 
= (Cco,.2Pe%ex8c0, 53 | oemaameemmee o.'CO2 = (Coo, .2P¢) 


) 9 | Ss (5, ms Deo.) 0 
+ 5y (Ccon#P 2x80, sti oy T NcOp.g 5y (Cco.eP a) 


A) a) 08,(D, + Deore) 0 
"7 Ao (Cco,,2P 20 sz8co, u aS ree ee: = (Cco,.ePs) 


p) hco,,g F) hco,,g r) heo,,g 
+3 (Coosa Oat. 04, + 2{ Coosa = Bey asa Cco,,sPe 7 Ber 
ofl Onl ehi(al 
mT Cco,,gPghco, Bex = ap (=) 1 Cco,, sPehico, Ugy =~ dy (= }+Cco, gPelico, Usz =— oe (=) 
= =* O 0 {Nco 4 
+68,(D, + Deo,2) 2 (Corse) =( = ) 
= =* 0 d ( Nco 2 
+48,(D, + Deo, }& (Coors) (2 


e) Oeics, Cco,,gP 
z oS .(D: ar Whey, a (Coo, Pg) 2| arrae = }: EE (OgBx + OSes a 9 8x) 


Foo,, gBxU ex , Foor. sByU gy i Foo,, eBzU gz WecO7, gDx Vex ASisehy Dy“ gy 


Ee 
SS sts ht 


o/s 


_ Feo, gp2¥gz 1, OS ox are Oey ee OD gz 
ip dp RAL rah Tear oy, Tiss ae GZ 


1 ae Ee poll OURO) ot \Fol OU Cire 
+Etpo| mn ae) Tk zea Cova TS P| cox) nde 


2 2 
+Nco,,g0xco, +Nco,,goyco, th Noo,, goz©CO “ae ap (Oe ba 02x) 


+Nco, ,goy == (8% - 92) +Nco,,¢02 == (82 - 92, )= © (4p484Cc0s.45C0, (4.20) 


Entropy Balance for Carbon Dioxide in the Water Phase 


S 2) LGBy. 
A) A) ) A w 60,,0 | a) 
5 (Cco,,wPwwx8c0, } ea = faith ements 1CO2.W 5, (Cco,,wPw) 


3 | S+(Dy + Deo.» 


e) 
ce z Nco,,w 5y (Coon) 


+5 (Cco,.nPwns8c0, os 


a) P) S.(Dw “Ue Deo,,. 0 
i: 5, (Cco2.wPwPwz8c0, ) eal Dr a = (Cco,.wPw)} 


?) hco, 0 hco () heo 
NT Gs See) -——hG a pre ed) ( = we By 
2 CO,,wPw ie ae a CO,,wPw TH et axl C02: Pp T z 


a) ps onal aes 
+Cco,, wPwhco, Dux 5— ax (= }+Cco, wPwhico, Dwy =~ oy (= J Cco,, wPwhco, Owe =- a (=) 


+S, (B. +Deo,, wa g (Coo, wPw | g [Fees 


ox\ T 
+98u(Dy + Deo, “= (Coda?) => 2 (eos 
T 
Chey, Nco,,w 
~ Su( Dy + Deo,» JS (Coo, wPw} a 2 eda 
+ —oepinbe (Sux8x + OG, + Cee +Nco »,.wox®CO, = NCOswor = (Px a ox) 


1 (2 2 
+Nco,,woy®CO, +Nco,,woy ap (Oey im Onn) a Nco,,woz8co, 


1 C) 
© Nco,,woz 57 (Dive ro 02,) = 5 (PSwPwCco, W800, ). (4.21) 


dn 


ell 


t 
t 


j 7 7 ww (oo 
Le acess 5 eet, 


. 


acl 


r lance for Oilin the Oil Ph 


9 | So (D. . Ds) 3 


) 
ae (CR ee ) sd Tsar eee an 


A) a 8.(D.+ Da. | a) 
+5, (CocPoPey8o) + oy a tee Cott, 


A) OS, (D. a Da.) 0 


) 
3F a LG psiase) - ay Tinea a ie ae [expe 


0 h 0 h 0 | h 
a aaa I ay as FT GY Oy, sant Ge 2.3 


1 da al ) 
=a =] Co Poh Oy dy & af Co, oPo h,Do2—— ae (= ) 

Bye Ghee 2 (Ty 
ox (@roa) ox (3 :) +08(D, ats D, ole (C oO eon ae ) 


=— = C 
iE 8,( De rf Doo) (en Po) = & (Ze + =22F? (85,8; + Ooky ar 90282) 


+—T (Bas) ae Oey fle Core auc Re 
OXX ox T oyy dy OZZ Oz Oxy Oy ox 


i 
ak T 
1 GUS OU nal O9ox , Iox 
~ ae [ + 22 + Fan = = ax) = 8 5, (PS oPoCo, aS. ; (4.22) 


Entropy Balance for Water 


Su(Du + Dew 
ct Geo oaeay) ah g 


Ox Ox 4t ox 


) 
+5, (Cu.wPwwy8w) oF eae 


+ al, (Co Aa es ea) a, 


0 h 0 h 0 [ he ) 
, peel —!C —~ 38 ee 4 NC ath eos, 
ar a Cue T Pw a ay ( w,wPw T wy aS ,wP T 


*y 


i: > on 
\ i ' ,% " a 7 
a deal ait ie} ae : : 


ye) 


pi i 
sf a 5 a "4 yy) and 


in 


+Cy, wPwh yD (=) ane Py walw wy & (= ) ne Cy, wPwhw Ow <(=) 


pa mnent De) < (7 ave (D.. +B... )2 aye tS 


=—* JE 2 Ww, ae 
; )2(2#)+5 [ome + OwBy 2 bees e| 


(28s0).1 1 aes po fl ) 1 OO wy 
cies C waxes | tas Uwyy ok += Tw77z + Twx 
Ox oy if Oz u i Ox 


L 

di 

1 Wy IWy,) 1 Wy, Wy). 9 

tpt ett Sa) ptm G+ Be) +5 (Cove ete) 


2s 
<= 
N 
= 
TS 
Ol 
= 
ahs 
S) 
sits 
ae 
< 
Be) 
ee 
vie 


) 
ae a) +55 (CusPePnv8we) 


7 | he Oual 
poe bE ae) + (CusPs ial v, 2 | +CvgPaliwed yx =(=) 

C h : ie 
an w,gPg oes T + Cy, gPehwedez Ze T 

Q : shits 

+S (Ds; +Dvs)—(C wa?) (= i) Ds (B, +B...) (Cc waPa)5y(E] 
(+S (D: +Das le (c w,gPg 2 
+ (Sy —8 Sy )tit + in—{( 95 = 8x)+(8 


a £ ($8 uP Cw, n8 + GS, prCuyena} Soe 


p 
y Cy ae E( [See + DoBy a 9 9282) 


ntr lance forth k Matrix 

Oskreo lied (kp OD 2 (Au ( T(z) oT =(3) 

oe Aad capac [a eal sel oa, | ee k om 

a (ae ot). 2 (Hs at) 4.3 ib Gy uF Ar Ox JOx\T i br Oy dy \T 
( 


Cian al ds 
alae) ete) Ps eG r. 
+(ns Saleen) 


a 


These partial differential equations represent the conservation of mass, 
momentum, total energy, and entropy for the non-isothermal and non-equi- 
librium immiscible displacement of oil by carbon dioxide and brine. They re- 
quire additional constitutive relationships, constraint equations, and initial 
and boundary conditions. Table 4.1 (on the next page) lists the constitutive re- 
lationships and constraints. Appendix A contains a possible set of boundary 
conditions. 


33 


ane 5— | 
<tpe-tioe re. = 


ws ie aie 


sien panty ache Gem ted 
niet Swit ee 


| oa) 
ah) WU es 


Table 4.1 - Constitutive Relationships and Constraints 
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59. Z=Z(x,y,z) fixed coordinates 


According to Perkins and Johnston*®, the dispersion tensor consists of 
molecular diffusion and convective terms. Bear and Buchlin?! conducted a 
thorough study on transport processes of materials in porous media and 
showed that for porous media that are isotropic with respect to dispersion and 
where an orthonormal coordinate system parallel to the directional flow is 


chosen, the dispersion tensor becomes 


QD, oO 0 
Q=| 0 Dr 0 |. 
® WG 


(4.25) 


Also, Newtonian fluids were assumed. Therefore, the following rela- 


tionships for the shear-stress components are true. 


TT, =—2 00x 2 09, , Wy 00, 
a ee ert Oyen ozs)? 
00 00 
t = — Jy 2. en OD, Oy Od, : 
ie OVE eT Oxeeeoy = 07Z 

oo, 2 {d6, od, ad 
a) y Pasta bl dite eae, Med ea A , 
ae pe 2 y( er 4 Be 
AY map rane | Ware fag y” 
Pe Ey al ae oy } 
Ty = Txz = (as) 
VA OE VA) OF LL ax dz ; 


For mass transfer across the inteface, the two film theory was assumed 
to be valid. The theory assumes that the phases are in equilibrium at the 
actual points of contact at the interface, and the bulk of each phase being well 


mixed; thus, the flux equations can be written as®? 


(4.26) 


(4.27) 


(4.28) 


(4.29) 


(4.30) 


(4.31) 
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Nig =k]a pS (PiCi =|) ice) (4.32) 
where, 
Nj) = mass transfer rate of component i in or out the j phase through 
the jl interface, 
Kk = local mass transfer coefficient of component i in j phase related to 
the jl interface, 
a,, = interfacial area between phase j and 1, 
SF 
cH = mass fraction of component i in j phase at the jl interface. 


= average mass fraction of component i in the bulk of j phase, 


Furthermore, Darcy’s law for flow in porous media was assumed to be 
valid. Thus, Darcy’s equation, which can be expressed as follows, applies. 


k . 
Oj= aye j+P jgVZ). (4.33) 
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4,2 - Derivation of Scaling Groups 


This section deals with the derivation of similarity groups which are 
used to design experiments for a non-isothermal and non-equilibrium carbon 
dioxide WAG process. There are two approaches for deriving the similarity 
groups: Inspectional Analysis and Dimensional Analysis. Inspectional Anal- 
ysis requires the variables in a set of equations which fully describe the whole 
process. In almost all instances, the equations describing a flow process are 
partial differential equations, and the relevant boundary and initial condi- 
tions need to be specified for the process to be completely described. Dimen- 
sional Analysis requires knowledge of all of the relevant variables influenc- 
ing the process. This method is employed to derive the similarity groups 
when the partial differential equations describing the flow process of interest 
are not known. Geertsma, Croes, and Schwarz?? showed that Dimensional 
Analysis often yields a larger set of similarity groups than inspectional anal- 
ysis, but the physical meaning of the similarity groups themselves is more 
apparent from Inspectional Analysis. For this reason, Inspectional Analysis is 
often preferred. Dimensional Analysis is always used in conjunction with 
Inspectional Analysis to ensure that important groups are not omitted. Both 
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of these methods have been used extensively by many other researchers to de- 
sign their physical models*?*"?”. 


4.2.1 - Inspectional Analysis 


The procedure to derive the similarity groups by this method is to first 
express the governing partial differential equations, initial and boundary 
conditions, constitutive relationships, and constraints in terms of the dimen- 
sionless variables and their reference quantities. In other words, it is neces- 
sary to divide each variable and or property by some characteristic reference 
quantity. The following example can clearly illustrate this step. The general 
property M can be written as M = MpMr, where Mp is the dimensionless 
form of the property M, and Mr some constant characteristic reference quan- 
tity. Replacing M with MpMk in the governing partial differential equations, 
as well as in the constraints, constitutive relationships, initial and boundary 
conditions, results in an equation with the form as shown below for the mass 
balance of carbon dioxide in the oil phase. Appendix B lists the dimension- 
less form of the momentum, energy and entropy balances for carbon dioxide 
in the oil phase. 


Mass Balance for ©O> in the Oil Phase in Dimensionless Form 
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For other components, their dimensionless balances can be written in a 
similar manner. The reference quantities in the braces represent the coeffi- 
cient for that term in the equation. All coefficients have the same units. To 
obtain similarity groups, the entire equation is divided by any one of these co- 
efficients. This gives the dimensionless forms of the equations and their sim- 
ilarity groups. Similarly, the dimensionless forms of other governing equa- 
tions are obtained in the same manner. They are included in Appendix B. 
Table 4.2 lists the similarity groups and their physical meanings. 
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Table 4.2 - Similarity Groups From Inspectional Analysis. 


Now, let xXp=L, yp=W, and zp=H 
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Group 
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SoRHoRPRDLCO,,oR 


korPor 


Physical Meaning 


porosity 


geometric scaling factor 

geometric scaling factor 

time scaling factor 

ratio of gravitational to viscous forces 
ratio of gas to oil density 

ratio of water to oil density 

gas-oil mobility ratio 

water-oil mobility ratio 

ratio of gas to oil saturation 

ratio of water to oil saturation 

ratio of initial water to initial oil saturation 
ratio of production to oil pressure 


ratio of oil-water capillary to gas-oil capillary pressure 


ratio of water to oil pressure 


ratio of CO> partial pressure to water vapour partial 
pressure in the gas phase 


ratio of dispersive forces to viscous forces 
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Table 4.2 - (Continued) 
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Physical Meaning 


ratio of transverse dispersion to longitudinal dispersion 
of CO in the oil phase 


ratio of longitudinal dispersion to transverse dispersion 
of CO? in the gas phase 


ratio of longitudinal dispersion to transverse dispersion 
of CO) in the water phase 


ratio of longitudinal dispersion to transverse dispersion 
of oil in the oil phase 


ratio of longitudinal dispersion to transverse dispersion 
of water in the water phase 


ratio of longitudinal dispersion to transverse dispersion 
of water vapour in the gas phase 


ratio of longitudinal dispersion to transverse dispersion 
of gas in the gas phase 


ratio of mass transfer rate of CO in or out of the oil 
phase through the oil/gas interface to viscous forces 


ratio of mass transfer rate of CO in or out of the oil 
and gas phases 


ratio of mass transfer rate of CO? in or out of the water 
and oil phases 


ratio of CO, concentration in the water and oil phases 


ratio of CO concentration in the gas and oil phases 


ratio of water vapour concentration in the gas and water 
phases 


ratio of buoyancy force to mass force 
ratio of buoyancy force to dragforce 


ratio of the rate of water condenses or evaporates to 
viscous forces 
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Table 4.2 - (Continued) 
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Physical Meaning 


ratio of CO2 enthalpy to oil enthalpy 
ratio of water vapour enthalpy to water enthalpy 


ratio of CO? internal energy to oil internal energy 


ratio of water vapour internal energy to water internal 
energy 


ratio of the energy stored in the rock to that in the oil 


ratio of amount of latent heat released/absorbed when 
water condenses/evaporates to water enthalpy 


ratio of the rate of water condenses or evaporates to 
water injection rate 


ratio of convective to conductive heat transfer 


ratio of the longitudinal dispersive entropy flux of CO2 
in the gas phase to CO? entropy 


ratio of entropy gained or lost by CO? due to convec- 
tive heat transfer to CO2 entropy 


ratio of water entropy to oil entropy 
ratio of CO, entropy to water vapour entropy 


ratio of the entropy of the rock matrix to water entropy 


ratio of the conversion of gravitational forces to en- 
tropy to oil entropy 


ratio of CO> injection rate to viscous forces 


ratio of water injection rate to gas injection rate 


ratio of the momentum of the injected CO} to the longi- 
tudinal dispersive momentum of CO? in the oil phase at 
the injection well 
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Table 4.2 - (Continued) 


Group Physical Meaning 
Mco,R 
away rma ratio of CO2 momentum to water momentum 
wR 
53 £co rb ratio of the kinematic energy of CO to conduction heat 
’ kiR TRA injR transfer at the injection well 
54. _ Lwrb ratio of the kinematic energy of water to conduction 
Kirr TrAinjr heat transfer at the injection well 
55 (Swor — SwR )rhtp ratio of the entropy released/absorbed when water con- 
GMa oes 0n oan denses/evaporates to water entropy 
56 8co,R Weo,rb ratio of the CO? entropy to the entropy created due to 
i kip A ink conduction heat transfer at the injection well 
57. Ur ratio of entropy gained or lost by the rock due to con- 
TRS:R ductive heat transfer to the rock entropy 
58. ToiR ratio of initial oil and reservoir temperatures 
TR 


4.2.2 - Derivation of the Relaxed Scaling Groups 


The aim is to obtain sets of similarity groups which can be used with 
the scaled physical model, which was previously designed and built by Rojas!. 
The derivation of these groups is found in Appendix C. 


Approach No. 1 


In this approach, it is assumed that both the reservoir and the model 
have the same porous medium, fluids, pressure drops, temperature, and ge- 
ometric similarity is satisfied. The advantage of approach is that it leads to 
the satisfaction of the scaling requirements for viscosity, density, solubility, 
diffusivity, equilibrium constants, mass transfer rate, condensation and/or 
evaporation of water, and other properties which depend on pressure and 
temperature. Using the same porous medium (same porosity, same perme- 
ability, same grain size, and same wettability) allows the residual saturations 
to be the same in the model and in the prototype. Therefore, viscous forces, 
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capillary forces, and diffusive forces are properly scaled. This leads to the fol- 
lowing set of relaxed scaling groups. 


« 2 
H L HorSorDce0,,0R ORSoRMoRL? Neo,,ogRborL Neo,,owr Nco,,gor 
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In addition, it should be noted that all the dimensionless properties 
must be the same function of their dimensionless variables for the model and 
the prototype. For a model reduced in length by a scaling factor “a” and 
employing the same fluids as the prototype, Apmax, Poi, Pprod, T, k, 9, Soi, Swi, 
and F must be the same in the model and the prototype. The parameters H, 
W, Weo,, Ww, must be reduced by “a” and t by “a?”. 


Approach No. 2 


In this approach, the reservoir and the model are assumed to have dif- 
ferent fluids and pressure drops, the same porous medium, and geometric 
similarity. This approach is rather similar to Approach No. 1. The only dif- 
ferences are the model oil and the experimental pressure and temperature 
conditions. The model oil is selected so that its viscosity at room temperature 
equals that of the prototype oil at reservoir temperature. The experimental 
model pressure and temperature are different from the reservoir pressure 
and temperature. 
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In this approach, the determination of the experimental pressure at 
atmospheric temperature is crucial, because it is hypothesized that the solubil- 
ities of carbon dioxide in oil at two different pressure and temperature condi- 
tions can be the same. In other words, the experimental pressure at 21°C is 
determined such that the solubility of carbon dioxide in the model oil at this 
pressure is equal to that in the reservoir oil at the field pressure and tempera- 
ture. 


In addition, in order to balance viscous forces (viscous x/viscous y) and 
viscous divided by gravitational forces while maintaining geometric similari- 
ty, the pressure drop is relaxed. It is proposed to use the same porous 
medium in the model as in the field. In other words, the model will have 
the same permeability, porosity, and saturation distributions as the prototype. 
The main weakness of this approach is that all properties (except solubility) 
which depend on pressure and temperature, and composition are not scaled 
because of the different pressure and temperature in the model and the proto- 
type. The resulting groups are as follows: 


L H opSortorl LorSorDco,, oR PorSRH PgR Pwr HorKgr Uorkwr 
ak ie Ww’ koRPortR PorKorFR PoR Te Por “Herkor “LwrKor 


Neo, ogRMorL” Nco,,goR Nco,,owR De05,2R Deo,,wR Dwerer Coo, eR 
Coo,,orPoRKoRPoR Nco,,0gR Nco,,woR Dco,,or Door | Dw,wR Cco,,0R- 
Cco,,wR Cwe,eR pe Swr Soin PCO, Pwr PprodR Poin Foo, ,gBR 
Co oR ‘ CywR Sor SoR “Swik Pwe,gR Por PoR  Pwik “Coo, RP gRER- 
Foo,.sBR _ Nco,,ogrhco R-Fr Aconr hwer Uco,R Uwer PrUip 
Foo,,gDR Cco,RPgRERDCO,,oRSoR erable eon mena Csrtion: 

(Rwer —Hwr)titgte SprDeoneRMco,stHer hor _ gRH COR Swe SR 
CywRPwRowrORSwR  TREcO,RPeRKgrFR " TrSco,R 8CO,RTR SoR  SwR SoR 
(wer — Swe )tigte 
Cy wRPwREwRORSwR 


In addition to these requirements, the dimensionless properties as 
functions of dimensionless variables should be the same in the model and 


Ms q”’ 


the prototype. For a model reduced in length by “a”, and considering the 


same fluids, 6, Soi, Swi and K must be the same in the model and the proto- 
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type. The parameters W, H, and t must be reduced by “a”. The parameters 
Wco, and Ww must be reduced by “a2”. 


4.2.3 - Dimensional Analysis 


The dimensional analysis approach to derive the similarity groups for 
a process is based on the Buckingham x-Theorem, which uses the Principle of 
Similarity. The procedure is to first select the relevant variables for the pro- 
cess. The similarity groups can be determined using the Buckingham z-Theo- 
rem. The general rule for this approach is that if there are n separate vari- 
ables and m primary quantities, then the set will be complete when there are 
(n-m) dimensionless groups. More details on this approach can be found in 
Ref. 98. The similarity groups derived by this method are given below. 
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sional analysis approach because gAp,, was taken as a variable. 


5 - EXPERIMENTAL APPARATUS and PROCEDURE 


This chapter presents a description of the apparatus, materials, and pro- 
cedures used in the present research. The first part describes the procedure 
for packing and saturating the model prior to conducting an experiment, and 
the second part gives details of how the diffusivity and solubility of carbon 
dioxide in oil, in the presence of nitrogen, were measured. A discussion of 
the procedure for conducting an immiscible WAG experiment is also pro- 
vided. 


5.1 - Experimental Apparatus 


Figure 5.1 shows a schematic of the apparatus used for the displace- 
ment experiments. As shown, the apparatus used in this study consists of the 
physical model, fluids and porous medium, fluid injection and production 
systems, and the data acquisition system. 


5.1.1 - Physical Models 


Two models: linear and two-dimensional, were used in the present re- 
search. The linear model was partially scaled while the two-dimensional 
model was fully scaled to the Aberfeldy reservoir in Saskatchewan. The simi- 
larity groups which were used to design the model can be found in Ref. 1. 
The linear model was built to act as a screening model for the two-dimen- 
sional model. It was 415 mm in length and 98 mm in diameter. Chevron- 
type fittings were used to seal the ends of the pipe, as well as to form the injec- 
tion and production ports. Figure 5.2 presents a schematic of the linear 
model. 


In contrast to the linear model, the two-dimensional model was more 
complex. Cross sections of the two-dimensional model are shown in Figure 
5.3. Much effort was expended in designing and fabricating it!. A brief de- 
scription of this model is given below. 


- Rectangular shape: 45.7 cm x 45.7 cm x 2.2 cm. 
- Three reinforcing members. 
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Figure 5.2 - Schematic of the Linear Model. 
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Reinforcing Members 
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Reinforcing Members 


a) Horizontal Cross-Section 
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b) Vertical Cross-Section 


Figure 5.3 - Cross-Sections of the Two-Dimensional Model. 
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- Maximum internal pressure: 10.0 MPa at 51°C. 

- Maximum deflections of walls at 10.0 MPa: < 0.01 mm. 

- Weight of model: 1.0 tonne. 

- Number of wells: 9. 

- Possible patterns to simulate: 5-spot, 9-spot, line drive. 

- The model can be rotated for horizontal, inclined, or vertical floods. 


In this research, in addition to conducting experiments at iso-thermal 
conditions, it was intended to perform experiments at non-isothermal condi- 
tions. The field data that were used to design the scaled non-isothermal ex- 
periments are given in Table 5.1, on the next page. 


5.1.2 - Fluids and Porous Medium 
Oil 


The oils used in all experiments were from the Aberfeldy and Battrum 
South fields in Lloydminster, Saskatchewan and from a heavy oil field in 
Oklahoma. 


Carbon Dioxide Gas 


Commercial carbon dioxide gas with a purity of up to 99.9%, purchased 
from Medigas, was used in all runs. 


Brine 


Field-simulated brines were used in all runs. 


Porous Medium 


For most of the runs, Ottawa Silica Sand from Ottawa, Michigan, was 
used to represent the field porous medium since it has a grain size similar to 
that of the reservoir sands (70-140 mesh). 
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Table 5.1 - Field Information 


Field name 

Date discovered 

Producing depth 

Net sand thickness 
Original reservoir pressure 
Current reservoir pressure 
Production pressure 
Reservoir temperature 
Oil gravity 


Oil density (@ surface temperature) 


Formation volume factor 
Initial oil saturation 
Current oil saturation 


Oil viscosity at bottom hole temperature 


Porosity 
Permeability 
Geometry 
Well spacing 


‘Nn 
i) 


X 

November, 1987 
912m 

Salam 

8560 kPa 

4799 kPa 

2600 kPa 

3/G 

L7G 2 2A BI 

950 kg 

1.024 rm3/sm3 
60.7% 

44.0% 

160 mPa.s 
USVYs 

0.1 to 2.5 darcies 
inverted 9-spot 
161,875 m2 
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5.1.3 - Fluid Injection and Production System 
Gas Injection 


Carbon dioxide was injected using a Matheson gas metering system. 
This system controlled and measured the gas entering the model. The 
Matheson Dyna-Blender helped to control the flow rate of gas. A gas com- 
pressor was also used to maintain a constant gas injection pressure. A 
Matheson totalizer provided the cumulative volume of gas injected into the 
model. 


Oil Injection 


A positive displacement Milroyal pump was employed to inject heavy 
oil into the model. 


Brine Injection 


Brine was injected by a constant rate screw-type piston pump. The pump 
flow rate was controlled by varying the pump speed. 


Fluid Production 


The effluent was collected in a glass cylinder at atmospheric conditions 
(101.325 kPa and 23°C). Oil and water, because their densities were greater 
than that of gas, were collected at the bottom of the cylinder while gas dis- 
placed a volume of water in the upright glass burette equal to the total vol- 
ume of gas produced. Since oil and water mixed with each other at the time 
of collection, they had to be heat-separated to determine the produced vol- 
umes of each. 


5.1.4 - Data Acquisition System 


The production pressure was controlled by a back-pressure regulator 
which was connected to the production end of the model. Two Heise pres- 
sure gauges were used to measure the injection and production pressures. 


590) 


matVs SPR OaH ay, 
. on). isheay et cdiat 

ICD. wily AL) BRE OO feo 1 

A .crangstep well im eng antes 4 ih iach at 


ome 


; aoe “ou ek ee ~ 
cont Sepoupetie ey“ i 2th amma elt Bebe i a8 bees 
r : thn . ie sbi 
: ae Ru _ " 
v, a i; - 7 ee a ~~ 
- wy nowt : fe e 
i] a - : ar : eae 


Rau ° ; i 
@°° ' nl Ae Dern gs ay teyo:, 


' ' ett sri *) aia yweerse ofe.* Solana £ = wel nn ie 
’ | er cry ae tina aaa \ 


—s 


-4ar. ban sirens +n debs byoaanty Ue a) Bay wade 
hitita HO wesliiien eat aeul i Wew Ts i! ’ c 
aay actly plays: alt Va: Patch eet acts » bokoliei, 
at Si) 19 berg vids bai tigen, wif) at hig 
ete Mr eee) is wath trea 108 wt fh: 4 fim a é rhe 
Loy Geeta, nial dail nh mee: ret od ot Uh ap 


Wh 


— misty gn 
ae _ ae 


5.2 - Experimental Procedures 


In this research, as mentioned previously, two models: linear and two- 
dimensional, were used to conduct experiments. The terms ‘Linear’ and 
'Two-Dimensional' are used for convenience only; in fact, flow in any physi- 
cal model is three-dimensional. Much effort was made to minimize the effect 
of gravity, which acts in the third dimension. 


For both models, the experimental procedures used were identical, ex- 
cept that dry packing was used for the linear model and wet packing for the 
two-dimensional model. The procedures are as discussed below. 


5.2.1 - Packing 
Linear Model 


Dry packing was used for the linear model. The packing procedure is 
relatively simple. After the bottom Chevron end cap was installed on the 
production end of the model, the model was inverted so that the open 
(injection) end was up and so that it was perfectly vertical. A level gauge was 
used to check if it was in the vertical position. An air vibrator was then 
strapped on the side of the model. Next, Ottawa sand was slowly poured into 
the model while it was being vibrated. In this way, a tight sand pack was 
achieved. Afterwards, the model was left vibrating for eight-to-ten hours. 
After vibration, the top Chevron end cap was mounted, and a vacuum pump 
was connected to the top end to evacuate air from the model while it was 
again being vibrated for another eight-to-ten hours. A vacuum was drawn at 
the top to achieve the best vacuum possible. At this point, the model was 
ready for pore volume determination. 


Two-Dimensional Model 


While the linear model was dry-packed, the two-dimensional model 
was wet-packed for convenience. Similar to the linear model, the two-di- 
mensional model was first inverted so that the open cavity was facing up. 
Next, an aluminium extension was temporarily mounted on the top of the 
model, and distilled water was added to the model. The purpose of the exten- 
sion was to maintain a 10-cm head of water above the sand level. An air 
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vibrator was clamped on the top of the model and activated, and Ottawa sand 
was slowly poured in until the sand level was about 2 cm above the head of 
the model. The model was then vibrated for at least eight-to-ten hours. 
Afterwards, the 10-cm head of water, air vibrator, and aluminium extension 
were removed, and the top flange was put on and bolted. Finally, the model 
was pressure-tested at about 6.0 MPa or higher to check for leaks. If no leak 
was detected, the model was now ready for pore volume determination. 


5.2.2 - Pore Volume Determination 
Linear Model 


After a vacuum was drawn on the model, a plastic tube from a cali- 
brated cylinder containing an initially known volume of brine was connected 
to the bottom end of the model and brine was drawn up into the model due 
to the pressure difference between the model and the atmosphere. By inject- 
ing water from the bottom of the model, a more accurate pore volume and a 
more uniform water saturation could be achieved. The difference between 
the initial and final volumes of brine yielded the pore volume of the model. 
The porosity was calculated by dividing the pore volume by the bulk volume 
of the model. 


Two-Dimensional Model 


For the two-dimensional model, the determination of the pore vol- 
ume was more time-consuming than that for the linear model. First, the 
model was rotated so that the flange side faced down. Next, brine with a re- 
fractive index of 1.3446 was injected at the bottom of the model using the con- 
stant rate screw-type piston pump, while distilled water was being produced 
and collected at the top of the model. Brine injection was continued until the 
refractive index of the produced water reached 1.3446. At this time, the model 
was believed to be 100% brine saturated, and the injection was stopped. For 
each sample of water collected, its refractive index was measured using a re- 
fractometer to estimate a gradual change from water to brine solution. The 
refractive indices of the first and last water samples were plotted versus the 
percent of brine in solution, since it was believed that the first sample con- 
tained 0.0% brine and the last 100.0% brine (Figure 5.4). From this plot, 
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Figure 5.4 - Two-Dimensional Model Fraction of Brine in Solution. 
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Figure 5.5 - Two-Dimensional Model Pore Volume Determination. 
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knowing the refractive index of each water sample, the percent of brine in so- 
lution could be found. To determine the pore volume of the model, the per- 
cent of brine in solution was plotted versus the cumulative volume of water 
produced, and the area under the curve was divided into two equal portions 
(Figure 5.5 on the previous page). The pore volume was the cumulative vol- 
ume at which area A equalled area B. 


5.2.3 - Permeability Determination 


For both the linear and the two-dimensional models, the permeability 
was measured using the same approach. Note that after the pore volume de- 
termination, the model was brine saturated; thus, brine was used as the fluid 
to measure the permeability of the sand pack. The horizontal permeability of 
the model was measured by flowing brine through the model in a horizontal 
position, at a specific flow rate and pressure differential. A known volume of 
water was collected at a given time and pressure differential, and the perme- 
ability was determined using Darcy's linear flow equation for the linear 
model and Muskat’s equation”? for the two-dimensional model. 


5.2.4 - Oil Saturation 


For both physical models, the procedures to saturate the model sand 
pack with oil were similar. First, the model was inverted so that the injection 
port was facing up and the production port facing down, and the model pres- 
sure was brought to the experimental pressure by injecting brine into the 
model with the production back-pressure regulator (BPR) closed. Also, oil 
had to be pressurized to the experimental pressure by activating the constant 
rate Milroyal pump, with the inlet valve closed, until the oil pressure was at 
least a little higher than or equal to the model pressure. Then it was injected 
into the model at a very slow rate by fully opening the inlet valve. Right after 
oil breakthrough occurred, injection was stopped and the volume of brine 
produced was recorded. This volume of brine was used to predict the initial 
oil saturation, as follows. 
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For the linear model: 


Brine Volume Produced - Oil Volume in Chevron - Type Caps 


Soi = 
Pore Volume 


x 100% 


For the two-dimensional model: 


go Brine Volume Produced 100% 
oF Pore Volume 


At this time, the model was believed to be oil-saturated and ready for 
an experiment. 


5.2.5 - WAG Process, Post Waterflood, and Blowdown 


The same procedure was used to conduct an experiment in both physi- 
cal models. To start an experiment, the volumes of gas and brine to be in- 
jected had to be calculated first. The calculation procedure is shown below. 


Since it was found by previous researchers!” that a total gas slug size of 
20% HCPV and a water volume four times the gas volume, both of which 
were divided into ten equal slugs, were optimal, they were used in all exper- 
iments in this study. 


Total Gas Volume @ Experimental Conditions = 0.20 x HCPV_[cm’] 


It had to be converted to its equivalent volume at the meter (standard) 
conditions. 


Total Gas Volume @ Meter Conditions 


MD @ Experimental Conditions (mol/em’) 


= 0.20 x HCPV x [cm?] 


MD @ Meter Conditions (mol/em’] 
where, 
MD = molar density calculated using the Starling Equation of State!©. 


Total Brine Volume = 4 x Total Gas Volume @ Experimental 
Conditions [cm/%]. 
Finally, the total gas and brine volumes were divided into ten slugs 
each. 
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After the preliminary calculations had been completed, the model was 
prepared to start an experiment. With the model in the horizontal position, 
a gas slug was first injected then followed by a water slug until ten slug pairs 
of gas and water had been injected, and the WAG process ended. 


The WAG process was followed by the "post-WAG waterflood". This 
waterflood was carried out only when, after the WAG process, the producing 
WOR was still below 20:1. The model was flooded with brine until the WOR 
reached 20:1, or higher, when the waterflood was terminated. The "blow- 
down" was commenced by first closing the injection valve and then slowly 
lowering the pressure to atmospheric pressure by releasing the production 
BPR. Subsequently, the model was left for at least eight hours to make sure 
all gas was produced. At this time, the experiment was terminated. 


After the termination of the experiment, the model was opened and 
the sand pack was removed and discarded. The model, as well as the injec- 
tion and production ports, were cleaned first with Varsol, then toluene, and 
readied for the next experiment. The data collected were analyzed to deter- 
mine various parameters indicative of the overall performance of the exper- 
iment. A typical run took a total of two weeks. 


5.2.6 - Data Processing 


The experimental data were processed using a previously developed 


computer program? 


. The program was based on the material balance of oil, 
water, and carbon dioxide mixture. The volume of fluids injected was calcu- 
lated by this program. It also computed the water-oil ratios (WOR), gas-oil ra- 
tios (GOR), oil recovery, the total volume of oil produced, oil produced-fluid 
injected ratio (OPFIR), carbon dioxide retention and carbon dioxide required 


to produce a unit volume of oil. 


The carbon dioxide material balance used the Starling equation of 
state!°0. The equation of state is as follows: 
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3 
c 
+ir(l +p” )exp(-rp”) (5.1) 
where, 
p = pressure (MPa) 
T = Temperature (K) 
p = molar density (kmol/m?) 


The constants for carbon dioxide in SI units are: 


A, = 0.176976 B, = 0.024588 C, = 2.451876E04 
Dg = 1.883482E06 E, = 2.631556E04 R = 0.008314 
a = 0.009434 b = 0.003784 c = 1.4197888E03 
d = 0.055761 a = 0.0000961229 y = 0.006421 


Newton-Raphson's method was applied to the above equation to de- 
termine the molar densities of nitrogen-carbon dioxide mixtures of various 
compositions. According to Starling!, the above equation predicts experi- 
mental density data with an average error of less than 1.0%. 


5.3 - Diffusivity Experiments 


This section is divided into two sub-sections. The first presents the 
mathematical analysis leading to the determination of the diffusivity coeffi- 
cient, and the second discusses the experimental technique. 


5.3.1 - Mathematical Analysis 
The one-dimensional unsteady-state diffusivity equation is 


a; _ gt FC 
ot 3 Ox? 


(5.2) 


subject to the following initial condition: 


C,(x,t=0)=0 for O<x<o (5.3) 
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and the following two boundary conditions (assuming an infinite length): 


Ca xe ie 0) Ge (5.4) 

C;(x =2,t > 0) =0 (5.5) 
where, 

C;.= equilibrium concentration at the interface. 


The solution to equation (5.2) subject to the above conditions is!®! 


X 
: (5.6) 
24/Dit 


The total mass of the gas component i, m,, absorbed after time t can be 


C; (x,t) = C;,erfe 


defined as 


m, = AJ C;(x,t)dx (5.7) 
0 


where, 
A = diffusional cross-sectional area. 


Inserting equation (5.6) into equation (5.7) gives 


= M ir age scen! . 
m, rine Bs ; cae (5.8) 


The solution to equation (5.8) given by other investigators”® 3% 37,99 is 


te 
= JAG. Dit (5.9) 
T 


The parameter, m, is mass of gas component i injected into the diffu- 
sion cell after time t to maintain a constant pressure during the duration of a 
diffusion experiment, which can be determined directly from the observed 
volume of gas injected into the cell, which can be converted to mass by using 
an appropriate equation of state. For carbon dioxide, the Starling equation of 
state is the most accurate. The interfacial equilibrium concentration, C;,, in 
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kg/m? of solution, which can be measured experimentally by taking a sample 
at the top of the oil column. The diffusional cross-sectional area, A, can be 
measured directly. By plotting m, vs. vt, a straight line with a slope equal to 


2AC, | DE can be obtained. The diffusivity coefficient, Dj, can be determined 


from the slope of the straight line. 


5.3.2 - Experimental Measurement 


The method of measuring the carbon dioxide diffusivity in oils pre- 
sented here is quite similar to the one designed by Pomeroy et al.?8, except for 
the completion of the experiment. In Pomeroy et al.*°’s method, at the end of 
the diffusion period, the diffusion cell was shaken or continuously rotated to 
obtain equilibrium and complete mixing of the gas and oil. Then, a sample 
was taken, and the equilibrium concentration was determined. In this 
method, at the end of the time allowed for the diffusion to occur, a sample at 
the top of the oil column was taken, and the volume of gas released from the 
oil was used to determine the equilibrium concentration. 


The diffusion cell used in this study is shown in Figure 5.6. It consisted 
of a stainless steel cylinder fitted with two flanges. The internal cross-sec- 
tional area of the cell was 32.17 cm? and the length was 122.0 cm. The cell was 
always placed in the vertical position during packing and cleaning, as well as 
during the actual experiments. The top flange was connected to a high pres- 
sure carbon dioxide cylinder. A Heise pressure gauge was also connected to 
the top flange to measure the pressure inside the cell during the experiment. 
The bottom flange was equipped with a two-way high pressure valve which 
permitted collection of oil samples for determining the concentration of dif- 
fused carbon dioxide in the oil at the end of the experiment. 


The procedure for conducting a diffusion experiment was as follows. 
First, the diffusion cell was evacuated for six hours using a vacuum pump 
connected to the top flange of the cell. Next, while the cell was still being 
evacuated, a plastic tube from a calibrated cylinder containing an initially 
known volume of oil was connected to the bottom valve, and oil was drawn 
up into the cell to obtain an oil column of at least 20 cm, and the bottom 
valve was closed. Carbon dioxide at the desired experimental pressure was 
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injected into the cell at the top flange, and the pressure was kept constant 
during the course of the experiment by injecting carbon dioxide. The volume 
of car bon dioxide injected and the time it was injected were recorded. As the 
injection of carbon dioxide continued, a straight line was obtained on the plot 
of m;, vs. vt. At this time, a sample at the top of the oil column was taken. 
The volume of gas released was trapped and measured. The experiment was 
then terminated. The cell was opened and cleaned to prepare for the next 
experiment. 


A sample plot of m, vs. Vt for one of the diffusivity experiments is 
provided in Figure 5.7, on the next page. The lower part of the curve near the 
origin where the data points are scattered indicates the period during which 
the interfacial equilibrium concentration is becoming established. The upper 
part of the curve where the curve rises rapidly indicates the diffusing gas has 
reached the bottom of the diffusion cell. When this has occurred, it means 
that the infinite boundary condition assumption is no longer tenable. The 
middle section (between A and B), where a straight line is established indi- 
cates the diffusion period which may be used to determine the diffusivity co- 
efficient. 


It should be noted that when carbon dioxide diffuses in oil, it causes the 
oil to swell. In this study, the effect of swelling was included in the calcula- 
tion of the equilibrium concentration at the interface by using Welker and 
Dunlop’s correlation!’. 
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Figure 5.7 - Mass of Carbon Dioxide Injected into Diffusion Cell vs. Square 
Root of Time. 
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6 - DISCUSSION OF EXPERIMENTAL RESULTS 


This chapter consists of two parts. The first part deals with the presen- 
tation and discussion of the work done to investigate the diffusivity of a gas 
phase into a liquid phase. The second part contains the discussion of the dis- 
placement experiments — isothermal and non-isothermal. 


I - Diffusivity of Carbon Dioxide and Methane 


Tables 6.1 and 6.2, respectively, present the 76 measurements made to 
investigate the diffusivities of carbon dioxide and methane into different oils. 
The raw data were plotted and are included in Appendix D. The conditions at 
which the measurements were made are given in Tables 6.1 and 6.2. 


6.1 - Effect of Pressure 


Figure 6.1 shows a log-log plot of the diffusivities of carbon dioxide 
versus pressure for an 1842 mPa.s oil at 294.15K. Figures 6.2 and 6.3, respec- 
tively, present the logarithm of the diffusivity of carbon dioxide in 1842 and 
3607 mPa.s oils vs. the logarithm of pressure at the temperature conditions at 
which carbon dioxide is gaseous. It is apparent in these two figures that as 
long as carbon dioxide is in a gas phase, its diffusivity increases linearly with 
increasing pressure. In Figure 6.1, there are two distinct regions: gaseous and 
liquid carbon dioxide diffusion. First consider gaseous carbon dioxide diffu- 
sion consisting of Experiments 3 to 7 conducted at 0.69 to 5.17 MPa. As is 
shown by the gaseous carbon dioxide diffusion line, the diffusivity of carbon 
dioxide into oil increased with increasing pressure, under the experimental 
conditions. Now, consider the liquid carbon dioxide diffusion region, consist- 
ing of Experiments 8 to 11 performed at pressures from 6.89 to 10.34 MPa. 
Clearly, liquid carbon dioxide diffusion remained unchanged and was much 
lower than gaseous carbon dioxide diffusion. The reasons can be as follows. 
When carbon dioxide is in a gas phase, molecular forces are weak, thus allow- 
ing carbon dioxide molecules to move quickly and freely into the oil. Hence, 
high carbon dioxide diffusion can result. On the other hand, at 6.89 to 10.34 
MPa and 294.15°K, carbon dioxide is in a liquid phase, which means that the 
molecular speed of the carbon dioxide molecules is low due to the strong 


66 


i ; - | a4 é 
“warwig af) Ai alauia sent et iter oor! Ye 
Hon se yeh il) aieghes ray 6) enh frown aa 
: 4 oe el “a aah sae Sane r) ay Sie ie Sta 


ms (acter) rew Gan larrre 


rob ot . yp reaanien AY Ge peer alert en | hes £8. valde? a 

pu io. oly oe wheat hela ebixiilb ries Mh al hag ARS acy oie oe: aadeorn - : ; 

te ROW TD! " Ci Veo *. ov tuhitsni ste hin ne eich wer of ~ 
beet Lith epee 4 - 4649 Ghats «me eonmcrtoniaetre ont f sie 


a: _—? 


. eure 
7.1 is = 
enor te cer wit ra % 
enlay. ive , oC: Ons te: tale yo’ Gel « awoda ia ® 
P50 bas are ee0 ta lio a. le CORE Gee teem 
iY on sorb A ane: > opbecteg: dy wilt jo whicaget or) soar i 


fe anh a@ihine of psvtwead 4 free Pier: i ocrlraget oath ine alla pe 
és tod! gehogf hey eiteerigh Ait Slrovety a shies teed 

y »/ ase, soya? Stared, wil wh; agg.2 nba itomib 0 ott 
be ced aver ua dunt Suvi oe Sed! tA og sé “ohwaty 3 uneeial 

my, pds ayieee “ taimavies teri shler > abewdl rrr ilk af sn a 
jiiedy BOLO FE Ses eg teat hota © ' oo |) citearme aniy- = as a 
pistes er wet eh LY: alo eat Meas ** Nis pial) nec apron 

AE ee bo wire ieee BR San PD tt itty “= : 
Low?.200 erat: Sets aR dens mare whol + att 
ai SL.01 ul PRA ane i went MBAt YY vom cP bes 


Maki: swe seul Iie ie ‘hs Pasig” fh smite A sei: 

ay IO! ort a) ae ; ae ~caraphag 
LS ita tent “ . 

Re ee as ery sien toh 

ed ae. ph ol Ue #'t toad 


ts ; ub Sach: opsopny 
s a sechanc awn 


Table 6.1 - Diffusivity of Carbon Dioxide. 


Experiment 
No. 


Oil Viscosity | Temperature | Pressure 
at 21°C and 
0.101 MPa 

(mPa.s) 


603 6.89 
1058 6.89 


1842 21.0 
1842 21.0 1.38 
1842 21.0 2.07 
1842 21.0 3.45 
1842 210 57 


1842 6.89 
1842 6.89 
1842 8.62 
1842 10.34 
1842 

1842 2.07 
1842 3.45 


1842 57.1 6.89 
1842 Dial 10.34 
1842 O71 15a 
1842 yt 17.24 
1842 
1842 65.4 2.07 
1842 65.4 3.45 
1842 65.4 6.89 
1842 65.4 10.34 
1842 65.4 1379 
1842 65.4 17.24 


1842 


1842 725 2.07 


1842 3.45 


Diffusivity 
(m?/s) 


2.5452E-08 
7 2974509 
1.2677E-09 
1.6529E-09 
2.0635E-09 
2.4311E-09 
2.7448E-09 
1.0025E-09 
1.0693E-09 
1.1857E-09 
1.0811E-09 
2.3648E-09 
3.8234E-09 
4.6159E-09 
6.1387E-09 
7.2715E-09 
8.4097E-09 
9.0767E-09 
2.8093E-09 
4.4365E-09 
5.5934E-09 
7.9272E-09 
8.3963E-09 
9.5944E-09 
1.1040E-08 
3.5671E-09 
5.4081E-09 
6.7305E-09 
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Table 6.1 - Cont’d. 


Experiment 
No. 


Oil Viscosity | Temperature | Pressure | Diffusivity 

at 21°C and 

0.101 MPa 
(mPa.s) 


1842 6.89 | 9.1780E-09 
1842 79.3 10.34 | 1.1575E-08 


WIN 
=) ||e) 


1842 79.3 13.79 
1842 79.3 17.24 
1842 93.2 
1842 93.2 2.07 
1842 3.45 
1842 93.2 6.89 
1842 93.2 10.34 


WwW | Q | Q 
QIN rR 


Go 


4 


16%) 
Oo 


OQ) 
N 


1S) Q 
00 ON 


1842 93.2 13.79 
1842 17.24 
3607 

3607 2.07 
|___ 3607 3.45 

3607 
3607 57.1 2.07 
3607 3.45 
3607 6.89 
3607 10.34 
3607 57.1 13.79 
3607 17.24 
3607 
3607 2.07 
3607 3.45 
3607 6.89 
3607 10.34 
3607 13.79 
3607 17.24 | 3.5797E-09 
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Table 6.1 - Cont’d. 


2.07, 


3607 oo. 
3607 93.2 


70 3607 17.24 | 5.0193E-09 
15402 4.9337E-10 


Table 6.2 - Diffusivity of Methane. 


Experiment 
No. 


Oil Viscosity | Temperature | Pressure | Diffusivity 

at 21°C and 

0.101 MPa 
(mPa.s) 


Pmt] oo | 3.151811 | 
: 
; 


13.79 | 4.1067E-10 
17.24 | 4.5441E-10 
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Figure 6.1 - Diffusivities of Gaseous and Liquid Carbon Dioxide. 
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Figure 6.2 - Diffusivity of CO, in an 1842 mPa.s Oil as a Function of Pressure. 
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molecular forces. As a result, the diffusivity of carbon dioxide is reduced. 
This is why low diffusion coefficients were obtained at 6.89 to 1034 MPa and 
294.15 K, at which carbon dioxide was liquid. 


This finding of low liquid carbon dioxide diffusion helps explain the 
low oil recoveries reported by Rojas!, who injected liquid carbon dioxide 
rather than gaseous carbon dioxide in some of his immiscible carbon dioxide 
WAG displacement experiments. 


6.2 - Effect of Increasing Temperature 


One of the objectives of this research was to investigate the effect of in- 
creasing temperature on the diffusivity of carbon dioxide. Figures 6.4 and 6.5 
depict the effect of temperature on the diffusivity of carbon dioxide into 1842 
and 3607 mPa.s oils at different temperatures. In these figures, the log-log plot 
of the diffusion coefficient in m?/s versus temperature in K is presented. A 
straight line was obtained, indicating that the diffusivity of carbon dioxide in- 
creases exponentially with increasing temperature. The reason for an in- 
crease in carbon dioxide diffusion with temperature increase is that increasing 
temperature causes an increase in the kinetic energy which accelerates the 
movement of the carbon dioxide molecules, consequently leading to an incr- 
ease in carbon dioxide diffusion. Also, at an elevated temperature, oil 
molecules are more spread out, the “holes” among the oil molecules are big- 
ger, the intermolecular forces between the oil molecules are weaker, and the 
viscosity of the oil is less, thus allowing the carbon dioxide molecules to dif- 
fuse through more easily. 


6.3 - Effect of Oil Viscosity 


Five experiments were performed using oils having different viscosi- 
ties to study the influence of oil viscosity on the diffusivity of carbon dioxide. 
They were Experiments 1, 2, 15, 46, and 71. The viscosities of the oils used in 
these five experiments were 603, 1058, 1848, 3607 and 15402 mPa.s, respec- 
tively. The pressure and temperature at which these five experiments were 
conducted were 6.89 MPa and 330.26°K, respectively. The results are shown in 
Figure 6.6 as the logarithm of the diffusivity versus the logarithm of the oil 
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Figure 6.4 - Diffusivity of CO, in an 1842 mPa.s Oil as a Function of Temperature. 
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Figure 6.6 - Diffusivity of Carbon Dioxide as a Function of Oil Viscosity. 
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Figure 6.7 - Diffusivity of Carbon Dioxide as a Function of Oil 
Molecular Weight. 
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viscosity. It is noted that the diffusivity-oil viscosity relationship is exponen- 
tial and that the diffusivity decreases with increasing oil viscosity. This is be- 
cause lower viscosity oils have weaker molecular bonds. 


6.4 - Effect of Oil Molecular Weight 


Besides oil viscosity, oil molecular weight was also found to have an 
effect on the diffusivity of carbon dioxide. Figure 6.7 (on the previous page) 
shows the logarithm of the carbon dioxide diffusivity vs. the logarithm of the 
oil molecular weight. The diffusivity coefficients plotted in this figure are 
those of Experiments 1, 2, 15, 46, and 71, which were discussed in Section 6.3. 
As shown, the diffusivity of carbon dioxide exponentially decreases with in- 
creasing oil molecular weight. 


6.5 - Diffusivity of Methane 


Besides studying the diffusivity of carbon dioxide, the diffusivity of 
methane into oil was also studied. Five experiments (Experiments 72 to 76) 
were performed using methane as the diffusing gas at 330.26K and 0.69 to 
17.24 MPa. The results are given in Table 6.2 and plotted in Figure 6.8. Like 
carbon dioxide diffusivity, the diffusivity of methane into oil increased with 
increasing pressure. In Figure 6.8, the diffusivities of carbon dioxide into an 
oil of identical viscosity (i.e., 1846 mPa.s) are also shown. Comparing the two 
straight lines in the figure reveals that more carbon dioxide than methane 
could diffuse into the oil under the experimental conditions. The reason is as 
follows. According to the Stokes-Einstein equation, for diffusing molecules 
large with respect to solvent molecules, the rate of diffusion depends in- 
versely on the radius of the diffusing molecule (i.e.,Qj, @ 1/V'/%)!. Ross 
and Hildebrand!%, Nakanishi, Voigt and Hildebrand!°4 showed that with 
dissolved gas molecules the diffusivity mainly depends on the cross-sectional 
area (V2/3) of the diffusing molecules. Based on these two theories, it can be 
said that diffusivity varies inversely with the molecular volume to the third 
power (Q)j,a 1/V*/%). Under experimental conditions, the molecular vol- 


umes of methane (eg., 763.75 cm?/mol at 3.45 MPa and 57.1°C) were much 
higher than those of carbon dioxide (eg., 686.20 cm?/mol at 3.45 MPa and 
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Figure 6.8 - Diffusivity of Methane in an 1842 mPa.s Oil. 
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57.1°C). Therefore, the diffusivity of methane was found to be lower than 
that of carbon dioxide. 


6.6 - Correlations of the Diffusion Coefficients 


The technique of correlating the experimental data employed in this 
work is the one first proposed by Wilke and Chang!, which is widely used. 
The main idea of this technique is to collapse all data points on one single 
straight line. More details on this technique can be found in Reference 105. 
From the discussions presented in Sections 6.1 to 6.4, it is obvious that the 
diffusivity of carbon dioxide in oils depends on pressure, temperature, oil 
viscosity, and oil molecular weight. As a result, an empirical equation which 
is similar to the one proposed by McManamey and Wollen.*!, that has a form 
as shown below, is adopted. 


Deo, oi = op" T'LIMWe. (6.1) 


where @, a, B, y, and € are the correlation parameters which are to be deter- 
mined. 


By linear regression analysis, the numerical values of o, a, B, y, and € 
were respectively found to be 5.7582 x 1092 0.3734, 5.6110, -0.6606, and -4.4851. 
Thus the final form of the correlation is | 


0.3734-75.6110 
* (9) |) 4h 
Deo, -oi = 95-7582 x 10° “yep a aes: (6.2) 
sie Uo MW, 


where, p=MPa, 
Wels, 
lo = mPa.s, 
MW, = g/mole. 


Figures 6.9 to 6.10 compare the experimental and calculated values. 
This correlation was found to predict the diffusivity of carbon dioxide in 
heavy oils with an average error of less than 6%. 
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Figure 6.9 - Comparison of Correlated and Measured Values for 
Kla-Da-Ing Oil. 
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Figure 6.10 - Comparison of Correlated and Measured Values for 
Della-Bell Oil. 


78 


« —— — 


————— — 


Mie TY 


acter ely > - tbr A ‘ Te te P| , 


——le) Pe —— ed 


7 
~~ es 
+ 
= 
ee eT as ; Fy 
~~... 7 
—\ Ui Vn 4 hs oT] 
~ a” 4 ~ a, i 
« = a] ~ 
— i “Tike, i PA a, : 
y'vs? se < he ma 


eA 
( Anm™ nla a: 


vo an 


II - Displacement Experiments 


The second part of this chapter presents a discussion of the displace- 
ment experiments conducted in previous and present studies. The first sec- 
tion presents correlations of the experimental data (included in Table 6.3 on 
page 80) reported by the former investigators at the University of Alberta! +>’ 
with the dimensionless groups. The second section presents a discussion of 
the experiments made in this study. Most of the experiments conducted in 
this study and previous studies!4”-” were completed So; (initial oil saturation) 
instead of Sor (residual oil saturation) because the process is apt to be em- 
ployed in reservoirs where primary recoveries are very low and waterflood- 
ing is inefficient. Four experiments were performed to investigate the possi- 
bility of using a carbonated waterflood in place of an immiscible carbon diox- 
ide WAG flood. Table 6.4 contains the results of these four carbonated water- 
flood experiments. Fifty experiments were made to study the gravity segrega- 
tion of the injected fluids. Seven experiments were conducted to observe the 
performance of the immiscible carbon dioxide WAG process under non- 
isothermal conditions. Two physical models, one linear and one two-dimen- 
sional, were used to conduct the experiments. The results reported in this 
study are summarized in Tables 6.5 and 6.6. The prefixes "V" and "H" in the 
run No. refer to vertical and horizontal, respectively. A typical run took a to- 
tal of two weeks or longer, depending on the type of run. For non-isothermal 
runs the results of which are summarized in Table 6.8, the model had to be 
heated to the experimental temperature and pressure conditions for at least 5 
to 7 days or even longer to make sure the temperature in the model was uni- 
form, as shown by a constant temperature on the temperature controller. 


Figure 6.11 provides a pictorial view of the different experiments con- 
ducted in the linear and two-dimensional models. The tabulated results of 
the experiments can be found in Appendix E. Appendix F contains the pro- 
duction histories of all experiments conducted. 
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Table 6.6 - Summary of Horizontal WAG Displacement Experiments in a Quarter of a 5-Spot System (Cont’d). 
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Figure 6.11 - Map of experiments conducted in this research. 
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6.7 - Correlation of Displacement Results 


This section presents the correlations of the experimental displacement 
results reported by Rojas!, Zhu‘, Dyer’, Prosper’, and in this study. The di- 
mensionless groups used in the correlations are those previously derived by 
Rojas!. Note that the groups derived by Rojas! are a sub-set of the two sets de- 
rived in this study. 


6.7.1 - Single Slug and WAG Correlations 
Figure 6.12 presents the correlations of the overall average gas-oil ra- 


DpH, 
skAp,, 


tios!4 with the dimensionless group (which is the ratio of viscous 


forces to gravitational forces) for single slug immiscible carbon dioxide dis- 
placements and immiscible carbon dioxide WAG displacements in a quarter 
of a 5-spot flood pattern. Let us first consider the upper line representing the 
overall average gas-oil ratios for single slug carbon dioxide injection experi- 
ments (Runs 6R-10R! and 12Z*) conducted using injection velocities ranging 
from 0.18 to 2.9 m/d, and a total carbon dioxide slug size of 20% HCPV. In 
these runs, carbon dioxide was continuously injected until 20% HCPV had 
been injected; then, water was injected to displace the carbon-dioxide swollen 
oil until the water-oil ratio reached about 20 to 1, which is the cut-off water- 
oil ratio. As is shown in the figure, the overall average gas-oil ratio linearly 


nN) 
increases from 22 to 508 sm?/sm? as ee increases from 0.0096 to 0.23, or the 
gkAp,, 


injection velocity increases from 0.18 to 2.9 m/d. In other words, there is an 
exponential relationship between the overall average gas-oil ratio and the ra- 
tio of viscous forces to gravitational forces for single slug immiscible carbon 
dioxide injection. It is clearly shown by the straight line that a higher gas in- 
jection rate results in a higher overall average gas-oil ratio. This is basically 
due to a higher mobility ratio at a higher injection rate. It should also be 
noted that when carbon dioxide is injected at a high rate, it does not have 
enough time to dissolve as much as it can to develop a complete phase equi- 
librium with the oil phase. Therefore, at a higher rate, a larger volume of 
carbon dioxide is produced during the flood. In contrast to this, a lower injec- 
tion rate will give carbon dioxide more time to dissolve and to develop a bet- 
ter phase equilibrium with the oil. As is indicated in the figure, at a viscous- 
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gravitational force ratio equal to 0.0096, the overall average gas-oil ratio is 22 
sm?/sm°, which is about one-twenty fifth of that at 0.23, for single slug im- 
miscible carbon dioxide injection. 


On the bottom of the same figure, there is a dashed line which repre- 
sents the overall average gas-oil ratios for experiments (Runs 16R!, 23R!, 
13Z4, 16Z4, 24Z4 and 25Z*) carried out utilizing the immiscible carbon dioxide 
WAG injection method. Similar to the single slug carbon dioxide injection, 
the overall average gas-oil ratio for the carbon dioxide WAG process rises to a 
higher value at a higher viscous-to-gravitational forces ratio. The inserted 
picture shows this more clearly. Observing the two gas-oil ratio curves, one 
sees that much less gas, even at a higher viscous to gravitational force ratio, is 
produced during the carbon dioxide WAG injection than during the single 
slug carbon dioxide injection. This means that injecting carbon dioxide in the 
WAG mode reduces the mobility of the gas phase, which consequently leads 
to an increased amount of carbon dioxide going into solution in the oil. 


6.7.2 - WAG Ratio Correlations 


It has been shown that injecting carbon dioxide in the WAG fashion 
yields better mobility control on the gas phase than injecting it in the single 
slug mode!. One can ask at what WAG ratio the carbon dioxide WAG injec- 
tion will give its best performance. Referring to Figure 6.13, where the over- 
all average gas-oil ratios of Runs 14R-18R! conducted using the same carbon 
dioxide slug volume but different WAG ratios! are plotted versus the prod- 
ucts of the viscous-gravitational forces ratio and WAG ratio, the dashed con- 
cave downward curve represents the overall average gas-oil ratio curve while 
the solid concave upward curve represents the oil recovery. The WAG ratios 
employed in the experiments were varied from 1:1 to 6:1. In order to make a 
fair comparison, the same injection velocity and oil were used in these runs. 
It is shown by the overall average gas-oil ratio curve that the volume of pro- 
duced gas decreases as the WAG ratio increases until the WAG ratio reaches 
4:1, where the gas-oil ratio is minimum. Then it rises back up when the 
WAG ratio increases to 5:1 then to 6:1. The reasons may be most likely as fol- 
lows. When a WAG ratio lower than 4:1 is utilized, the slug volume of water 
is not enough to reduce the mobility of the carbon dioxide gas phase. When a 
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Figure 6.13 - WAG Ratio Correlation (Data from Ref. 1). 
20 


3 


Legends 
© 2% HCPV (GTD1’) 


o 10% HCPV (GTD3’) 
¢ 20% HCPV (H2D8) 


= 30% HCPV (GTD4’) 
a 61% HCPV (TD1°) 


3 


Overall Average Gas-Oil Ratio, sm /sm 
a 


2-D Flood at 2.5 MPa 
4:1 WAG 


— 
i>) 


U1 


: 2 3 Pe A Oe 2 ad SEG 07 8:9 2 aes? aes 
0.001 i 

Delo Voo, 

gkAPp,, 


Figure 6.14 - Total Carbon Dioxide Slug Size Correlation (Data from Refs. 5 & 7). 
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WAG ratio higher than 4:1 is used, the water saturation in the flushed zone is 
higher while the corresponding oil saturation is lower, compared to when a 
4:1 WAG ratio is used. This seems to affect the mass transfer rate of carbon 
dioxide to oil. It is known that the mass transfer rate of carbon dioxide to oil 
in porous media depends on the oil saturation, i.e., the lower the oil satura- 
tion, the lower the rate of mass transfer of carbon dioxide to oil. Experimental 
work done by Denoyelle and Bardon’ verifies this. Thus, the lower is the oil 
saturation, the lower is the interfacial area available for mass transfer between 
carbon dioxide and oil. This plays a very important role in mass transfer. It is 
because a lower/higher interfacial area will lead to a lower/higher mass trans- 
fer rate. As a result of this, when the next carbon dioxide slug is injected fol- 
lowing the preceding water slug, more carbon dioxide will go into the water 
phase, or remain in the gas phase, than goes into the oil phase, thus resulting 
in the production of large gas volumes. 


In terms of displacement efficiency, the 4:1 WAG ratio is the most effec- 
tive. As is shown by the oil recovery curve, the highest recovery was ob- 
tained with the 4:1 WAG ratio. 


Summarizing, the use of a 4:1 WAG ratio in the immiscible carbon 
dioxide displacement process yields the most effective mobility control of the 
gas phase and the highest displacement efficiency. 


6.7.3 - Total Carbon Dioxide Slug Size Correlations 


Five experiments (Runs GTD1’, GTD3’, GTD4’, TD1°, and H2D3) were 
performed using slug sizes of 2, 10, 20, 30, and 61% HCPV at 2.5 MPa and 21°C. 
These runs were conducted using a 4:1 WAG ratio. Also, the same injection 
rate and oil were utilized in these runs. The results of these runs were corre- 
lated with the product of the ratio of viscous-to-gravitational forces and slug 
size!” and are shown in Figures 6.14 and 6.15. Figure 6.14 provides a correla- 
tion of the overall average gas-oil ratios with the product of the ratio of vis- 
cous-to-gravitational forces and slug size. As indicated, there is an exponen- 
tial relationship between the gas-oil ratio and the dimensionless group 
B eH Vco, 


, and the gas-oil ratio increases when the slug size increases. This is 
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Figure 6.15 - Recovery-Slug Volume Correlation (Data from Refs. 5 & 7). 
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Figure 6.16 - Total Slug Size Correlation for Low Pressure Operations 
(Data from Ref. 7). 
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because, under a given pressure and temperature condition only a certain 
volume of carbon dioxide can dissolve in the oil and the excess remains in 
the gas phase, which will then bypass the oil. The higher is the slug size, the 
higher are therefore the gas volume and the producing gas-oil ratio. 


The correlation of recovery with the same dimensionless group as 
mentioned above was also made and is shown in Figure 6.15. The recovery 
increases with a larger slug size and reaches its maximum at a slug size equal 
to 20% HCPV, then flattens out when the slug size is greater than 20% HCPV. 
The same correlations were performed as well for Runs ENE DE, {RE Be 
GTD8’, GTD10’, H2D26 and H2D31-33 conducted at 1.0 MPa. Figure 6.16 pre- 
sents a correlation of the overall average gas-oil ratios with the product of the 
ratio of viscous-to-gravitational forces and slug size for these runs. The total 
carbon dioxide slug sizes used in these runs were varied from 10 to 30% 
HCPV, in an increment of 10% HCPV. The same feature observed in Figure 
6.14 can also be observed in Figure 6.16. As is shown in Figure 6.16, at 1.0 
MPa, increasing the volume of carbon dioxide used in the immiscible carbon 
dioxide WAG process increases the producing GOR. The explanation men- 
tioned the preceding paragraph applies here also. 


Figure 6.17 provides the correlation of recovery with the same dimen- 
sionless group for Runs GTD5’, GTD7’, GTD8’, GTD10’, H2D26 and H2D31-33 
conducted at 1.0 MPa. The figure shows three different straight lines for the 
oil recoveries of runs utilizing total carbon dioxide slug sizes of 10, 20, and 
30% HCPV. It is indicated in Figure 6.17 that the recovery at 1.0 MPa is 
optimal when the volume of carbon dioxide utilized in the immiscible WAG 
process is 20% HCPV. 


In addition to correlating the results of the experiments conducted uti- 
lizing a quarter of a 5-spot flood pattern, correlations were also done for those 
conducted utilizing linear core floods. Figure 6.18 presents a correlation of 
the overall average gas-oil ratios with the product of the ratio of viscous-to- 
gravitational forces and slug size for Runs LC16°, LC17°, LC30-33 and LC36’ 
conducted at 2.5 MPa, in the linear model. The total carbon dioxide used in 
these runs were 5, 20, 40, 64, and 89% HCPV. In these experiments, an oil 
with a viscosity of 1055 mPa.s was used; except in a few runs, oils with viscosi- 
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Figure 6.17 - Recovery-Slug Volume Correlation for Low Pressure 
Operations (Data from Ref. 7). 
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Figure 6.18 - Total Slug Volume Correlation for a Linear Coreflood 
(Data from Refs. 5 & 7). 
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Figure 6.19 -Recovery-Slug Volume Correlation for a Linear Coreflood 
(Data from Refs. 5 & 7). 
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ties of 1230 mPa.s and 1130 mPa.s were used. Similar to what already ob- 
served, the producing GOR increases with increasing the total carbon dioxide 
slug size for the experiments conducted in the linear model. The correlation 
of recovery with the same dimensionless group for linear coreflood experi- 
ments performed at 2.5 MPa is shown in Figure 6.19. The curve in the figure 
reveals that there is no increase in oil recovery when the slug size greater 
than 20% HCPV is used. This is similar to the observation made on Figure 
Gi: 


In conclusion, the immiscible carbon dioxide process has its best per- 
formance when it is employed in the WAG mode, with a WAG ratio of 4 to 1 
and a total slug size equal to 20% HCPV. 
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II.1 - Isothermal Displacement Experiments 
6.8 - Carbonated Waterflood 


As has already been mentioned in Chapter 3, the idea of flooding an oil 
reservoir with carbonated water is not new. The method was first tested in 
the laboratory and tried in the field in the 50’s and 60’s. In this study, several 
experiments were made by injecting carbonated water, instead of injecting 
carbon dioxide alternately with water. The experimental results for these 
runs are summarized in Table 6.4 on page 82. 


6.8.1 - Carbonated Waterflood vs. Immiscible WAG flood 


Run CWFI1 was conducted by injecting carbonated water at 2.5 MPa and 
21°C, in a linear model. The carbonated water was prepared by mixing brine 
with carbon dioxide at experimental pressure and temperature conditions 
until equilibrium was reached, which could be ascertained by observing the 
constant pressure on the pressure gauge. Before the carbonated water was 
injected, the solubility of carbon dioxide in brine was measured to be 16.25 
sm?/sm? at 2.5 MPa and 21°C. The pH of the carbonic acid formed due to the 
chemical reaction taking place between carbon dioxide and brine was also 
measured using a pH meter. It was measured to be 5.1. 


The experiment was conducted by continuously injecting carbonated 
brine until the producing water-oil ratio reached 20:1; then injection was 
stopped to start blowdown. Figures 6.20 depicts the production history of this 
run. 


The oil recoveries of this run at each phase of production were 51.71% 
of oil recovered in the carbonated brine flood phase and 4.82% in the 
blowdown phase, thus giving a total recovery of 56.53%. The GOR curve 
shows that during the early part of the flood no carbon dioxide gas was 
produced, which is very good in terms of controlling the mobility of the gas 
phase. After gas breakthrough occurred at 0.72 PV, the production of gas 
climbed and reached its maximum value at 228 sm?/sm° at the end of the 
flood, which is at 1.73 PV as indicated in the figure. This very high producing 
GOR is undesirable in any enhanced oil recovery method where the injected 
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Figure 6.20 - Production History of Run CWF1. 
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gas plays a dominant role. 


In order to compare the performance of this recovery method, the GOR 
of this run (i.e., Run CWF1) was plotted on the same graph as that of Run 
LC31’, which was done using an immiscible carbon dioxide WAG flood uti- 
lizing a 4:1 WAG ratio and a total slug size of 20% HCPV at 2.5 MPa. Figure 
6.21 shows the comparison. Examining the two GOR curves reveals some in- 
teresting features. At early time, the same volumes of gas were produced in 
both runs; but after gas breakthrough, much more gas was produced in Run 
CWF1 than in Run LC31’. Gas breakthrough occurred earlier in Run CWF1 
than in Run LC31’, at 0.72 PV in Run CWF1 and at 1.02 PV in Run LC31’. As 
indicated in the figure, at 1.0 PV both runs seemed to have nearly the same 
GOR; but after 1.0 PV, the GOR in Run CWF1 steeply and continuously rose 
while it gradually decreased in Run LC31’. On average, after gas break- 
through the GOR of Run CWF1 was about 10 times higher than that of Run 
LC31’. It can also be seen easily that Run CWF1 had a shorter flood life than 
Run LC31’. 


Figure 6.22 depicts the oil recoveries of Runs CWF1, LC31’, and LC34’, 
which was conducted by injecting brine only. It is clearly shown in the figure 
that the recovery curve for Run CWF1 nearly lies on that for Run LC31’ and 
that both runs have nearly an identical recovery. The recoveries for Runs 
CWF1 and LC31’ were, respectively, 56.53 and 57.1%. In terms of oil recovery, 
the two runs are comparable; but, in terms of carbon dioxide requirement, 
which is the volume of carbon dioxide required to produce one cubic meter of 
oil at standard conditions, the two runs are not comparable. In Run CWF'1, it 
required 63.6 sm? to produce one standard cubic meter of oil while it required 
only 10.29 sm? to produce the same volume of oil in Run LC31’. This indi- 
cates that in order to recover as much oil as the immiscible WAG process, the 
amount of carbon dioxide required for the carbonated brine flood should be 6 
times higher. The very high GOR and carbon dioxide requirement in Run 
CWF1 is due to the slow mass transfer process from carbon dioxide to oil. In 
Run CWF1, carbon dioxide had to diffuse out of the non-diffusing liquid wa- 
ter phase before contacting the oil, whereas in Run LC31’ carbon dioxide was 
in direct contact with the oil. 
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Figure 6.22 - Comparison of Recoveries of Runs CWF1, LC31’, and LC34’. 
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Figure 6.22 also shows a comparison of Run CWF1 with Run LC347 (a 
waterflood run). Examining the two recovery curves reveals that both runs 
had nearly the same flood life, but Run CWFI1 had a higher recovery than 
Run LC34’, which had a recovery of 41.45%. Comparing the two recoveries, 
approximately 36% more oil was produced in Run CWF1 than in Run LC34’. 


Summarizing, the carbonated brine flood method is as effective as the 
immiscible WAG method when the volume of carbon dioxide used in the 
former is 6 times higher than that in the latter. 


6.8.2 - Effect of Carbon Dioxide Requirement 


To study the effect of the carbon dioxide requirement on flooding an oil 
reservoir with carbonated brine, Run CWF2 was performed in a scaled model, 
utilizing a quarter of a five-spot pattern. In this run, the carbonated brine was 
prepared by mixing 20% HCPV of carbon dioxide with brine at a 4 to 1 ratio 
and at 1.0 MPa and 21°C. After the carbonated brine was injected, water was 
injected to bring the WOR to the limiting WOR. The production history of 
this run is depicted in Figure 6.23. 


Figure 6.24 depicts a comparison of Run CWF2 with Run GTD6’, 
which was conducted employing the immiscible carbon dioxide WAG process 
with a total carbon dioxide slug size equal to 20% HCPV and a 4:1 WAG ratio. 
It is shown in Figure 6.24 that from 0 to 1.04 PV, the GOR of Run CWF2 is 
higher than that of Run GTD6’, that after 1.04 PV the reverse is true, and that 
gas breakthrough in Run CWF2 occurred at 0.28 PV, which is earlier than that 
in Run GTD6’. Moreover, Run CWF2 had a shorter flood life than Run 
GTD6’ because the limiting WOR was reached earlier in Run CWF2, as 
shown in the figure. 


Figure 6.25 presents the recoveries of the three runs CWF2, GTD6’ and 
21a!3. Run 21a! was done by injecting brine only. The figure clearly demon- 
strates that at 1.0 PV 37% of the oil was recovered by the carbon dioxide WAG 
flood, 32% by carbonated brine flood, and 30% by waterflood. This shows that 
up to 1.0 PV, a carbonated brine flood is no better than a waterflood. It is also 
shown in the figure that, up to nearly 0.6 PV, the recovery curve of Run 
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Figure 6.23 - Production History of Run CWF2. 
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Figure 6.25 - Comparison of Recoveries of Runs CWF2, GTD6’, and 21a’. 
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CWF2 almost overlaps that of Run 21a!3. Examining the three curves re- 
veals that during the duration of the flood more oil was recovered in Run 
GTD6’ than in Runs CWF2 and 21a!. The total recoveries for the three 
Runs CWF2, GTD6’ and 21a! were respectively 51.3, 39.2 and 32.4%, which 
indicates that the WAG flood recovered 12.1% more oil than the carbonated 
brine flood and that the oil recovered by the carbonated brine flood followed 
by a waterflood up to 1.0 PV is almost the same as that by a waterflood alone. 


In short, a carbonated waterflood followed by a waterflood is not com- 
parable to an immiscible WAG flood when the same amount of carbon diox- 
ide is used. The former is no more effective than a waterflood. 


6.9 - Gravity Segregation 


In this section, the results of the experiments conducted to investigate 
the effect of gravity segregation on the performance of an immiscible carbon 
dioxide WAG process are discussed. 


A - Linear Corefloods 
6.9.1 - Vertical WAG Injection 


Two vertical runs, VLC1 and VLC2, were performed by, respectively, 
injecting the carbon dioxide WAG at the bottom and top of the model to 
study the effect of gravity segregation of the injected fluids. The results of the 
two runs are given in Table 6.5 (page 83), along with those for the others. 
Figure 6.26 shows the producing GOR’s of the two runs, along with that for 
the horizontal WAG injection run’ (GTD6) conducted at identical conditions 
and utilizing identical experimental parameters. 


The figure shows that a smaller volume of carbon dioxide was pro- 
duced when carbon dioxide was injected at the top than at the bottom of the 
model. This is because when carbon dioxide is injected at the bottom, the 
buoyancy forces (acting in the same direction as the viscous forces) cause the 
carbon dioxide to rise to the top, resulting in early carbon dioxide break- 
through and the production of high carbon dioxide volume. The comparison 
of the GOR’s reveals that at the same carbon dioxide injection rate (i.e. 308.0 
cc/h) more gas was produced in the horizontal WAG injection run’ (GTD6) 
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Figure 6.26 - Producing Gas-Oil Ratios of Runs VLC1, VLC2, and GTD6’. 
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than in Runs VLC1 and VLC2 and that, in the horizontal WAG injection run, 
carbon dioxide production started after about 0.25 PV of fluids injected, show- 
ing that at the same injection rate the viscous forces acting in the horizontal 
direction (Run GTD6’) were greater than the sum of the buoyancy forces and 
viscous forces both acting the vertical direction (Runs VLC1 and VLC2). 


To examine the gravity segregation effect on water slugs, it is necessary 
to refer to Figure 6.27, where the producing WOR’s of the three runs are 
shown. In contrast to what was observed with the carbon dioxide slugs, 
injecting water slugs at the top of the model led to early water breakthrough 
and high WOR’s, as indicated in the figure. This was basically due to a higher 
water density compared to oil. Furthermore, when carbon dioxide dissolved 
in oil, in addition to reducing the oil viscosity it reduced the oil density, 
which induced a greater oil-water density difference, which consequently en- 
couraged early water breakthrough. 


The immiscible carbon dioxide WAG process for different model posit- 
ions: top (Run VLC2), bottom (Run VLC1), and horizontal (Run GTD6’) has 
an effect on the oil recovery as well. As is shown in Figure 6.28, the highest 
oil recovery was obtained when the WAG injection was conducted at the bot- 
tom of the model. The reason is that when carbon dioxide was injected at the 
bottom, gravity helped to induce the mass transfer rate of carbon dioxide by 
solution and diffusion into the oil, causing a significant reduction in the vis- 
cosity and density of the oil. When water slugs were injected following the 
injection of carbon dioxide slugs, due to its greater density compared to that of 
carbon dioxide-oil mixture, water resided at the bottom and pushed the car- 
bon dioxide-oil mixture upwards. Gravity also helped to stabilize the water 
front, resulting in a high displacement efficiency. As a result of these two ad- 
vantages, a high recovery was achieved in Run VLC1. 


In Run VLC2, where the immiscible carbon dioxide WAG was injected 
at the top rather than at the bottom, gravity kept the injected carbon dioxide at 
the top due to its lower density compared to that of oil, causing an adverse ef- 
fect on the mass transfer rate of carbon dioxide by solution and diffusion into 
oil. When water was injected, due to its higher density, it quickly found a 
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Figure 6.27 - Producing Water-Oil Ratios of Runs VLC1, VLC2, and OTD. | 
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Figure 6.28 - Comparison of Oil Recoveries of Runs VLC1, VLC2, and 
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way to finger through and pushed the oil upwards, resulting in the produc- 
tion of oil at high WOR’s. 


6.9.2 - Continuous Carbon Dioxide Injection 


Runs VLC3 and VLC5 were conducted to demonstrate more clearly 
segregation of the injected fluids. Continuous injection of carbon dioxide was 
carried out until 20% HCPV of carbon dioxide was injected, which was then 
followed by the injection of water to displace oil. 


Figure 6.29 shows the GOR curves of the two runs. Examining the 
GOR curve of Run VLC5, where carbon dioxide injection was conducted at 
the top of the model, reveals that gravity kept carbon dioxide at the top, in- 
ducing the formation of a gas cap which pushed the oil downwards. This 
can be observed by looking at the volume of oil produced after injecting 20% 
HCPV carbon dioxide. As indicated by the recovery curve of Run VLC5 in 
Figure 6.30, after 20% HCPV carbon dioxide was injected, approximately 2.0% 
of the oil was recovered, while none was produced when injecting the iden- 
tical volume of carbon dioxide at the bottom (Run VLC3). This volume of 
oil produced was relatively small compared to the volume of carbon dioxide 
injected. However, it clearly demonstrates that in the immiscible carbon 
dioxide WAG process, as time goes on, carbon dioxide will segregate and rise 
to the top to form a gas zone which pushes the reservoir oil downwards. 
The GOR curve of Run VLC3 shows a totally different trend. As shown in 
Figure 6.29, after the injection of 20% HCPV carbon dioxide and 0.95 total 
fluid PV, a very small amount of carbon dioxide was produced, showing that 
gravity segregation of the injected carbon dioxide does not occur right away. 
As more water was injected, the water pushed the carbon dioxide upward. 
Since carbon dioxide occurred as free gas, it could travel upward at a much 
higher rate than oil, resulting in the increased production of carbon dioxide. 
This clearly shows that water slugs injected following the carbon dioxide 
slugs in the WAG process induce the gravity rise of carbon dioxide to the top 
of the reservoir. 


The continuous water injection following carbon dioxide injection in 
Run VLC5 helps to illustrate the gravity segregation of the water slugs in the 
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Figure 6.29 - Comparison of the Producing GOR's of Runs VLC3 and VLC5. 
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Figure 6.30 - Comparison of the Oil Recoveries of Runs VLC3 and VLCS. 
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WAG process. It is indicated by the WOR curve of Run VLC5 (where water 
was injected at the top) in Figure 6.31 that high volumes of water were col- 
lected during the course of the experiment, which is thought to have a detri- 
mental effect on the immiscible WAG displacement process. To clarify this, 
the GOR’s, WOR’s, and recovery of Run VLC5 are plotted on the same graph, 
as shown in Figure 6.32. It is demonstrated by this figure that at the same 
pore volume of fluids injected much more water was produced than carbon 
dioxide and oil, which shows that in the immiscible WAG process the slugs 
of water push carbon dioxide and oil to the top of the reservoir, which causes 
the oil to re-saturate the zone already swept by the carbon dioxide gas. The 
same phenomenon can also be observed by considering Figure 6.33, where the 
GOR’s, WOR’s and recovery of Run VLC3 are plotted. 


Based on the phenomena observed in these two runs and those de- 
scribed in Section 6.9.1, it can be said that the immiscible carbon dioxide WAG 
displacement process performs the best when carbon dioxide and water are in- 
jected at the bottom. 


6.9.3 - Effect of Pressure on Gravity Rise of Carbon Dioxide 


The effect of pressure on the gravity rise of carbon dioxide was ob- 
served by Conducting Run VLC6 at 2.5 MPa. The experimental parameters 
used in this run, as included in Table 6.5 (page 83), were similar to those used 
in Run VLC1 in order to make a fair comparison. The tabulated experimental 
data of this run are provided in Appendix E. Figure 6.34 depicts the produc- 
tion history of this run. 


Figure 6.35 details the GOR’s of Runs VLC1 and VLC6 at each time flu- 
ids were injected. As already mentioned, Run VLC1 was carried out at 1.0 
MPa. An observation of the GOR curves of both runs reveals that pressure 
has the effect of retarding the gravity rise of the injected carbon dioxide gas. 
Gas breakthrough occurred at 0.48 PV in Run VLC1 (conducted at 1.0 MPa), 
whereas it occurred at 0.88 PV in Run VLC6 (conducted at 2.5 MPa). This is 
basically due to higher carbon dioxide solubility in Run VLC6. Other factors 
contributing to the earlier gas breakthrough in Run VLC1 can be as follows. 
When carbon dioxide dissolves in oil, it swells the oil and forms dispersed 
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Figure 6.31 - Comparison of the Producing WOR's of Runs VLC3 and 
VIEGS. 
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Figure 6.32 - Performance of Run VLC5 (Continuous Injection @ Top). 
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Figure 6.33 - Performance of Run VLC3 (Continuous Injection at Bottom). 
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NOTE: Average Run Conditions: Vertical WAG Flood at Bottom, 2.5 MPa, 21°C 
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Figure 6.34 - Production History of Run VLCO06. 
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Figure 6.35 - Effect of Pressure on Gravity Rise of Carbon Dioxide Gas. 
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bubbles of carbon dioxide gas, which travel faster than the oil, at a velocity 
equal to the velocity of the oil plus the rise velocity due to the buoyancy of the 
bubbles. It is known that the size of the bubbles formed depends on the pres- 
sure, assuming that spherical gas bubbles of a single and constant size are 
formed. The lower is the pressure, the bigger are the bubbles. The larger bub- 
bles will rise faster because the buoyancy forces grow faster than the drag 
forces. The buoyancy forces depend on the bubble volume (or V « r?), while 
the drag forces depend on the bubble surface area (A ar”). That is why gas 
breakthrough occurred earlier in Run VLC1 than in Run VLC6. Moreover, 
having been conducted at 1.0 MPa, the density of carbon dioxide injected in 
Run VLC1 was lower, which also contributed to a faster rate of gas rising. 


After 1.0 PV of fluids had been injected, more gas was produced in Run 
VLC6 than in Run VLC! as is shown in Figure 6.35. This is because a higher 
volume of carbon dioxide gas at standard conditions was injected in Run 
VLC6. 


In short, the gravity rise of the carbon dioxide gas injected in the im- 
miscible WAG process is affected by the operating pressure. Gravity rise oc- 
curs faster at a lower operating pressure. Therefore, it can be speculated that 
in the miscible carbon dioxide displacement process conducted at a pressure at 
which there is only one fluid phase present, the gravity rise of carbon dioxide 
is less severe than in the immiscible one. 


6.9.4 - Effect of Gas Injection Rate 


Two runs, VLC7 and VLC9, were conducted in the linear model by con- 
tinuously injecting carbon dioxide gas at the top of the model until gas break- 
through at the bottom. The injection rates used in the two runs were 0.984 
and 0.492 m/d, respectively. The pressure and temperature at which they 
were conducted were, respectively, 1.0 MPa and 21°C. After gas breakthrough, 
the total volume of carbon dioxide injected was recorded and water was then 
injected at the bottom, while the model was still in the vertical position. 


Figure 6.36 shows the GOR histories of the two runs. An interesting 
feature of this plot becomes evident by considering the volume of carbon 
dioxide injected until gas breakthrough. In Run VLC7, utilizing a 0.894 m/d 


eas 


cine 2 ta die als re 
“ul leo nud ey 4 3 (oa 
eva ek nO ahiaaah & Sarre) es pri Leal 
op ja TAsietSo. bee cate 
che) Meeyigl oT aaicind sith 7 
Li waegl apie ory 
Ly +00) rae 
yO ait? AM e 


an 
- 


p peanee es 

at AOL 108 AF a Eh 8 le ites batniaog be si3 
ee ee seeder dia se ) 
atilehy 249 save teat a od ne Ty nti 

_ 

ast oiled bisa nag oesite Hh sod bac abl OE at 

in sarc abet ee pegh at ower al im IY et D ie 

Lt if Dim tr salar Gashrwes4 is ary ohixnihy nites. uot ual 


rt 
ae Gt ; 


fal oA al be wget may biped ebayer ent io way dik) art se 

or ate Ori) SNe peteryes. ah yt eres 4 vexony DAW siding ; 

ifs vinings sei 4) nh? coun: aaa reno 28 oe 

5 ob pay 0 ie Copia wvaity iranieosyeih ebixoly moches slelioatin He 
we ee yp Se a a eu 

1) (ORIG mm writ} vie nee Aneegerity arg binfi sey giao we ost ms 

) we a acai ses i st 9 


: a 7 wien cohol alee: 


‘store ent ane ie ers ove ye tea 
Aneet 12g 7 sid ig! sabe ee hep its rtentri 


4 


sere). at aul “eit th 


oth weet vai 


bate 7 
» ee 7 hoki) — 1T. 
Gay deer ae aes | >. 
yA) ef 


, . F , 7 * 7 nin we i" 
ar saqy Yee nowt Pr Di SOOM AK Te Ory 
AOR AS ve 1G) AS Noeanw Tahoe sb Ser oe yh oh 


2 aa 
a esi elite aX a a: 
Bel? ie (No em 


at : eae rte 


ry) : 
Thee 


3 


25 


12 


-o- Run VLC07 (GOR) 
Oo Run VLCO9I (GOR) 
-e- Run VLCO07 (WOR) 


—#- Run VLCO9 (WOR) 


Legends 


10 


20 


water — 


Gas 
Breakthrough 


ee i ee ‘O1VeY [IO-A9}e&M 


it's) i=) 
rr rf ite) 


\o st N 


us / ue ‘oney [IO-seD 


4.0 a0 


3.0 


2.0 
Fluid Injected, Cum. PV 


1.0 


0.0 


Effect of Gravity on Injection Rate. 


Figure 6.36 - 


injection rate, gas breakthrough was noted when 0.8 PV of carbon dioxide was 
injected while it was noted when 1.0 PV of carbon dioxide was injected in 
Run VLC9 utilizing a 0.492 m/d injection rate. This indicates that at a lower 
rate, the residence time is longer, permitting more gravity segregation to oc- 
cur. The figure also shows that the GOR of Run VLC7 is approximately twice 
that of Run VLC9, at breakthrough. This indicates that in a downdip dis- 
placement, when carbon dioxide is injected at the top, a higher injection rate 
will result in a longer mixing zone between carbon dioxide and oil. 


Also presented in Figure 6.36 are the GOR’s of both runs when water 
was injected at the bottom. As shown, both GOR curves have nearly the 
same trend and there is little difference between the two even though the wa- 
ter injection rate employed in Run VLC7 was twice that used in Run VLC9. 
This shows that over the range of the water injection rate employed to dis- 
place the carbon dioxide-swollen oil in the immiscible process the gravita- 
tional forces dominate the segregation of the gas. 


The effect of gravity on the rise of the injected water can be viewed in 
Figure 6.36. As indicated in the figure, the WOR’s of Run VLC7 were higher 
than those of Run VLC9, showing that gravity had less effect at a water injec- 
tion rate of 0.984 m/d than at 0.492 m/d. In other words, had gravity the same 
effect at a water rate of 0.984 m/d as it had at 0.492 m/d, then the two WOR 
curves would overlap. Furthermore, it also showed that at 0.984 m/d the vis- 
cous forces played a more dominant role than the gravitational forces. 


The volume of oil displaced in each run is shown in Figure 6.37. At 
gas breakthrough, 5.1% of the oil was displaced by carbon dioxide in Run 
VLC7, whereas 9.5% of the oil was displaced in Run VLC9. Overall, as shown, 
more oil was displaced in Run VLC9 than in Run VLC7, because more carbon 
dioxide was injected until breakthrough in the former run. 


6.9.5 - Effect of Inverting the Core 


Run VLC10 was done a little differently from those already mentioned. 
In this experiment, carbon dioxide was injected at the top until breakthrough, 
then the core was inverted to injected water at the bottom. It was performed 
at 1.0 MPa and 21°C, utilizing an injection rate of 0.492 m/d. In other words, 
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Figure 6.37 - Effect of Gravity on the Volume of Oil Displaced. 
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the experimental conditions and parameters utilized in this experiment were 
identical to those utilized in Run VLC9, except that the model was inverted to 
inject water at the bottom. 


The GOR for this run is plotted in Figure 6.38. As shown, by inverting 
the model, the GOR at the first volume of water injected was 2.75 sm?/sm?, 
which is 2.5 times less than that of Run VLC9 which was done without in- 
verting the model. The reason for the lower GOR at the first volume of water 
injected is as follows. After carbon dioxide was injected until breakthrough, 
the carbon dioxide saturation at the top was higher than that at the bottom. 
When water was injected following inversion of the model, the gas (which 
was at the bottom) now at the top was produced. This indicates that it still 
took time for a larger volume of carbon dioxide at the bottom (which had 
been initially at the top) to rise to the top even though with the help of the 
driving forces provided by the water injected at the bottom. It also shows that 
the solubility of carbon dioxide in oil plays a role in retarding the rise of car- 
bon dioxide by virtue of gravity segregation. This is to say that had there been 
no carbon dioxide solubility in oil, all carbon dioxide would have risen to the 
top immediately once the model was inverted and water injected. 


B - Two-Dimensional Floods 
6.9.6 - Effect of Rate 


The immiscible carbon dioxide WAG process was proved to be success- 
ful in recovering oils from thin, deep and moderately heavy oil reservoirs 
underlain by a bottom water layer!®. A number of experiments conducted in 
a quarter of a five-spot system utilizing oils of different viscosities were made 
in this study to investigate the effect of injection rate. They were Runs H2D1 
to H2D25. All these experiments were conducted at 2.5 MPa, 21°C and utiliz- 
ing a carbon dioxide slug volume of 20% HCPV, a 4:1 WAG ratio and injec- 
tion rates varying from 0.78 to 3.81 m/d. The oils used in these experiments 
had viscosities from 603 to 3607 mPa.s. This viscosity range is normally en- 
countered in thin, deep, and moderately heavy oil reservoirs in Alberta and 
Saskatchewan. Table 6.6 (pages 84 to 86) contains a summary of these experi- 
ments. 
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Figure 6.38 - Effect of Inverting the Model. 
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Figure 6.39 presents the results of the experiments expressed as the oil 
DpH 
skAp,, 


recovery vs. the dimensionless group . It is shown in Table 6.6 (pages 85 


and 87) that for all oils tested the recovery climbs to its maximum as the injec- 
tion rate increases to 2.54 m/d and falls off when the injection rate exceeds 
2.54 m/d, indicating that 2.54 m/d is the optimal injection rate. The reasons 
could be as follows. When carbon dioxide is injected at a lower rate, it has a 
longer time to diffuse into the oil phase. The reverse is true when carbon 
dioxide is injected at a higher rate. Therefore, the higher is the volume of 
carbon dioxide that goes into solution in the oil, the greater is the oil viscosity 
reduction. The greater oil viscosity reduction always leads to a higher oil 
recovery. In the immiscible carbon dioxide WAG process, carbon dioxide and 
water are generally injected at the same rate. This means that when the 
injection rate of carbon dioxide is low, the injection rate of water is as well 
low. When water is injected at a low rate, the gravity effect is large. The vis- 
cous-to-gravitational force ratio is small. As a result, the water injected, in- 
stead of displacing oil, segregates at the bottom of the model or flows down- 
ward vertically to the bottom, which causes the displacing front to be nearly 
flat or horizontal. Therefore, only portions of oil near and at the bottom of 
the model are removed, resulting in a poor volumetric sweep and hence low 
oil recovery. On the other hand, when water is injected at a high rate, it will, 
instead of displacing the oil, bypass the oil, leading to a poor volumetric 
sweep and hence low oil recovery. This is why a lower oil recovery was ob- 
tained when an injection rate lower or higher than 2.54 m/d was employed. 


Figure 6.40 shows a similar plot for field conditions, assuming the labo- 


ratory recovery equals the field recovery; but the results were expressed as the 
Qyto ae. 
oil recovery vs. the dimensionless group a ARE A similar correlation was 
4kgAp,, 
done for runs conducted at 1.0 MPa. Figure 6.41 contains such a correlation. 
Similar to the observation made on Figure 6.39, the volume of oil produced 
increases as the injection rate increases as shown in Figure 6.41. The volume 
of oil produced is maximum at the injection rate of 2.54 m/d. Combining this 
observation and the one made on Figure 6.39 reveals that the optimal injec- 


tion rate for the two experimental pressures: 2.5 and 1.0 MPa is 2.54 m/d. 
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To investigate further the effect of carbon dioxide injection rate and to 
determine the optimal carbon dioxide injection rate at low pressures, a num- 
ber of experiments were performed by injecting carbon dioxide gas at different 
injection rates. They were Runs H2D26 to H2D29. The gas injection rates 
employed in the four runs were, respectively, 0.26, 0.52, 2.6, and 5.2 m/d, 
while the water injection rate was maintained at 2.6 m/d. The volume of 
carbon dioxide injected in these runs was 20% HCPV, and the WAG ratio was 
4:1. The instantaneous GOR’s of the four runs are plotted in Figure 6.42. The 
effect of gas injection rate becomes evident by observing the GOR curves in 
the figure. Gas breakthrough occurred earliest at the highest injection rate, 
i.e., at 5.2 m/d, and more gas was produced in Run H2D29 than in Runs 
H2D26, H2D27 and H2D28. Observing the GOR curves of Runs H2D26 and 
H2D27 raises a very interesting feature about the effect of gas injection rate on 
the solubility of carbon dioxide in the oil. There was very little or almost no 
production of gas in these two runs until 0.9 PV of fluids were injected, indi- 
cating that the low injection rates (0.26 and 0.52 m/d) used in these two runs 
helped retard gas breakthrough and gave the carbon dioxide enough time to 
dissolve as much as it could at the pressure and temperature conditions in 
the oil and thus to establish a close (but not complete) phase equilibrium with 
the oil. It should be noted that, based on experience, in order for carbon diox- 
ide to reach complete phase equilibrium with the oil, mechanical mixing of 
carbon dioxide and oil is needed for at least 3 weeks or longer; whereas, in 
Runs H2D26 and H2D27, the flood time was about 3 days. After the occur- 
rence of gas breakthrough, the GOR’s of both runs were almost similar, show- 
ing that injecting carbon dioxide at 0.26 or 0.52 m/d produced the same effect. 


Figure 6.43 compares the oil recoveries of the four runs. There is very 
little distinction in the oil recovery when carbon dioxide was injected at 0.26 
(Run H2D26), 0.52 (Run H2D27), or 2.6 m/d (Run H2D28). A lower recovery 
was noted when carbon dioxide was injected at 5.2 m/d (Run H2D29). 


To sum up, the injection rate of carbon dioxide has an important effect 
on the performance of the immiscible carbon dioxide WAG process. That is, 
when it is very low, it helps to delay gas breakthrough and reach more com- 
plete phase equilibrium between carbon dioxide and oil than when it is high. 
The experimental results show little or no effect of gravity on the gas 
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injection rate, particularly at low rates, because the experiments were con- 
ducted in short time spans. It is also obvious from the experimental results 
that, based on the oil recovery, the immiscible carbon dioxide WAG process 
performs best when carbon dioxide is injected at a rate less than or equal to 
the optimal water injection rate. This is true in the laboratory. In the field, 
because the reservoir thickness and well spacing are much larger than those 
in the model and the flood can take many years, injecting carbon dioxide at a 
rate lower than the optimal water injection rate will result in gravity segrega- 
tion of the carbon dioxide, which will cause gas tonguing, consequently caus- 
ing an adverse effect on the performance of the process. As such, for field ap- 
plications, it is recommended that carbon dioxide be injected at a rate equal to 
the optimal water injection rate. 


6.9.10 - Effect of Time 


The effect of time on the gravity segregation of carbon dioxide was 
investigated by conducting Runs H2D34 to 37, which were performed in a 
two-dimensional model utilizing a carbon dioxide gas slug volume of 5% 
HCPV at 2.5 MPa and 21°C. Other parameters used in these runs are provided 
in Table 6.6 (page 86). These experiments were done differently from those 
already discussed. After 5% HCPV of carbon dioxide was injected, a certain 
amount of time was allowed to let carbon dioxide “soak” into the oil, then 
water was injected to displace the oil. The soak times for the runs were, re- 
spectively, 0, 3, 4.83 and 10 days. Note that the same injection rate was em- 
ployed in all runs. 


Figure 6.44 shows the GOR’s of all four runs as a function of the vol- 
ume of fluids injected. A close observation of the GOR at breakthrough re- 
veals an interesting aspect of the effect of soak time on gravity segregation of 
the injected carbon dioxide gas: increasing the soak time increased the break- 
through GOR. The breakthrough GOR’s for soak times of 07377483, and 10 
days were respectively 0.06, 0.16, 0.27, and 0.33 sm3/sm%. A more obvious 
demonstration of the effect of soak time on gravity segregation of the injected 
carbon dioxide gas is shown by Figure 6.45, where the breakthrough GOR’s are 
plotted versus soak time. As noted, the breakthrough GOR increased by 3, 4.5, 


and 5 fold when the soak time increased to 3, 4.83, and 10 days, respectively. 
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Figure 6.45 - Effect of Soak Time on Gas Breakthrough. 
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The breakthrough GOR increases with soak time because as more time is al- 
lowed more carbon dioxide gas segregates and rises to the top by virtue of the 
density difference between carbon dioxide and oil, thus leading to a higher 
GOR at breakthrough. Based on this observation, it is speculated that gravity 
segregation of carbon dioxide will continue until no more segregation is pos- 
sible, assuming enough time is allowed. At that time, because the concentra- 
tion of carbon dioxide on the top is higher than that on the bottom, carbon 
dioxide will diffuse downward. Then, the GOR at breakthrough will be con- 
stant with soak time. 


The problem of gravity segregation becomes more pronounced with 
time. Based on this laboratory observation, it can be said that for field applica- 
tions where the flood can be as long as many years, the gravity segregation of 
the carbon dioxide gas will reduce the contact between carbon dioxide and oil, 
thus leading to a drastic reduction in the amount of carbon dioxide going into 
solution in the oil as the carbon dioxide moves farther away from the injec- 
tion well. 


6.9.11 - Gravity vs. Transverse Diffusion 


Two experiments, V2D1 and V2D2, were performed to investigate the 
role of transverse diffusion, normal to the vertical direction, on delaying the 
gravity rise of carbon dioxide. The gases used in Runs V2D1 and V2D2 were, 
respectively, carbon dioxide and nitrogen. The reason for the choice of nitro- 
gen in the second experiment was because it was to investigate only the grav- 
ity rise of the injected gas without any mass transfer involved. Nitrogen has 
a very low solubility and diffusivity in oil, compared to other gases. 


The two experiments were done as follows. With the two-dimensional 
model in the vertical position, 5% HCPV of carbon dioxide or nitrogen was 
injected at the bottom. Then, the model was left undisturbed until gas was 
detected at the top. The time when the gas was detected was recorded. 
Finally, water was injected at the bottom to complete the experiment. Table 
6.7 (page 87) summarizes the results of the two runs. 
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In Run V2D1, carbon dioxide was noted at the top after 65.4 days, while 
nitrogen was noted after 37.6 days in Run V2D2. The longer time in Run 
V2D1 indicates that the transverse diffusivity of carbon dioxide in the direc- 
tion normal to the longitudinal direction helped delay the gravity rise of the 
carbon dioxide gas travelling at a velocity equal to the sum of the diffusive 
velocity due to diffusion and the convective velocity due to gravitational 
forces. 


II.2 - Non-Isothermal Experiments 


In addition to conducting experiments under isothermal conditions, it 
was also intended to conduct experiments under non-isothermal conditions. 
As has already been defined, non-isothermal in this study means fluids at the 
atmospheric temperature are injected into a reservoir at a higher tempera- 
ture. Hence, phase changes and heat transfer from the reservoir fluids to the 
injected fluids will take place. 


6.10 - Unscaled Experiments 
6.10.1 - Effect of Temperature 


To investigate the effect of temperature, Run H2D38 was performed at 
37°C and 2.5 MPa in the two-dimensional model. Note that this experiment 
was not a scaled experiment. The 37°C temperature was arbitrarily selected. 
A carbon dioxide volume of 20% HCPV divided into ten equal slugs and a 4:1 
WAG ratio were utilized in this run. The oil used had a viscosity of 1058 
mPa.s at 21°C and atmospheric pressure. Both carbon dioxide and water were 
injected at 2.54 m/d. In short, the experimental parameters utilized in this 
run were identical to those used in Run H2D8 conducted at 21°C. The only 
difference between the two was the experimental temperature. The tabulated 
experimental data for this run can be found in Appendix E. 


Figure 6.46 details the producing GOR’s and production history of Run 
H2D38. The compositions by mole percent of the gas produced in this run 
were measured to be on average 0.39 propane, 6.96 water, 1.41 methane, and 
91.24 carbon dioxide, while the produced gas in Run H2D8 was found to be 
mostly carbon dioxide and negligible amounts of water and methane. This 
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Figure 6.46 - Production History of Run H2D38. 
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shows that the temperature affected the composition of the gas produced in 
Run H2D38. It caused some of the injected and reservoir water to vaporize 
and made the light components of the oil more volatile. 


Figure 6.47 presents a comparison of the GOR’s of Run H2D38 with 
those of Run H2D8. As shown, gas breakthrough occurred earlier in Run 
H2D8 (21°C) than in Run H2D38 (37°C). This is due to the lower gas-oil mo- 
bility ratio encountered in the latter run. The viscosity of the oil in Run 
H2D8 was higher than that in Run H2D38 because the latter was conducted at 
37°C while the former was conducted at 21°C. At 37°C and atmospheric pres- 
sure, an oil with a viscosity of 1058 mPa.s at 21°C and atmospheric pressure 
was measured to be 327 mPa.s. Appendix G contains the experimentally mea- 
sured viscosities of different oils as functions of temperature. Empirical cor- 
relations based on Walther’s equation were also made and are included. 
Even though gas breakthrough occurred earlier in Run H2D8, the GOR curve 
for Run H2D38 is higher than that for Run H2D8, as shown in the figure. 
There were two factors which caused the higher gas production in Run 
H2D38. The first was the higher temperature which reduced the solubility of 
carbon dioxide in oil. The second might be as follows. Due to the presence of 
water vapour, propane, and methane in Run H2D38 (conducted at 37°C), 
when carbon dioxide was injected, it mixed faster with the water vapour, 
methane, and propane than with the oil because mass transfer occurs much 
faster from gas to gas than from gas to liquid. Once carbon dioxide mixed 
with these gases, its solubility in oil decreased. However, the diffusivity of 
carbon dioxide was higher at 37°C, but this did not contribute much because 
diffusion is a very slow mass transfer process compared to solution. 


Figure 6.48 compares the oil recovery of Run H2D08 with that of Run 
H2D38. As shown, the recoveries were 54% and 45% for Runs H2D38 and 
H2D8, respectively. This clearly indicates that temperature affects the dis- 
placement efficiency of the immiscible WAG process. 


6.10.2 - Effect of Oil Viscosity 


It was noted in the preceding section that more oil was recovered at a 
higher temperature because the viscosity of a 1058 mPa.s oil at 21°C reduced to 
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Legends 
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Figure 6.47 - Comparison of GOR's of Runs H2D8 and H2D38. 
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Figure 6.48 - Comparison of Oil Recoveries of Runs H2D8 and H2D38. 
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327 mPa.s at 37°C. Consequently, it was intended to conduct an experiment 
using an oil with a viscosity equal to 1058 mPa.s at 37°C. Therefore, an oil 
with a viscosity of 5200 mPa.s at 21°C was chosen because it viscosity became 
1058 mPa.s at 37°C. Run H2D39 was conducted using this oil. The experi- 
mental parameters employed in this run were exactly the same as those used 
in Runs H2D8 and H2D38. Note that it was done at 37°C. Appendix E con- 
tains the tabulated experimental data for this run. 


The instantaneous GOR’s of Run H2D39, together with those of Run 
H2D8, are shown in Figure 6.49. It is shown in the figure that gas break- 
through occurred earlier in Run H2D39 and that the GOR curve for Run 
H2D39 is higher than that for Run H2D8, even though the two runs were 
conducted using oils with the same viscosity at two different temperatures. 
The earlier gas breakthrough and higher instantaneous GOR’s occurring in 
Run H2D39 were due to the lower gas viscosity and solubility of carbon diox- 
ide in oil at a higher temperature, which did not cause the same viscosity re- 
duction as it did in Run H2D8. 


The recoveries of the two runs are shown in Figure 6.50. The curve 
with a lower trend represents the recovery history of Run H2D39. It is clearly 
shown in the figure that the total recovery of Runs H2D39 is 4% lower than 
that of H2D8, which was basically, as has been recently mentioned, due to the 
lower carbon dioxide solubility at a higher temperature. 


The observation made above reveals that in order to correlate a non- 
isothermal experiment to an isothermal one or vice-versa, not only the oil 
viscosity must be considered but also the solubility of carbon dioxide in the 
oil. 


6.10.3 - Effect of Carbon Dioxide Solubility 


Based on the observation made in the preceding section, it was in- 
tended to conduct an experiment at 37°C and at a predetermined pressure 
such that the solubility of carbon dioxide in oil at this pressure and tempera- 
ture condition was identical to that at 21°C and 2.5 MPa, ignoring the effects of 
water vapour and hydrocarbon gases at the higher temperature. To deter- 
mine this experimental pressure, Chung, Jones, and Nguyen’s correlation! 
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Figure 6.49 - Comparison of GOR's of Runs H2D8 and H2D39. 
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Figure 6.50 - Oil Recovery Comparison of Runs H2D8 and H2D39. 
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for carbon dioxide solubility in oils was used. At 3.14 MPa and 37°C the solu- 
bility of carbon dioxide was identical to that at 2.5 MPa and 21°C. Hence, Run 
H2D40 was performed at 3.14 MPa and 37°C utilizing the same experimental 
parameters as used in Run H2D8. The same oil used in Run H2D39 was used 
in this run. Appendix E contains the tabulated results of this run. 


The instantaneous GOR’s of Run H2D40 are plotted in Figure 6.51, 
along with those of Run H2D8. Observing the GOR’s of both runs raises an 
interesting point. As indicated, gas breakthrough occurred at the same time 
in both runs. Also shown in the figure is that the GOR curve for Run H2D40 
has a higher trend than that for Run H2D8. This is because a higher carbon 
dioxide volume at standard conditions was utilized in Run H2D40. Another 
reason which is believed to cause the higher GOR’s in Run H2D40 is the 
lower solubility of carbon dioxide due to the higher concentration of water 
vapour at 37°C. 


Figure 6.52 compares the volume of oil displaced in the two runs. The 
two recovery curves nearly overlay, demonstrating that almost the same oil 
recoveries were obtained in the WAG, post-waterflood, and blowdown phases 
in both runs, and consequently the total recoveries of the two differ by less 
than 2% (45.4% for Run H2D8 vs. 43.7% for Run H2D40). This recovery dif- 
ference is small and can be neglected. The almost similar oil recoveries re- 
ported in the two runs indicate that the same viscosity reduction was approx- 
imately achieved in both runs. 


In short, combining the observations made in the preceding section 
and this, it can be concluded that in order to correlate an isothermal experi- 
ment to non-isothermal conditions, or vice-versa, the oil viscosities and 
carbon dioxide solubilities in the two runs be equal. 


6.10.4 - Effect of Slug Size 


Based on the observation made in Section 6.10.1, it is hypothesized that 
because the solubility of carbon dioxide is lower at a higher temperature, a 
smaller volume of carbon dioxide should be used instead. Run H2D41 was 
conducted using a total carbon dioxide slug size of 10% HCPV. In order to 
make a fair comparison with Run H2D339, the same oil with identical physical 
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Figure 6.51 - Effect of Similar Carbon Dioxide Solubility at Two 
Different Pressure and Temperature Conditions. 
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Figure 6.52 - Comparison of Oil Recoveries of Runs H2D8 and H2D40. 
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properties and the same injection rate and experimental parameters used in 
Run H2D39 were used in Run H2D41. The tabulated experimental results of 
this experiment can be seen in Appendix E. 


Figures 6.53 and 6.54 show the comparisons of Runs H2D41 and H2D339. 
As shown in Figure 6.53, similar to what was observed in Section 6.7.3, a 
higher producing GOR curve corresponds to a higher volume of carbon diox- 
ide injected. In Figure 6.54, a comparison of the recoveries of the two runs 
reveals that a higher oil recovery was obtained with a higher volume of car- 
bon dioxide injected. The total volume of oil recovered in Run H2D41 (10% 
HCPV CO?) was 31.2% while it was 41.1% in Run H2D39 (20% HCPV CO2). As 
explained previously, a lower recovery was obtained when a 10% HCPV of 
carbon dioxide was injected because the volume of carbon dioxide utilized 
was insufficient to cause the maximum oil viscosity reduction. 


6.11 - Scaled Experiments 
6.11.1 - Experimental Design 


The purpose of this section is to use the scaling criteria derived in 
Chapter 4 (page 18) to design the scaled experiments according to the field data 
presented in Table 5.1 (page 52). The first parameter to be determined is the 
scaling factor, a, which is the ratio of the length in the prototype to that in the 
model. It is usually determined by the maximum physical size available for 
the experiment. In this study, it is desired to design the scaled experiments in 
such a way that they can be conducted in the existing scaled physical model. 
The length of the existing model, as mentioned in Chapter 5 (page 47), is 0.457 
m; the scaling factor, a, is therefore 


_ field well spacing — 201.17 m 


a= = = 440.2. 
model well spacing 0.457 m 


Model Thickness: it is calculated by the group (#) as follows: 
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Figure 6.53 - Effect of Slug Size on GOR's at Non-Isothermal Conditions. 
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Figure 6.54 - Effect of Slug Size on Oil Recovery at Non-Isothermal 
Conditions. 
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Thus, the desired model representing the reservoir will have dimen- 
sions of 0.457 mx0.457 mx0.012 m. As it is very costly to fabricate a model 
with these dimensions, the existing scaled physical model with dimensions of 
0.457 m x 0.457 m x 0.022 m is used instead. The existing model is 0.010 m 
thicker that the desired model. 


Fluid Properties: each dimensionless fluid property should be the same func- 
tion of dimensionless pressure and temperature for the model as it is for the 
field. To satisfy this requirement, it is best to use the reservoir fluids in the 
model; thus, fluid properties such as density, viscosity, solubility, and diffu- 
sivity are automatically scaled. This is true if the model pressure and temper- 
ature are identical to the reservoir pressure and temperature. 


For Approach 1, the fluid properties are scaled because the reservoir 
fluids are used in the model; moreover, the pressure and temperature are the 
same in the model and the prototype. As for Approach 2, only two fluid 
properties are scaled: viscosity and diffusivity, because the viscosity of the 
model oil and the solubility of carbon dioxide in the model oil at the labora- 
tory temperature and experimental pressure are chosen to be respectively the 
same as those in the reservoir oil at the reservoir temperature and pressure. 
The viscosity of the reservoir oil at 37°C is 160 mPa.s; therefore, the viscosity 
of the model oil must be 160 mPa.s at 219°C. 


Model Pressure: for Approach 1, there is no need to determine the model 
pressure because it is chosen to be equal to the field pressure. As for App- 
roach 2, it is necessary to determine a model pressure such that the solubility 
of carbon dioxide at 21°C and the model pressure equals that at field tempera- 
ture and pressure, i.e., at 37°C and 4.8 MPa, respectively. By using Chung et 
al.’s correlation!, the model pressure is 3.58 MPa. 


In addition , it is necessary to select the pressure drop in the model 
such that the gravitational-to-viscous forces ratio is scaled. This can be done 


as follows: 


(p = Pprod eee ~ (AP o8H) ae 
(p FE: cee) (Ap og 8H) Field 
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Re-arranging gives 


(AP oF) ae 


(P—Pprod 4) =(P— Pprod Ju * 
Prod ] Model prod /Field (APogH)..., 


Therefore, the model pressure drop, by virtue of this expression, is 
reduced by a factor “a”. Given the densities of oils and carbon dioxide at the 
two pressure and temperature conditions, the true model pressure can be 
determined as shown below. 


Pco, = 83.5 kg/m? at 3.58 MPa and 21°C 
Pco, = 111.0 kg/m? at 4.8 MPa and 37°C 
Poit = 950.0 kg/m? at 3.58 MPa and 21°C 
Poi = 943.8 kg/m? at 4.8 MPa and 37°C 


and the field production pressure is 2.6 MPa. The expression for the model 
pressure can be written as 


PModel = (Pprod agg + 9-00987 


If the model production pressure is chosen to be 3.58 MPa, then the ac- 
tual average pressure of the model will be 3.59 MPa, which is not very differ- 
ent from the pressure determined using Chung et al’s correlation! for pre- 
dicting solubility of carbon dioxide in viscous crudes. 


Permeability Determination: for both approaches, the permeability is the 
same in the model and the prototype because the reservoir sand is used to 
represent the field porous medium. 


Injection Rate: to calculate the injection rate of carbon dioxide and water, the 
mal see : 

group ORM is used. The field current injection rate, being unavailable, 
PerKgrPprH 

the optimal injection rate determined in Section 6.9.6 (page 122) is used in the 


scaled experiments. 


nee i 


ei a | ; 
in nbieoth adig tee alo ts “ident orl} a; 
7L/ANgIh a oh at Ba reer 
fal 


_ 7 


- 


7 


7 
te: ae 
PREbas «1M 8.8 Set gal QT ona 2 
PTE Ste oT RE Ie Soh gpk G.Oee 
POS hee eM Reis ee 


Ye a teh mabe oe ee St ad sacar mobtsubonty hl 
bi PETIT SKE at 


FHOLO+ iced or a 


o) wet SFY el ed eeei> ob eapenstg adibsctona etna ey 
ban hs doitw 2 Tit ORS od May abo oft §> svsaeng o: ans 
vit Tomiie ole i ye elie ae 
aalsiig Sineodht, €k siiiales rake re 


ite oi f' i 2 _acbateart red tot capitadiamtetogh 
ai als >! Lani lowreten enh aipeoad exec a ber 
panty 


a r | 
& ssi Bite ell ae iiss il ec 
V sltawei griie’ sige ote 


af} 
= 
~ 
? 
td 


6.11.2 - Discussion of The Results of The Scaled Experiments 


Two scaled experiments were performed according to the design dis- 
cussed in the preceding section to investigate the possible application of the 
immiscible carbon dioxide WAG recovery technique in the pre-waterflooded 
reservoir under consideration (moderately heavy oil reservoir in Saskatche- 
wan) for which reservoir description is provided in Table 6.8 (page 147). They 
were tertiary Runs H2D42b and H2D43b. Run H2D42b was conducted using 
the design for Approach 1 and Run H2D43b was conducted using Approach 2. 
The total volume of carbon dioxide injected in both runs was 20% HCPV ina 
2:1 WAG mode. The 2:1 WAG ratio was selected because a previous study on 
this reservoir showed that utilizing a 2:1 WAG ratio is as effective as utilizing 
a 4:1 WAG ratio, for the reservoir under study. The temperature and pres- 
sure conditions at which Runs H2D42b and H2D43b were conducted were re- 
spectively 37°C at 4.8 MPa and 21°C at 3.58 MPa. Even though they were con- 
ducted at two different temperatures, the viscosities of the oils used in the 
two runs were similar, and the solubilities of carbon dioxide in the oils were 
supposed to be the same as well. Hence, a fair comparison could be made, as- 
suming there was no oil density difference. 


Figure 6.55 presents a comparison of the instantaneous GOR’s of the 
two runs. It is shown in the figure that the GOR curve of Run H2D42b (37°C) 
is higher for the most part than that of Run H2D43b (21°C). This is because a 
higher volume of water vapour was present in Run H2D42b. The water 
vapour from the analysis of the produced gas in Run H2D42b was 3.88 
mole%, while it was 0.81 mole% in Run H2D43b. This higher concentration 
of water vapour prevented carbon dioxide from going into the oil. That is, it 
reduced the carbon dioxide solubility, causing carbon dioxide to remain in a 
free gas phase which would then bypass the oil. Another interesting feature 
can be spotted when comparing the late time GOR’s of the two runs. The late 
time GOR of Run H2D42b (non-isothermal at 37°C), as shown in the figure, is 
3 times higher than that of Run H2D43b, which indicates that in the blow- 
down stage when the pressure is depleted, the higher temperature enhances 
the release of the carbon dioxide gas dissolved in the oil. 


The comparison of the recoveries of the two runs, according to the two 
recovery curves in Figure 6.56, reveals that the same volume of oil was 
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Figure 6.55 - GOR's of Run H2D42b vs. Those of Run H2D43b. 
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Figure 6.56 - Recovery of Run H2D42b vs. That of Run H2D43b. 
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recovered in the WAG stage (from 0 to 0.3 PV) and the post-water recoveries 
of the two runs differed by less than 0.5% HCPV. Overall, the total recoveries 
differed by less than 1.5% HCPV. This total recovery difference is very small, 
indicating that both approaches can be used to predict the performance of the 
immiscible carbon dioxide process under non-isothermal conditions. 


6.12 - Reproducibility of the Experimental Results 


An important aspect of this study concerns the reproducibility of the 
experimental results. 


Repeatability of the experiments was examined further by conducting 
Runs CWF3, CWF4, VLC4, VLC8 and H2D44b to check respectively the results 
of Runs CWF1, CWF2, VLC3, VLC7 and H2D43b. Run VLC11 was performed 
to check also Run VLC7. The results of these experiments were also tabulated 
and are contained in Appendix E. 


Figures 6.57, 6.59, 6.61, 6.63, and 6.65 present the reproducibilities of the 
GOR’s of Runs CWF1, CWF2, VLC3, VLC7, and H2D43b, respectively. They 
show that the reproducibilities of the GOR’s of these experiments were fairly 
good, because the overall errors involved were small-in the range of 15 to 30 
%. This range of error is normally encountered in any experiment involving 
the flow of a gas phase, because gases do not flow at a constant rate. Conse- 
quently, it is very difficult to have a constant gas flow rate. 


Figures 6.58, 6.60, 6.62, 6.64 and 6.66 present the reproducibilities of the 
recoveries of Runs CWF1, CW2, VLC3, VLC7 and H2D43b, respectively. 
Unlike those of the GOR’s, the reproducibilities of the oil recoveries were 
good, the errors being in the range of 5 to 10%. 
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Figure 6.57 - GOR Repeatability of Run CWF1. 
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Figure 6.58 - Recovery Repeatability of Run CWF1. 
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Figure 6.59 - GOR Repeatability of Run CWF2. 
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Figure 6.60 - Recovery Repeatability of Run CWF2. 
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Figure 6.61 - GOR Repeatability of Run VLC3. 
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Figure 6.62 - Recovery Repeatability of Run VLC3. 
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Figure 6.65 - GOR Repeatability of Run H2D43b. 
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7 - SUMMARY and CONCLUSIONS 


This investigation was devoted to several aspects of the immiscible 
(subcritical) carbon dioxide WAG process for the recovery of moderately vis- 
cous oils. In this study, experimental data reported by Rojas!, Zhu‘, Dyer”, 
Wilson’, and in this work were correlated with the dimensionless similarity 
groups derived in this work. A method supported by a mathematical model 
was developed to measure experimentally the diffusivity of a gas into a liq- 
uid. Seventy-six experiments were conducted to measure the diffusivities of 
carbon dioxide and methane in various oils. An empirical correlation for 
predicting the diffusivity of carbon dioxide in oil was developed using the 
data obtained. Vertical and horizontal displacement experiments were con- 
ducted in two scaled models to examine gravity segregation of the injected 
fluids. Horizontal displacement data were correlated with the dimensionless 
group describing the ratio of viscous-to-gravitational forces. In addition, a 
non-isothermal and non-equilibrium mathematical model including phase 
change and interfacial mass transfer was developed. Two sets of similarity 
groups for the non-isothermal immiscible carbon dioxide WAG displacement 
process were derived. Non-isothermal displacement experiments were also 
designed and performed. 


Based upon the experimental observations, the following conclusions 
can be reached. 


r ioxide Diffusion into Oi 


1. Diffusivities of carbon dioxide and methane increase with increasing 
pressure, as long as both gases are in the gas phase (2.3648E-09 at 0.69 MPa vs. 
6.1387E-09 at 6.89 MPa). 


2. Increasing temperature and/or decreasing the oil viscosity increases 
the diffusivity of carbon dioxide into oil (6.1387E-09 at 57.1°C vs. 7.5272E-09 at 
65.4°C). 


3. The diffusivity of carbon dioxide is affected by the molecular weight 
of the oil. It decreases as the oil molecular weight increases (2.5452E-08 in a 
516.73 g/mol oil vs. 4.9337E-10 in a 737.59 g/mol oil). 
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4. Carbon dioxide, having a smaller molecular size than methane, can 
diffuse faster in oils than methane (6.1387E-09 for carbon dioxide vs. 8.4217E- 
11 for methane). 


rrelation of Experimental D 


1. Dimensionless groups were derived for a comprehensive mathemat- 
ical model of the immiscible carbon dioxide flooding process. 


2. The groups can be used to correlate the previous and present exper- 
imental data. 


ravi r ion of rbon Dioxi 


1. Based on the results of the vertical displacement experiments con- 
ducted, it can be concluded that gravity plays an important role in the immis- 
cible carbon dioxide WAG displacement process. It causes the injected carbon 
dioxide to rise to the top, which results in the formation of a gas zone which 
is believed to both finger through and push the oil down towards the bottom, 
which would then be displaced by the injected water. Gravity segregation of 
the injected water slugs encourages the rise of carbon dioxide to the top and 
pushes the oil upwards, which will then re-saturate the zone already swept by 
carbon dioxide. This volume of oil will never be recovered, consequently 
causing the loss of oil recovery. 


2. In horizontal floods, transverse diffusion of carbon dioxide in the 
horizontal direction normal to the horizontal longitudinal direction can help 
delay the gravity rise of the gas. 


Non-lsothermal Experiments 


1. Experiments at an elevated temperature showed distinct effects on 
the mechanism of the process, mainly due to evaporation of water and mix- 
ing of water vapour with carbon dioxide, thereby reducing its solubility in oil. 


2. Experimental results show that if two oils have the same viscosity 
and gas solubility at two different temperatures, the overall performance of 
the immiscible carbon dioxide WAG process will be very similar in the two 
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cases, thereby indicating that the experimental results can be extended to 


other temperatures. 
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8 - RECOMMENDATIONS for FUTURE RESEARCH 


The following studies are recommended to extend the scope of this re- 
search. 


1. The effect of a bottom water layer on the performance of the immis- 
cible carbon dioxide process should be examined. This can be done by creating 
a water layer underneath the oil zone. This study should be carried out be- 
cause the field reservoir is underlain by a water layer and the experimental 
results will more accurately predict the true field performance of the process. 


2. The improvement in the relative permeability to oil afforded by car- 
bon dioxide at immiscible conditions should be studied. 


3. A study should also be conducted to investigate the effect of gas 
trapping which occurs due to gravity segregation. 
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APPENDIX A 


INITIAL AND BOUNDARY CONDITIONS 


A.1 - Initial and Boundary Conditions 
No fluid flow across top and bottom of reservoir 
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APPENDIX B 
DIFFERENTIAL EQUATIONS IN DIMENSIONLESS FORM 


B.1 - Momentum Balance for COz in the Oil Phase in Dimensionless Form 
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B.2 - Total Energy Balance for Carbon Dioxide in the Oil Phase in 
Dimensionless Form 
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APPENDIX C 


DERIVATION OF THE RELAXED SCALING GROUPS 


For the two approaches, the method of deriving a relaxed set of similar- 
ity groups is similar. Based on the assumption made in each approach, the 
terms corresponding to these assumptions are deleted from the governing 
partial differential equations. Each equation is then divided by one of its own 
remaining coefficients to yield the dimensionless form of the equation. The 
coefficients represent the relaxed set of similarity groups which can subse- 
quently be reduced to their simplest form. The constitutive relationships, 
constraints, and initial and boundary conditions are treated in a similar man- 
ner. 


C.1 - Approach No. 1 


For approach no. 1, the effects of gravity and transverse and longitudi- 
nal dispersion are assumed to be neglible. For the sake of simplicity the fol- 
lowing model for diffusion in porous media was adopted??. 
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C.2 - Approach No.2 


For approach no. 2, with the exclusion of the pressure drop term from 
the governing equation the mass balance for carbon dioxide in the oil phase 
becomes 
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APPENDIX D 


Tabulated Data of Diffusivity Experiments in Graphical Form 
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Figure D1 - Mass of CO} Injected vs. Square Root of Time for Diffusivity 
Experiments No. 1, 2, 3, and 4. 


1500 


1500 


203 


Eg 


m,, Mass of Gas Injected, g 


m;, Mass of Gas Injected, g 


0.12 0.30 


Experiment No. 6 


0.10 Experiment No. 5 ' 0.25 Ces r0sa5 ke/m° 
Cie = 0.4236 kg/m ; 
5D 
os 
0.08 2 0.20 
a 
k= 
1?) 
0.06 (5 0.15 
al 
© 
i?) 
S 
0.04 S 0.10 
= 
0.02 0.05 
0.009 0.00-4 
0 500 1000 1500 0 500 1000 
vt, s” vt, s”2 


15 


Experiment No. 7 
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Figure D2 - Mass of CO, Injected vs. Square Root of Time for Diffusivity 
Experiments No. 5, 6, 7, and 8. 
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Figure D3 - Mass of CO) Injected vs. Square Root of Time for Diffusivity 
Experiments No. 9, 10, 11, and 12. 
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Figure D4 - Mass of CO, Injected vs. Square Root of Time for Diffusivity 
Experiments No. 13, 14, 15, and 16. 
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Figure D5 - Mass of CO, Injected vs. Square Root of Time for Diffusivity 
Experiments No. 17, 18, 19 and 20. 
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Figure D6 - Mass of CO, Injected vs. Square Root of Time for Diffusivity 
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Figure D7 - Mass of CO» Injected vs. Square Root of Time for Diffusivity 
Experiments No. 25, 26, 27, and 28. 
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Figure D8 - Mass of CO; Injected vs. Square Root of Time for Diffusivity 
Experiments No. 29, 30, 31, and 32. 
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Figure D9 - Mass of CO, Injected vs. Square Root of Time for Diffusivity 
Experiments No. 33, 34, 35, and 36. 
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Figure D10 - Mass of CO} Injected vs. Square Root of Time for Diffusivity 
Experiments No. 37, 38, 39, and 40. 
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Figure D11 - Mass of CO, Injected vs. Square Root of Time for Diffusivity 
Experiments No. 41, 42, 43, and 44. 
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Figure D12 - Mass of CO} Injected vs. Square Root of Time for Diffusivity 
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Figure D13 - Mass of CO; Injected vs. Square Root of Time for Diffusivity 
Experiments No. 49,50, 51, and 52. 
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Figure D14 - Mass of CO} Injected vs. Square Root of Time for Diffusivity 
Experiments No. 53, 54, 55, and 56. 
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Figure D15 - Mass of CO? Injected vs. Square Root of Time for Diffusivity 
Experiments No. 57, 58, 59, and 60. 
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Figure D16 - Mass of CO, Injected vs. Square Root of Time for Diffusivity 
Experiments No. 61, 62, 63, and 64. 
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Figure D17 - Mass of CO, Injected vs. Square Root of Time for Diffusivity 
Experiments No. 65, 66, 67, and 68. 
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Figure D18 - Mass of CO2 Injected vs. Square Root of Time for Diffusivity 
Experiments No. 69, 70, 71, and 72. 
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Figure D19 - Mass of CH4 Injected vs. Square Root of Time for Diffusivity 
Experiments No. 73, 74, 75, and 76. 
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APPENDIX F 


Production Histories of All Experiments Conducted 
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NOTE: Average Run Conditions: Carbonated Waterflood, 2.5 MPa, 21°C 


Model Parameters: Average Injection Rate = 308 cc/hr, [, 
9 = 36.0%, k = 11.3 darcies, S, = 96.2 %, Swe = 3.8 % 
Model Type: Linear 
[1.608 moles of CO.] 

Figure F1 - Production History of Run CWF3. 
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NOTE: Average Run Conditions: Carbonated Waterflood, 1.0 MPa, 21° 
Model Parameters: Average Injection Rate = 308 cc/hr, L, = 1058.0 mPa.s, 


g = 37.5 %,k = 10.9 darcies, S, = 90.1 %, Sy. = 9.9 % 


Model Type: Linear 
[20% HCPV of CO, (0.143 moles) Mixed with Water at 4:1 Ratio] 


Figure F2 - Production History of Run CWF4. 
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NOTE: Average Run Conditions: Vertical WAG Flood at Bottom, 1.0 MPa, 21°C 
Model Parameters: Average Injection Rate = 308 cc/hr, 1, = 1053.0 mPa.s, 


9 = 35.9 %, k = 11.33 darcies, S, = 95.4 %, Swe = 4.6 % 


Model Type: Linear 


[0.20 HCPV CO, @ 1.0 MPa (0.089 moles) 4:1 WAG, 10 Slugs] 


Figure F3 - Production History of Run VLC1. 
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NOTE: Average Run Conditions: Vertical WAG Flood at Top, 1.0 MPa, 21°C 
Model Parameters: Average Injection Rate = 308 cc/hr, 1, = 1058.0 mPa.s, 


9 = 35.5 %, k = 11.12 darcies, S, = 90.7 %, Swe = 9.3 % 


Model Type: Linear 
[0.20 HCPV CO, @ 1.0 MPa (0.087 moles) 4:1 WAG, 10 Slugs] 
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Figure F4 - Production History of Run VLC2. 
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NOTE: Average Run Conditions: Continuous Injection at Bottom, 1.0 MPa, 2150 
Model Parameters: Average Injection Rate = 308 cc/hr, 1, = 1058.0 mPa.s, 


g = 34.74%, k = 9.18 darcies, Sy = 94.35%, Swe = 5.65% 


Model Type: Linear 
[0.20 HCPV CO, @ 1.0 MPa (0.088 moles)] 


Figure F5 - Production History of Run VLC3. 
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NOTE: Average Run Conditions: Continuous Injection at Bottom, 1.0 MPa, 21°C 
Model Parameters: Average Injection Rate = 308 cc/hr, lL, = 1058.0 mPa.s, 


g = 35.50%, k = 10.23 darcies, S, = 91.96%, Sy, = 8.04% 


Model Type: Linear 
[0.20 HCPV CO, @ 1.0 MPa (0.088 moles)] 


Figure F6 - Production History of Run VLC4. 
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NOTE: Average Run Conditions: Continuous Injection at Top, 1.0 MPa, 21°C 
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Model Parameters: Average Injection Rate = 308 cc/hr, [L, = 1058.0 mPa.s, 


§ = 35.38%, k = 11.34 darcies, S, = 91.82%, Sy, = 8.18% 


Model Type: Linear 
[0.20 HCPV CO, @ 1.0 MPa (0.088 moles)] 


Figure F7 - Production History of Run VES: 
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NOTE: Average Run Conditions: Continuous CO), Injection at Top, 1.0 MPa, Dine 


Model Parameters: Average Injection Rate = 308 cc/hr, , = 1055 mPa.s, 
@ = 35.38%, k = 10.02 darcies, S, = 95.27%, Sy, = 4.73% 


Model Type: Linear 
[0.904 HCPV CO, @ 1.0 MPa (0.391 moles), Water Injected @ Bottom] 


Figure F8 - Production History of Run VLC7. 
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NOTE: Average Run Conditions: Continuous CO) Injection at Top, 1.0 MPa, 213G 
Model Parameters: Average Injection Rate = 308 cc/hr, 1, = 1055 mPa.s, 


g = 35.45%, k = 9.44 darcies, S, = 95.55%, Swe = 4.45% 
Model Type: Linear 


[1.703 HCPV CO, @ 1.0 MPa (0.748 moles), Water Injected @ Bottom] 


Figure F9 - Production History of Run VLC8. 
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NOTE: Average Run Conditions: Continuous CO) Injection at Top, 1.0 MPa, 21°C 
Model Parameters: Average Injection Rate = 154 cc/hr, U1, = 1053.0 mPa.s, 


g = 35.38%, k = 10.58 darcies, S, = 96.18%, Sy, = 3.82% 
Model Type: Linear 


[1.015 HCPV CO, @ 1.0 MPa (0.484 moles), Water Injected @ Bottom] 


Figure F10 - Production History of Run VLC9. 
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NOTE: Average Run Conditions: Continuous CO, Injection at Top, 1.0 MPa, 21°C 
Model Parameters: Average Injection Rate = 154 cc/hr, 1, = 10550 mPa.s, 


g = 35.38%, k = 11.02 darcies, S, = 5.09%, Sy, = 4.91% 


Model Type: Linear 
[1.341 HCPV CO, @ 1.0 MPa (0.578 moles), Water Injected @ Bottom] 


Figure F11 - Production History of Run VLC10. 
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NOTE: Average Run Conditions: Continuous CO, Injection at Top, 1.0 MPa, 21°C 
Model Parameters: Average Injection Rate = 308 cc/hr, U, = 1053.0 mPa.s, 


@ = 35.56%, k = 10.75 darcies, S, = 95.88%, Sy. = 4.12% 
Model Type: Linear 


[0.935 HCPV CO, @ 1.0 MPa (0.448 moles), Water Injected @ Bottom] 


Figure F12 - Production History of Run VLC11. 
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NOTE: Quarter of A 5-Spot 
Model Parameters: Average Injection Rate =92.42 cc/hr, 1, = 603 mPa.s, 


g = 40.1 %, k = 12.0 darcies, Sp; = 93.3 %, Sy. = 6.7 % 


~s 


3.0 


40 


WwW 
i=) 


3 


' sm?/sm 


No 
So 


Water-Oil Ratio 


— 
i=) 


S S So 
BS n [o.) 


oS 
bo 
Oil Produced-Fluid Injected Ratio, sm?/sm° 


0.0 


[0.20 HCPV CO, @ 2.5 MPa & 21°C (0.450 mol), 4:1 WAG, 10 Slugs] 


Figure F13 - Production History of Run H2D1. 
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NOTE: Quarter of A 5-Spot 
Model Parameters: Average Injection Rate = 1.55 m/d, , = 603.0 mPa.s, 


9 = 39.7 %, k = 14.5 darcies, Sq; = 90.7 %, Swe = 9.3 % 


[0.20 HCPV CO, @ 2.5 MPa & 21°C (0.423 mol), 4:1 WAG, 10 Slugs] 
Figure F14 - Production History of Run H2D2. 
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NOTE: Quarter of A 5-Spot 
Model Parameters: Average Injection Rate = 2.54 m/d, , = 603.0 mPa.s, 


9 = 39.8%, k = 13.2 darcies, So; = 86.9 %, Swe = 13.1% 


[0.20 HCPV CO, @ 2.5 MPa & 21°C (0.406 mol), 4:1 WAG, 10 Slugs] 
Figure F15 - Production History of Run H2D3. 
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[0.20 HCPV CO, @ 2.5 MPa & 21°C (0.408 moles), 4:1 WAG, 10 Slugs] 


Figure F16 - Production History of Run H2D4. 
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NOTE: Quarter of A 5-Spot 
Model Parameters: Average Injection Rate = 451.23 cc/hr, L, = 603 mPa.s, 


9 = 42.8 %, k = 12.3 darcies, So; = 81.2 %, Swe = 18.8 % 


[0.20 HCPV CO, @ 2.5 MPa & 21°C (0.407 mol), 4:1 WAG, 10 Slugs] 
Figure F17 - Production History of Run H2D5. 
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NOTE: Quarter of A 5-Spot 
Model Parameters: Average Injection Velocity = 0.78 m/d, 1, = 1058 mPa.s, 


= 37.4 %, k = 12.18 darcies, S,; = 88.1 %, Sy, = 11.9 % 


[0.20 HCPV CO, @ 2.5 MPa & 21°C (0.403 mol), 4:1 WAG, 10 Slugs] 
Figure F18 - Production History of Run H2D6. 
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Model Parameters: Average Injection Velocity = 1.55 m/d, 1, = 1058 mPa.s, 


g = 41.6 %, k = 13.9 darcies, S,j = 88.3 %, Swe = 11.7 % 


[0.20 HCPV CO, @ 2.5 MPa & 21°C (0.445 mol), 4:1 WAG, 10 Slugs] 
Figure F19 - Production History of Run H2D7. 
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Model Parameters: Average Injection Rate = 2.6 m/d, |, = 1058 mPa.s, 


9 = 40.5 %, k = 12.5 darcies, So, = 90.7 %, Swe = 9.3 % 


[0.20 HCPV CO, @ 2.5 MPa & 21°C (0.432 mol), 4:1 WAG, 10 Slugs] 


Figure F20 - Production History of Run H2D8. 
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9 = 38.5 %, k = 14.5 darcies, So, = 89.9 %, Swe = 10.1 % 
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[0.20 HCPV CO, @ 2.5 MPa & 21°C (0.407 mol), 4:1 WAG, 10 Slugs] 


Figure F21 - Production History of Run H2D9. 
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NOTE: Quarter of A 5-Spot 
Model Parameters: Average Flow Velocity =3.81 m/d, 1, = 1058.0 mPa.s, 


9 = 38.1 %, k = 12.3 darcies, So; = 89.1 %, Swe = 10.9 % 


[0.20 HCPV CO, @ 2.5 MPa & 21°C (0.407 mol), 4:1 WAG, 10 Slugs] 
Figure F22 - Production History of Run H2D10. 
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Model Parameters: Average Injection Rate = 0.78 m/d, H, = 1842 mPa.s, 


9 = 41.3 %, k = 12.7 darcies, Soi = 93.2 %, Swe = 6.8 % 


(0.20 HCPV CO, @ 2.5 MPa & 21°C (0.460 moles), 4:1 WAG, 10 Slugs] 


Figure F23 - Production History of Run H2D11. 
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Model Parameters: Average Injection Rate = 1.55 m/d, H, = 1842 mPa.s, 


¢ = 43.7 %,k = 14.1 darcies, S,; = 91.4 %, Swe = 8.6 % 


[0.20 HCPV CO, @ 2.5 MPa & 21°C (0.462 mol), 4:1 WAG, 10 Slugs] 
Figure F24 - Production History of Run H2D12. 
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Model Parameters: Average Injection Rate = 2.54 m/d, LW, = 1842 mPa.s, 


g@=41.1 %,k =12.5 darcies, S,; = 87.8%, Swe = 12:2:% 


[0.20 HCPV CO, @ 2.5 MPa & 21°C (0.419 mol), 4:1 WAG, 10 Slugs] 


Figure F25 - Production History of Run H2D13. 
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Model Parameters: Average Flow Velocity = 3.17 m/d, L, = 1842 mPa.s, 


9 = 39.5 %, k = 13.4 darcies, So, = 90.1 %, Swe = 9.9 % 


[0.20 HCPV CO, @ 2.5 MPa & 21°C (0.418 mol), 4:1 WAG, 10 Slugs] 


Figure F26 - Production History of Run H2D14. 
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NOTE: Quarter of A 5-Spot 
Model Parameters: Average Flow Velocity = 3.81 m/d, H = 1842 mPa.s, 


g = 39.8 %, k = 13.1 darcies, So; = 90.3 %, Swe = 9.7 % 


[0.20 HCPV CO, @ 2.5 MPa & 21°C (0.421 mol), 4:1 WAG, 10 Slugs] 
Figure F27 - Production History of Run H2D15. 
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NOTE: Quarter of A 5-Spot 
Model Parameters: Average Flow Velocity = 0.78 cc/hr, UL, = 3295 mPa.s, 


g = 39.9 %, k = 12.7 darcies, So; = 90.0 %, Swe = 10.0 % 


[0.20 HCPV CO, @ 2.5 MPa & 21°C (0.421 mol), 4:1 WAG, 10 Slugs] 
Figure F28 - Production History of Run H2D16. 
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Figure F29 - Production History of Run H2D17. 
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g = 39.1 Jo, k = 12.9 darcies, So; = 89.5 %, Swe = 10.5 % 


[0.20 HCPV CO, @ 2.5 MPa & 21°C (0.411 mol), 4:1 WAG, 10 Slugs] 
Figure F30 - Production History of Run H2D18. 
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Model Parameters: Average Flow Velocity = 3.17 m/d, 1, = 3295 mPa.s, 


g = 38.5 %, k = 13.6 darcies, So, = 93.1 %, Swe = 6.9 % 
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[0.20 HCPV CO, @ 2.5 MPa & 21°C (0.405 mol), 4:1 WAG, 10 Slugs] 


Figure F31 - Production History of Run H2D19. 
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NOTE: Quarter of A 5-Spot 
Nel Parameters: Average Injection Rate = 3.81 m/d, Ly = 3295 mPa.s, 


9 = 38.1 %, k = 13.6 darcies, So; = 94.0 %, Swe = 6.0 % 


[0.20 HCPV CO, @ 2.5 MPa & 21°C (0.420 mol), 4:1 WAG, 10 Slugs] 
Figure F32 - Production History of Run H2D20. 
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Model Parameters: 


Average Injection Rate = 0.78 m/d, [, = 3607 mPa.s, 


g = 44.3 %, k = 12.6 darcies, Soi = 92.1 ToS wee e 


[0.20 HCPV CO, @ 2.5 MPa & 21°C (0.479 mol), 4:1 WAG, 10 Slugs] 
Figure F33 - Production History of Run H2D21. 
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Model Parameters: Average Injection Rate = 1.55 cc/hr, 1, = 3607 mPa.s, 


g = 39.5 %, k = 12.4 darcies, So; = 91.5 %, Swe = 8.5 % 


[0.20 HCPV CO, @ 2.5 MPa & 21°C (0.424 mol), 4:1 WAG, 10 Slugs] 


Figure F34 - Production History of Run H2D22. 
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Model Parameters: Average Injection Rate = 2.54 m/d, 1, = 3607 mPa.s, 


9 = 39.7 %, k = 12.4 darcies, So, = 89.0 %, Sw = 11.0 % 


[0.20 HCPV CO, @ 2.5 MPa & 21°C (0.415 mol), 4:1 WAG, 10 Slugs] 


Figure F35 - Production History of Run H2D23. 
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Model Parameters: Average Flow Velocity = 3.17 m/d, H, = 3607 mPa.s, 


9 = 40.2 %, k = 12.8 darcies, S,; = 90.3 %, Swe = 9.7 % 


(0.20 HCPV CO, @ 2.5 MPa & 21°C (0.426 mol), 4:1 WAG, 10 Slugs] 


Figure F36 - Production History of Run H2D24. 
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Model Parameters: Average Flow Velocity = 3.81m/d, LL, = 3607 mPa.s, 
g = 41.0 %, k = 12.9 darcies, S,; = 90.6 %, Sy, = 9.4 % 


[0.20 HCPV CO, @ 2.5 MPa & 21°C (0.436 mol), 4:1 WAG, 10 Slugs] 


Figure F37 - Production History of Run H2D25. 
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NOTE: Quarter of A 5-Spot 
Model Parameters: Average Injection Rate = 308 cc/hr, 1 = 1058 mPa.s, 
g = 36.1 %, k = 22.0 darcies, Soi = 92.3 %, Swe = 7-7 % 
CO, Injected at 1/10 of Water Rate 
[0.20 HCPV CO, @ 1.0 MPa & 21°C (1.12 mol), 4:1 WAG, 10 Slugs] 
Figure F38 - Production History of Run H2D26. 


368 


[ao es 
t at ee} 
wd 
a. - 
w " » 
| * 
, @ . 
a) 
ra r 
: x 
_ =. * '\ : 
~~ an qo-< alitad se Peakypunitien 0.0) 
va 
“ 
i 
Yio t ro 
! 
} 
/ Ii} 
. : 
* f r a 
a j 
s > WAR 
a Hi | ‘ 
oh ee 
Nel Es 
il i [ 
bot 
t = | ; % t 
-_ m4} L~ ; 7 4 
ery \ 
i lS Sys % 
e Oa at 0.0 


1 el borwial teott | . 
joq? 7AM wed SAM 
~ Bl « 2) aise Fle abd) enti qa eet hom 
2 OS ok es | 2 aiid OSS = Lee 
do A ree HW OE is tenosipl SOU 
legate St Aw th aboot iTS 2 MAO ® CO VION Of01 : 
QOS nei lo amet notenbend - Si yt’ 


Gas-Oil Ratio, sm?/sm? 


Oil Recovery, %HCPV 


— 
So 


= 
in 


40 


WwW 
i>) 


N 
So 


20 


Legends 
-o- GOR 
-s WOR 


— 
Nn 


3 


. sm?/sm 


jan 
(=) 


Water-Oil Ratio 


Ns 


/sm 


Ww 
8 
S 
S 
es 
i) 
oo 
3 


ee 0.6 


0.4 


0.2 


Oil Produced-Fluid Injected Ratio, sm 


0.0 
0.0 0.5 1.0 tS 2.0 i) 
Fluid Injected, Cum. PV 


NOTE: Quarter of A 5-Spot 
Model Parameters: Average Injection Rate = 308 cc/hr, L, = 1055 mPa.s, 
g = 40.0 %, k = 14.1 darcies, S,; = 91.1 %, Sy, = 8.9 % 
CO, Injected at 1/5 of Water Rate 
[0.20 HCPV CO, @ 1.0 MPa & 21°C (0.154 mol), 4:1 WAG, 10 Slugs] 
Figure F39 - Production History of Run H2D27. 
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NOTE: Quarter of A 5-Spot 
Model Parameters: Average Injection Rate = 308 cc/hr, , = 1055 mPa.s, 
g = 42.1 %,k = 14.5 darcies, So; = 84.9 %, Sy, = 15.1 % 
CO, Injected at Water Rate 
[0.20 HCPV CO, @ 1.0 MPa & 21° C (0.178 mol), 4:1 WAG, 10 Slugs] 
Figure F40 - Production History of Run H2D28. 
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NOTE: Quarter of A 5-Spot 
Model Parameters: Average Injection Rate = 308 cc/hr, U,, = 1055 mPa.s, 
g = 39.5 %, k = 16.6 darcies, S,; = 90.3 %, Sy, = 9.7 % 
CO, Injected at 2 times of Water Rate 
[0.20 HCPV CO, @ 1.0 MPa & 21°C (0.167 mol), 4:1 WAG, 10 Slugs] 
Figure F41 - Production History of Run H2D29. 


Water-Oil Ratio, sm?/sm 


3 


Oil Produced-Fluid Injected Ratio, sm /sm? 


Dae 


7 - 
7 
ee] 
4 
: i. . : 
“| on” \ 
' f " 
j | 
J 5 FI \ 
' & VA 
f- % 
’ 
Ww 
; leas 
‘ 
| ; 
~ 
j jf 
| 
Yo 
i / 4 “ r 
| ’ 4 «@ 
a ‘ 
d ' , a ve 
_ eye ees 0 eg v 
; $9 
4 
ry 
a > 


“a ° id 
al a 
an are 


se : Se eng ey nye 
¢ af oi 20 - aR 
vq it et, bth] 


8g8-2.A torte TOM 

ha BOT = iI qivan OOF = ote npitoaysl ogarerA jeroteriiawnT 

RIL = yh PCM 2S geisha = 4 PURE we 
ate 9 -ateW Yo sagan S30 batae{al OD 


Aguila 4M Mom TOE > is Titantron! Mi mein 


_— 7 7 — 


Gas-Oil Ratio, sm?/sm? 


Oil Recovery, ZHCPV 


20 


15 Legends 
-o- GOR 
-s WOR 
10 
5 


40 O 0.8 


WwW 
Oo 
t) 
e 


0.6 


0.4 


0.2 


0.0 
0.0 0.5 1.0 iS 2.0 Pie) 
Fluid Injected, Cum. PV 


NOTE: Quarter of A 5-Spot 
Model Parameters: Average Flow Velocity = 0.831 m/d, ., = 1058 mPa.s, 
g = 39.5 %, k = 16.6 darcies, S,; = 91.0 %, Sy, = 9.0 % 
CO, Injected at 5 times of Water Rate 
[0.20 HCPV CO, @ 1.0 MPa & 21°C (0.087 moles), 4:1 WAG, 10 Slugs] 
Figure F42 - Production History of Run H2D30. 
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NOTE: Quarter of A 5-Spot 
Model Parameters: Average Flow Velocity = 1.29 m/d, L, = 1055 mPa.s, 
g= 41.4 %, k = 13.3 darcies, So; = 86.4 %, Sy_ = 13.6 % 
CO, Injected at Water Rate 
[0.20 HCPV CO, @ 1.0 MPa & 21°C (0.160 mol), 4:1 WAG, 10 Slugs] 
Figure F43 - Production History of Run H2D31. 
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NOTE: Quarter of A 5-Spot 


Model Parameters: Average Flow Velocity = 0.831 m/d, H, = 1055 mPa.s, 


g = 40.6 %, k = 13.3 darcies, S,; = 91.3 %, Sy, = 8.7 % 

CO, Injected at Water Rate 

[0.20 HCPV CO, @ 1.0 MPa & 21°C (0.171 mol), 4:1 WAG, 10 Slugs] 
Figure F44 - Production History of Run H2D32. 
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Model Parameters: Average Flow Velocity = 2.08 m/d, 1, 
g = 38.5 %, k = 12.5 darcies, S,; = 89.2 %, Sy, = 10.8 % 
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[0.20 HCPV CO, @ 1.0 MPa & 21°C (0.144 mol), 4:1 WAG, 10 Slugs] 


Figure F45 - Production History of Run H2D33. 
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NOTE: Quarter of A 5-Spot 
Model Parameters: Average Flow Velocity = 2.54 m/d, 1, = 1058 mPa.s, 


g = 39.2 %, k = 14.1 darcies, S,; = 90.6 Jo, Swo = 9-4 % 


[0.05 HCPV CO, @ 2.5 MPa (0.104 moles), Soak time = 0 days] 
Figure F46 - Production History of Run H2D34. 
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NOTE: Quarter of A 5-Spot 
Model Parameters: Average Flow Velocity = 3.17 m/d, 1, = 1058 mPa.s, 


g = 38.5 %, k = 14.0 darcies, So; = 91.6 %, Swe = 8.4 % 


0.0 0.5 


[0.05 HCPV CO, @ 2.5 MPa (0.103 moles), Soak time = 3 days] 
Figure F47 - Production History of Run H2D35. 
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Model Parameters: Average Flow Velocity = 2.54 m/d, Hy = 1058 mPa.s, 


§ = 38.6 %, k = 10.7 darcies, Soi = 91.7 %, Swe = 8.3 % 


[0.05 HCPV CO, @ 2.5 MPa (0.104 moles), Soak time = 4.83 days] 
Figure F48 - Production History of Run H2D36. 
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Model Parameters: Average Flow Velocity = 2.54 m/d, u, = 1058 mPa.s, 


9 = 39.4 %, k = 11.3 darcies, S,; = 89.1 %, Sw. = 10.9 % 


[0.05 HCPV CO, @ 2.5 MPa (0.104 moles), Soak time = 10 days] 
Figure F49 - Production History of Run H2D37. 
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Model Parameters: Average Injection Rate =308.0 cc/hr, HW, = 1058. mPa.s, 


9 = 37.7 %, k = 4.4 darcies, S,; = 87.1 %, Swe = 12.9 % 


[0.20 HCPV CO, @ 2.5 MPa & 37°C (0.354 mol), 4:1 WAG, 10 Slugs] 
Figure F50 - Production History of Run 121939. 
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NOTE: Quarter of A 5-Spot 


Model Parameters: Average Injection Rate =308.0 cc/hr, i, = 5200 mPa.s, 


g = 37.6 %, k = 3.9 darcies, S,; = 86.9 %, S,,, = 13.1% 


[0.10 HCPV CO, @ 2.5 MPa & 37°C (0.177 mol), 4:1 WAG, 10 Slugs] 
Figure F51 - Production History of Run H2D40. 
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NOTE: Quarter of A 5-Spot 
Model Parameters: Average Injection Rate =308.0 cc/hr, , = 1058 mPa.s, 


g = 37.6%, k = 3.2 darcies, S,; = 86.9%, Sy, = 13.1% 


[0.20 HCPV CO, @ 3.14 MPa & 37°C (0.354 mol), 4:1 WAG, 10 Slugs] 
Figure F52 - Production History of Run H2D41. 
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Model Parameters: Average Injection Rate =308.0 cc/hr, 1, = 603 mPa.s, 


9 = 37.3%, k = 6.2 darcies, So; = 86.8%, Swe = 13.2% 


[0.20 HCPV CO, @ 4.8 MPa & 37°C (0.375 mol), 2:1 WAG, 10 Slugs] 


Figure F53a - Production History of Run H2D42a. 
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Model Parameters: Average Injection Rate =308.0 cc/hr, 1, = 603 mPa.s, 


9 = 37.3 %, k = 6.2 darcies, S,; = 40.5 %, Swe = 59.5 % 


[0.20 HCPV CO, @ 4.8 MPa & 37°C (0.370 moles), 2:1 WAG, 10 Slugs] 


Figure F53b - Production History of Run H2D42b. 
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Model Parameters: Average Injection Rate =308.0 cc/hr, U, = 280 mPa.s, 


¢ = 37.0 %, k = 11.9 darcies, S,; = 87.3 %, S,,, = 12.7 % 


[Waterflood @ 3.58 MPa & 21°C] 
Figure F54a - Production History of Run H2D43a. 
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9 = 37.0%, k = 11.9 darcies, So; = 40.8%, Sw. = 59.2% 


[0.20 HCPV CO, @ 3.58 MPa & 21°C (0.279 mol), 2:1 WAG, 10 Slugs] 


Figure F54b - Production History of Run H2D43b. 


386 


Be) 


res cjote tae; 


6 ain 


hag 1+ ae tc 


- a 


I a 
babel stor 


i Eat Ss 
f.8e-_ yo ' Ae _ 
eae oe. 
as : 
ae , 
e : 
— snail 


Water-Oil Ratio, sm?/sm? 


Oil Recovery, ZHCPV 


30 
20 
10 
(ee fe 
0.0 0.5 1.0 15 2.0 IRS 3.0 35 
50 0.4 
40 ‘ Recovery 
0.3 
30 
0.2 
20 k 
0.1 
10 
0< 0.0 
0.0 0.5 1.0 i 2.0 2.5 3.0 ey 


Fluid Injected, Cum. PV 


NOTE: Quarter of A 9-Spot 
Model Parameters: Average Injection Rate =308.0 cc/hr, L, = 282 mPa.s, 


g = 37.3 %, k = 8.1 darcies, So; = 81.0 %, Sy. = 19.0 % 


[Waterflood @ 3.58 MPa & 21°C] 
Figure F55a - Production History of Run H2D44a. 
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[0.20 HCPV CO, @ 3.58 MPa & 21°C (0.288 mol), 2:1 WAG, 10 Slugs] 


Figure F55b - Production History of Run H2D44b. 
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APPENDIX G 


Viscosity-Temperature Relationship for Different Oils. 
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Figure G1 - Viscosity-Temperature Relationship for Battrum South Oil No.1. 
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Figure G2 - Viscosity-Temperature Relationship for Battrum South Oil No.2. 
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Figure G3 - Viscosity-Temperature Relationship for Kla-Da-Ing Oil. 
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Figure G4 - Viscosity-Temperature Relationship for Della-Bell Oil. 
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Figure G5 - Viscosity-Temperature Relationship for Epping Oil. 
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Figure G6 - Viscosity-Temperature Relationship for Senlac Oil. 
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Figure G7 - Viscosity-Temperature Relationship for South Aberfeldy 
Oil No.1. 
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Figure G8 - Viscosity-Temperature Relationship for South Aberfeldy 
Oil No.2. 
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